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FOREWORD 


An important aspect of estimating the aerodynamic characteristics of 
V/STOL aircraft in transition between hover and wingborne flight is an under- 
standing of the flow field induced by the efflux of the lifting-propulsive 
device. The influence of this efflux dominates the aerodynamic interference 
effects for aircraft that use high-disk- loading devices such as turbofan or 
turbojet lift engines. At the present time, approximately a dozen organiza- 
tions are working independently on this problem. Several fundamentally dif- 
ferent approaches are being applied toward understanding and solving the 
problem. Recognizing this activity, NASA decided to hold a Symposium on the 
Analysis of a Jet in a Subsonic Crosswind on September 9 and 10, 1969, at the 
Langley Research Center to present the currently available results . 

The objective of this meeting was to provide an opportunity for persons 
working in the area to present and discuss their work and, in many cases, 
their opinions. The presentations were made in two sessions: 

Session I - Experimental Results 

Session II - Empirical and Analytical Descriptions of the Jet and Induced Flow 

This publication includes 15 papers presented at the symposium. The 
attendees are listed at the end. 

It should be emphasized that this publication essentially represents a 
compilation of papers presented by authorities from various governmental and 
private organizations. In order to expedite publication, this document has 
been printed from copy provided by the authors. 

The Langley Research Center wishes to express its appreciation to all 
authors and to the organizations they represent, for their substantial con- 
tributions to this program and for the timeliness of their responses both to 
the program activities and to the symposium. 
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SUMMARY 


A summary indicating the general trends of jet— wake— induced aerodynamic 
effects on V/STOL aircraft in transition flight is given. The jet wake 
including the adjacent flow field is described qualitatively, and a quantita- 
tive measurement of the circulation in the wake is presented. These measure- 
ments indicate that the circulation may be somewhat greater than that generated 
by potential flow around a circular cylinder. Finally, an extensive list of 
related reports has been included. 


INTRODUCTION 


The purpose of this paper is to summarize results of previous experimental 
investigations of the aerodynamic interference effects experienced by V/STOL 
aircraft in the transition speed range. (Additional summary papers are pre- 
sented in refs. 1 to 6.) These effects have been the subject of a large number 
of experimental investigations (refs. 7 to 18). Most of this research effort 
has been directed toward the investigation of the forces and moments induced on 
the aircraft by the interaction of the vertical jets with the free-stream air. 
In addition, a description of the character of the jet wake is presented. 


SYMBOLS 


A 

b 

c 

D 

L 


area, ft 2 (m 2 ) 
wing span, ft (m) 
wing chord, ft (m) 
jet diameter, ft (m) 
lift, lb (N) 
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V /v. 

°°/ J 

X,T,Z 


x,y,z 


increment of lift interference caused by the flow from a jet, lb (N) 

increment of pitching-moment interference caused by the flow from a 
jet, ft-lb (m-N) 

free-stream dynamic pressure, lb/ft^ (N/m 2 ) 
jet dynamic pressure, lb/ft^ (N/m 2 ) 

thrust, lb (N) 

free-stream velocity, fps (mps ) 
jet velocity, fps (mps) 

effective velocity ratio 
coordinate axes 

Cartesian coordinates, ft (m) 

nondimens ional circulation strength, T/2DV 

circulation strength, ft^/sec (m 2 /sec) 

maximum circulation strength for potential flow around half of a 
circular cylinder, SDV^, ft^/sec (m 2 /sec) 

total circulation for a wing with an elliptic lift distribution, 
^L/pV^b, ft 2 /sec (m 2 /sec) 

flap deflection, deg 

jet deflection from free-stream direction, deg 

natural jet coordinate normal to |-t] plane, ft (m) 

natural jet coordinate in the lateral direction, measured parallel 
to Y-axis, ft (m) 

angular cylindrical coordinate, rad 

natural jet coordinate along the axis of the jet, ft (m) 
density, slugs/ft^ (kg/rn.5) 
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AERODYNAMIC INTERFERENCE EFFECTS 


During transition flight of v/STOL aircraft, the jets issuing from the 
aircraft are swept rearward by the free-stream flow and are rapidly rolled up 
into a pair of vortices (fig. l). The wakes induce suction pressures on the 
fuselage and a distribution of downwash over the aircraft. This downwash is in 
effect an induced twist on the wing and tail and an induced camber over the 
length of the airplane . 

The general trend of these jet-induced effects is illustrated in figure 2. 
There is usually a loss in lift which tends to increase with increasing forward 
velocity. The loss in lift is about the same with the tail off the aircraft 
and with the tail on. There is an increment of pitching moment in transition 
flight which tends to increase nose-up with increasing velocity. Because of 
the change in downwash in the vicinity of the tail, there is an additional 
increment of pitching moment induced when the tail is on. 

Recently, an investigation (ref. 18) was undertaken to determine the effect 
of the jet-exit location on the interference on a wing. This investigation used 
a model with a simple fuselage and an unswept, untapered wing with an aspect 
ratio of 6 and a 30-percent-chord slotted Fowler- type flap. This model was 
mounted on a sting-supported strain-gage balance. Two jets, one on each side 
of the fuselage, were mounted independently of the wing at about the 25-percent 
semispan station. The jet exits were positioned at each of the various longi- 
tudinal and vertical locations shown by the plus marks in figure 3- Results of 
this investigation show that negative lift interference was experienced when 
the jet exits were located ahead of the wing. Favorable interference effects, 
however, are encountered with the jets below the wing and aft cf the wing mid- 
chord. The fact that the interference effects are most favorable for locations 
closest to the flap indicates that the jet is probably helping the wing and 
flap achieve their full lift potential. 


VISUALIZATION OF FLOW PHENOMENA 
OF A JET IN A CROSS WIND 


The following description of the flow phenomena involved in the transition 
problem is based primarily on flow visualization. The shape of the jet wake 
when influenced by the free stream is presented in figure 4. These data are 
from an experimental investigation (ref. 19) in which the wake exits from the 
nozzle at a deflection angle to the free stream. (See refs. 20 to 32 for 

other investigations and refs. 33 "to 52 for analyses of the jet in a cross 
wind. ) The free stream deflects the wake back until it tends toward the free- 
stream flow. As a result of this investigation, the empirical equation given 
in figure 4 was obtained. It describes the path of the jet as a function of 
the effective velocity ratio (the square root of the ratio of the free-stream 
dynamic pressure to jet dynamic pressure) and the deflection angle of the jet. 
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A more detailed look at the jet in a cross wind is shown in figure 5 where 
the flows induced in and around the jet are shown. This photograph, which was 
obtained by using water-tunnel flow visualization at the Office National D'Etudes 
et de Recherches Aerospatiale (O.N.E.R.A. ), shows a flat plate with a jet 
exhausting from it. Near the leading edge of the plate are orifices through 
which colored milk is emitted. When the colored filament on the center line 
approaches the jet exit, it divides and is swept around the jet. This indi- 
cates that there is a stagnation point near the front of the jet. The visible 
portion of the jet wake indicates that some of the flow from the free stream 
is sucked into the jet. The filament adjacent to the jet passes beside the jet 
and is then induced upward into the turbulent wake region behind the jet. Even 
the filaments farthest from the center line of the exit are sucked toward the 
jet and into the wake region. In the wake region, a considerable amount of 
entrainment by the jet can be seen. 

In figure 6, two photographs of oil flow on the surface of the plate are 
presented. The lower one is a closeup near the jet exit. A line trails diago- 
nally downstream from the jet and bounds the wake region. A careful inspection 

of the stream lines on the surface shows that the free-stream flow is deflected 
around the periphery of the jet and into the wake boundary. In the wake region, 
the flow near the jet is entrained forward by the jet. Farther downstream of 
the jet exit, the wake flow is deflected away and carried off with the boundary 
of the wake on the surface of the plate. This wake boundary appears to be a 

second pair of vortices which are similar to a Karman vortex street. This vor- 

tex street is probably formed around the potential core of the jet. 

The pressure distribution which is induced on the surface of the plate is 
presented in figure 7* For these particular data, which are from reference 
the effective velocity ratio was 0.25. The region of positive pressure ahead of 
the jet indicates a force augmenting the jet thrust. The region of negative 
pressure adjacent to and aft of the jet indicates a force opposing the jet thrust 
behind the jet. This pressure distribution accounts for the nose-up interference 
pitching moment described earlier. The area covered by the negative pressure 
is noticeably larger than the area covered by the positive pressure. As a 
result there is the lift loss. 

Figure 8 is a photograph of the cross section of the jet wake in a water 
tunnel. The air bubbles are flowing around the jet wake; the white region is 
the jet. The primary features are the swirl of the two elements of vorticity 
in this jet cross section and the fact that the cross section is kidney shaped. 


STRENGTH OF THE VORTICES IN A JET WAKE 


An investigation was conducted recently in the Langley 300-MPH 7- by 
10-foot wind tunnel to determine experimentally the strength of the contra- 
rotating vortices in the wake of a jet in a cross flow by use of a vortex meter. 
The meter consisted of a l/2-inch-diameter paddle wheel which spins in a rota- 
tional flow field. 
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The first phase of the investigation was calibration of revolutions of the 
vortex meter against circulation strength of the flow field. This was accom- 
plished by measuring the strength of the trailing vortex sheet from a 
rectangular-planform wing with an aspect ratio of 6. (See ref. 18 for a descrip- 
tion of the wing.) The vortex meter was placed about l/2 inch aft of the 
trailing edge of the wing at several spanwise locations, and the revolutions 
were measured. In addition, the lift of the wing was measured with a six- 
component strain-gage balance. Then, this measured lift and a calculated lift 
distribution for the wing were used to determine the magnitude and distribution 
of the circulation in the wake of the wing as a function of spanwise location. 

(The measured values are superimposed on the calculated distribution of vortex 
strength in figure 9 where the circulation ratio is presented as a function of 
spanwise location. Inboard of the 30 percent location, the measurements were 
too low to be accurate , and outboard of the 95 percent spanwise location, the 
meter was measuring the change in circulation distribution caused by the rolling 
up of the wing-tip vortex. However, the correlation is good for the portion of 
the span between these two locations, and the measurements over this portion 
were used to establish the vortex-meter calibration. 

Next, the circulation in a jet was measured. For these measurements a 
1- inch-diameter nozzle was mounted through a ground board normal to the free 
stream. The measurements were made in planes normal to the jet axis as shown 
in figure 10. The most complete data were obtained for an effective velocity 
ratio of 0.25 at a distance of 6 jet diameters along the jet path. 

These results, which are presented in figure 11, indicate a nondimens ional 
circulation of approximately 2 for each vortex. The form of the nondimens ional 
circulation is obtained from analyses of the jet which use the circulation 
obtained in potential flow around a circular cylinder (ref. 37 )♦ This is done 
because the cross-sectional shapes predicted by this model are similar to the 
observed shapes. In this model the strength of the vortex sheet describing 
the flow around a circle is 


r(e) = — 2V 00 sin e 

Integrating this circulation around a semicircle gives the maximum strength of 
each vortex as 



The nondimens ional circulation is defined as 



The experimental results indicate that the circulation is approximately twice 
the value predicted by the potential-flow model. 
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It should be indicated that these results are not precise, and the pre- 
liminary nature of the experiment should be emphasized. There was difficulty 
in measuring small amounts of circulation because the meter was calibrated over 
a wide range of rpm. Hence, the sensitivity of the measurement was too coarse 
at low rpm. This led to error in determining the total circulation because the 
regions of low levels of circulation were large in area, and as a result, their 
contribution to the total vortex strength was more than expected. Further, 
since the vortex meter was rigidly mounted on a wand, there was some difficulty 
in alining it to the local flow. It was noted that when the meter was out of 
alinement by about 10° or greater, it tended to oscillate about an equilibrium 
position in the absence of vorticity. This oscillation showed up as revolutions 
on the meter. These results describe only a single test condition and will 
require additional verification before extrapolating to other conditions. 


CONCLUDING REMARKS 


The general trends of the effects of lifting jets on the aerodynamics of 
V/STOL aircraft in the transition speed range have been indicated. The charac- 
ter of the fluid flow in the jet and the adjacent flow field has been described 
qualitatively, and finally, a quantitative measurement of the circulation in the 
wake has been presented. These measurements indicate that the circulation may 
be somewhat greater than that generated by potential flow around a circular 
cylinder. An extensive list of related reports has been included. 
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Figure 1.- Jet wake from a V/STOL aircraft in transition flight rolls up 
into a pair of contrarotating vortices. 
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Figure 2.- General trend of jet-induced lift loss and pitching moment in 

transition flight. 
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Figure 5-- Photograph of the flow induced around and into a jet exhausting 
normal to the free stream. (Photograph from O.N.E.R.A. mot ion -picture 
film no. 575 entitled "Flows With Large Velocity Fluctuations.") 
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Figure 6.- Photograph of oil flow on the surface of a plate which has a 
jet exhausting normal to the plate and to a free stream (O.N.E.R.A. 
photograph). Lower photograph is closeup of upper. 
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Figure 8.- Cross section of the wake of a jet exiting perpendicular to the free 

stream in a water tunnel. This cross section is approximately 6 nozzle diameters 
downstream of the nozzle exit along the jet path. (Photograph from O.N.E.R.A. 
motion-picture film no. 575 entitled "Flows with Large Velocity Fluctuations. ) 
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Figure 10.- Location of normal plane where circulation contours were 

measured. = 0.250. 
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Figure 11.- Jet circulation contours at the |/D = 6 location. The posi- 
tion of the jet center line is indicated. In addition, a schematic 
sketch of the circle swept out by the l/2-inch vortex-meter paddle blades 
is shown. 
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SUMMARY 


A. description of the physical characteristics of the sub- 
sonic jet entering a stream with a non-zero, cross-component is 
presented. Three zones, in which the dynamics of the flow are 
significantly different, are defined. The structure of the 
mean and turbulent motion is discussed with reference to the 
various entrainment processes. 


INTRODUCTION 


Probably the first stimulus for an understanding of the 
mixing which arises in a jet of fluid discharging into a cross- 
flowing stream occurred because of its application to the 
emission of effluents from a chimney stack. Only a cursory ob- 
servation of a plume is needed to show the remarkably rapid 
diffusion which can take place in the natural wind with adiabatic 
or nearly adiabatic thermal lapse rates. It is clear from this 
that the momentum fluxes of the opposing streams create an ef- 
ficient mechanism of entrainment of free-stream fluid. Indeed, 
the principle is so basic that it is used to promote fluid 
mixing in a large number of mechanical systems. 

Lately there has been a resurgence of interest in the de- 
tailed features of the fundamental problem. This is partly be- 
cause of the increasing awareness of the public and government 
about our need to understand the role of turbulent dispersion 
of waste effluent. But there has been stimulation, as well, 
because of the application to jet thrust systems in some air- 
craft. The deflected jets of a V/STOL plane for example, in- 
terfere with the aerodynamic performance. This aspect has re- 
ceived some attention recently by Williams and Wood (ref. 1) 
and Margason (ref. 2). It should perhaps be noted that a much 
earlier paper by Ribner (ref. 3) anticipated this interest with 
an examination of the effect of a jet upon the stability of a 
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jet aircraft. 

Considering the basic problem of a jet in a cross-wind, 
Ruggeri et al (refs. 4,5,6 and 7), in a series of early papers, 
explored the effect of the addition of heat and change of ori- 
fice shape upon the flow. The significance of their work to 
the present discussion is limited, however, since the wind 
tunnel used to contain the cross-stream, was narrow. This es- 
sential two-dimensionality would be expected to limit the 
lateral spread of the jet and, as will be seen later, the effect 
of the flow around the jet edges is important to the entrain- 
ment. The first relevant experimental work on the problem was 
due to Jordinson (ref. 8), who made some simple measurements 
of the trajectory of the round jet. He defined the centreline 
as the point of maximum velocity measured with a pitot-tube 
and pitch meter at successive vertical stations. As well, 
Jordinson surveyed the cross-section of the jet. The result- 
ing contours were very revealing in that they showed the edges 
of the jet to be sheared by the cross-stream and distorted in 
the down-stream direction. The distinctive kidney-shape was 
indicative of a secondary mode of entrainment as will be dis- 
cussed in some detail in subsequent sections. 

A second experimental study by Gordier (ref. 9) carried 
out a series of tests over the same general range of velocity 
ratios as Jordinson but used a jet discharging into a water 
tunnel. The results defining the trajectory were consistent 
with previous work and this was the first indication that a 
Reynolds number similarity existed. The examination of the 
structure of the turbulence within the deflected jet system was 
undertaken by Keffer and Baines (ref. 10). Using hot wire 
anemometer techniques the mean of the turbulent intensities as 
well as the mean velocity quantities were taken. The results 
showed that within the range of velocity ratios used by Jordin- 
son and Gordier, a self-similarity of the mean flow field could 
be defined if one chose to scale the velocity terms by a ref- 
erence velocity difference along a natural system of co-ordi- 
nates intrinsic to the jet flow (see figure 1). The success 
of this scaling was only limited for the mean turbulent inten- 
sities, however. The rather high level of turbulence intensity 
precludes us from considering the data to be more than a quali- 
tative specification of magnitude. 

The problem has been extended further by a number of in- 
vestigators both experimentally and by phenomenological model 
equations. No attempt is made here to be comprehensive but a 
few examples are mentioned to show the present trend of the 
research. Pratte and Baines (ref. 11) examined the far field 
behaviour of the jet, that is, the region well beyond the zone 
of maximum curvature. Somewhat surprisingly, it appears that 
the proper scaling function which one should use to collapse 
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the trajectories on a single functional curve depends upon the 
velocity ratio rather than the momentum ratio. Fan (ref. 12) 
has introduced the effect of a stratified environment upon the 
growth and trajectory of what should be termed a plume. Ana- 
lytically, he has suggested a mathematical model based upon the 
drag of the plume in the cross-stream. Margason (ref. 13), in 
a series of wind tunnel tests considered the effects of large 
variation of the angle of entry of the jet to the main stream, 
as did Platten and Keffer (ref. 14). In the former study em- 
pirical fits were made to the trajectories found by a number of 
workers with rather uniform results. Platten and Keffer pro- 
posed a model equation to account for the two basic modes of 
entrainment which exist in the deflected jet. A solution of 
the equations gave good fits to the experimental trajectories 
over a limited range of injection angles. 

As a result of the work briefly outlined above, a distinct 
picture of the physical nature of the jet is beginning to emerge 
There are some important areas of disagreement, not, happily, 
in the experimental results, but in the mathematical models 
which are used to interpret them. Most significant perhaps, 
is the controversy about the mechanism of entrainment of ambient 
fluid by which the jet grows. This paper will inevitably in- 
corporate a personal bias into the argument and the emphasis 
will be placed upon the structure of the turbulence within the 
flow rather than the phenomenological replacement of the tur- 
bulence. It would appear that while the salient features are 
reasonably well understood, most of the intimate and therefore 
interesting details still command our intensive and enlightened 
interpretation . 


SYMBOLS 

d - orifice diameter 

R - velocity ratio 

U - axial jet velocity 

U c - maximum U at any cross-section 

Uj - initial jet velocity 

U - free stream velocity 

o 

U v - vortex entrainment velocity 

u — turbulent velocity 


x,y,z 


Cartesian co-ordinates 


£,n,C - curvilinear co-ordinates 

£*,x* - distances measured from virtual source 

0 - local jet angle 

0_. - angle of jet injection 

PHYSICAL DESCRIPTION 


In the discussion which follows we reduce the problem to 
its simplest form. This becomes an examination of a round jet 
of fluid discharging from an infinite plane into an ambient, 
miscible fluid of the same density, which has a non-zero, cross- 
stream component of velocity. The influence of the surrounding 
geometry, such as the presence of a chimney stack or the fuse- 
lage of an aircraft, will be ignored. The significant para- 
meter will thus be the velocity ratio R defined as the exit 
jet velocity U-? divided by the constant and uniform cross- 
stream velocity U Q . It will further be assumed that the 
initial absolute velocity in the jet is large enough that the 
developed flow, which takes place upon the mixing of the streams , 
will be turbulent. A second parameter is introduced, the angle 
of injection 0j . Clearly, a complete variation of 0j through 
180° would include the special cases of a jet in an opposing 
stream and the jet in a co-flowing stream. Generally, however, 
we shall restrict our interest to moderate deviations of 0j 
from the 90° direction. A schematic of the jet is shown in 
figure 1. 

It is convenient to think of the jet as arbitrarily sepa- 
rated into three regions: 

(1) the source flow, where the effect of orifice geometry, 
plate and/or tunnel boundary layer and potential core may be 
important , 

(2) the curvilinear zone, where the jet flow is fully de- 
veloped from the standpoint of the turbulence structure and the 
entrainment mechanisms proceed in a relatively straightforward 
manner and 

(3) the far region, where the jet flow, having approached 
its asymptotic state, is almost entirely dissipated and, on a 
mean flow basis, indistinguishable from the main stream. 


The Source Flow 

Relatively little work has been done in this region. Es- 
sentially the flow at exit bears a close resemblence to the 
free jet in that it possesses a core of potential fluid. The . 
term potential is relative. The fluid emerging from the orifice 
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is considered to have a negligible amount of turbulence compared 
to that which subsequently develops in the fully mixed region of 
the jet. It is probable that the flow in this zone is Reynolds 
number dependent, as was found by Baines (ref. 15) for the free 
jet, although not enough experimental data are available to make 
a definitive statement. It has been shown (refs. 10, 11, and 14) 
that the length of the potential core varies directly with R, 
but again the specific functional trend has not been reported 
in the literature for a large range of velocity ratios. 

Keffer and Baines (ref. 10) indicated that if substitution 
is made for the core region and a virtual origin or point source 
defined by extrapolating the developed jet flow back through the 
plane of emission, the decay of centreline velocity will plot 
universally on a single functional plot. This implies a kine- 
matic or Reynolds number similarity for all deflected jets, at 
least for the range of flows tested. At very low velocity ra- 
tios, however, of the order of 4 or less, the flow loses this 
characteristic. A premature deflection of the core region by 
the pressure field occurs and the trajectory centreline, extra- 
polated back to the virtual origin, does not go through the 
physical origin for the flow. For a jet discharging at right 
angles to the stream this means that the trajectory at the 
source is not tangent to the x-axis. This can lead to diffi- 
culties in the mathematical modelling of the flow. 

In addition to the displacement or shift effect it was 
found that, at these very low velocity ratios, the jet momentum 
relative to the free stream is not sufficient to allow the flow 
to escape the boundary layer of the exit plane. This introduces 
a complication for, unless particular care is taken to eliminate 
the boundary layer, the results from different laboratories will 
be influenced by the different magnitude of boundary layer en- 
trainment. Since this introduces an additional and relatively 
uncontrolled parameter into the system and we will not consider 
these low values of R further. It will be assumed that the 
velocity ratio is sufficiently high that the jet flow can escape 
the influence region of the surface. It should perhaps be noted 
that a related problem is of some significance. The performance 
of the jet when deflected by a gradient wind is certainly rele- 
vant. For reasons of simplicity, however, most studies to date 
have been concerned with the uniform wind profile. 

The previously mentioned experiments of Ruggeri et al con- 
sidered the effect of shape change of the discharge opening. 

This can be expected to affect the nature of the source flow. 

Its significance must be negligibly small in the developed re- 
gion of the flow, however, where it is possible to replace the 
initial conditions by a virtual origin. For an orifice flow, 
a contraction co-efficient is generally necessary to reduce the 
data to a convenient base, compatible say with a sharply con- 
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verging nozzle. If a length of tube is used to eject the jet 
fluid, the exit profile will be markedly non-uniform, the exact 
shape depending upon the length of tube or the stage of develop- 
ment of the tube flow. It will thus be necessary to define an 
appropriate momentum flux for the jet by means of a correction 
co-efficient multiplying the square of the average velocity. 


The Curvilinear Region 

As the jet issues from the exit plane the relative veloci- 
ty between the jet and stream fluid creates a turbulent mixing 
layer around the periphery. This mixing layer grows in the di- 
rection of the jet flow, eroding the potential zone of fluid in 
a manner similar to the free jet. In the present case, the 
growth is more rapid and the potential core length much shorter. 
For velocity ratios of 4 , 6 and 8 the core length is of the 
order 1, 2 and 3 orifice diameters respectively. For the free 
jet, the core length is about 6. At the end of the core, the 
jet fluid is fully turbulent and the lateral velocity profile, 
taken typically in a streamwise traverse, is more or less 
Gaussian. We define this as the beginning of the curvilinear 
zone and the flow is said to be fully developed at this point. 

At present, the behaviour of the flow in the lee of the 
jet, conveniently termed the 'wake', is not well understood. 

It is possible that the analogy of Jordinson has some validity. 
He suggests that the jet can be considered to act as a porous 
cylinder with suction, the sink effect being necessary to ac- 
count for the entrainment of free-stream fluid. Although of- 
fering some convenience this model cannot be expected to tell 
us much about the true nature of the entrainment processes. 
Clearly the mechanisms must be different in the stagnation re- 
gion at the front of the jet and in the 'wake* . 

The cross-sectional profiles of Jordinson and Keffer and 
Baines (for example see figure 2) have shown the strong lateral 
deflection of the jet sides by the shearing of the cross-stream. 
Turbulent entrainment of free stream fluid will occur at the 
sides because of the relative motion. As it is swept around 
into the wake region, the fluid will be incorporated into the 
main jet system. The result is a helical circulation pattern. 
The generation of the circulation will continue as long as the 
jet has a component of velocity normal to the free stream. 
Appropriately this helicity has been interpreted as a counter- 
rotating, vortex pair and it can be shown from analytical argu- 
ments (ref. 14) that the vorticity generation will reach a max- 
imum within this curvilinear region. A plot of the rate of 
change of circulation velocity is shown in figure 3 as a func- 
tion of deflection angle. It might be suspected that the 
helicity would play a major role in the entrainment of free- 
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stream fluid and the analytical model of Platten and Keffer to- 
gether with the experimental results, appears to bear this out. 

In the frontal regions of the jet, the flow is partially 
stagnated. As the jet deflects, there will be a component of 
motion of the free— stream fluid parallel to the mean direction 
of the jet flow. This situation is just that required for a 
free jet type of entrainment. Although the details of how a 
turbulent front expands into a non-turbulent region of fluid 
are not at all well understood at present, heuristically we _may 
argue in terms of a phenomenological model. A requirement for 
entrainment is simply the existence of relatiye motion between 
streams of fluid. An instability develops which breaks up into 
turbulent mixing and the high velocity fluid, or jet flow, dif- 
fuses into or entrains the outer stream. We would thus use, 
in lieu of the precise location of the turbulent -non-turbulent 
interface, a simple specification of the jet boundary. Con- 
veniently this could be the point where the jet velocity excess 
has decreased to an arbitrary level, say 10%. 

From a physical point of view we thus have two mechanisms 
of entrainment, the ordinary free-jet flow type which can be 
handled by the conventional phenomenological approach and the 
entrainment supplied by the helicity in the 'wake region of 
the jet. Clearly a successful analytical model must in some 
way incorporate these two processes. 

The developed jet flow in this zone displays some of the 
characteristics associated with the simple free-jet. . The 
measurements of Keffer and Baines of the mean flow distribution 
show a self-similarity of the axial velocity vector when taken 
along a curvilinear lateral co-ordinate, n. The scaling was 
achieved by using the excess velocity, (U-U 0 ) and dividing by 
the local centreline velocity excess, (U c -Uo) • The collapsed 
data shown in figure 4, resembles the conventional Gaussian 
distribution plot although the scatter is marked. The func- 
tional shape is not really surprising. With this technique 
of stretching the data small differences tend to be obscured. 
Nevertheless, the results show no unexpected anomalies. 

The variation of centreline velocity ,. i . e . along the jet 
trajectory, shows a universality which again tends to. resemble 
the free jet (figure 5). The decay of mean velocity is more 
intense, however, as would be expected with the higher rates 
of entrainment and the variation of mean velocity at the centre- 
line of the jet does not have the simple inverse linear rela- 
tionship with distance characteristic of the free jet. It is 
significant that when the differences of the source as a func- 
tion of velocity ratio are accounted for by allowing the posi-. 
tion of the virtual origin to vary, the. mean flow characteristics 
are strongly self-similar for this curvilinear region. 
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The turbulence measurements display less self-similarity 
when the same method of plotting is used (figure 6). As with 
all simple free turbulent shear flows , the structure of the 
turbulence requires a longer period than the mean flow to reach 
an equilibrium state. In the present situation, the intense 
mixing and decay prevents this and an equilibrium situation is 
not reached. In fact, since the region is really an interim 
state for the whole of the deflected jet flow, it is probably 
incorrect to speak in terms of an equilibrium or asymptotic 
state for the curvilinear zone. 

From the above remarks we can summarize the main features 
of the flow in this second region. The curvilinear zone clearly 
has most of the characteristics of simple jet flows, for example, 
entrainment, axial velocity decay, self-similarity and turbulence 
structure. We can choose to interpret the flow as a special case 
of the simple jet if an appropriate system of curvilinear co- 
ordinates is defined. The most distinguishing feature is the 
presence of the ’wake* region and the attendent circulatory flow 
which exists within it. Although it is possible to exclude this 
'wake' from the basic jet flow analysis, the consequences of the 
extra entrainment cannot be ignored. 


Far Field Region 

Beyond the zone of maximum curvature, entrainment has re- 
duced the jet velocity excess to a point where it cannot be 
measured by conventional techniques and it becomes necessary to 
define the jet boundary for these conditions by the use of a 
tracer such as low level heat or smoke. We may think of this as 
the onset of the final region of the deflected jet flow. Pratte 
and Baines (ref. 11) were able to observe the far field behaviour 
visually and it showed a double-valued power law variation for 
the centreline trajectory when scaled in terms of the velocity 
ratio. Their functions took the form, 


— * (X) n 

dR MR ; 


where the exponent ' n* was 1.0 up to an £/dR of 2 and 0.33 be- 
yond that. Their range of velocity ratios extended from 5 to 
35. The earlier results of Keffer and Baines, which accounted 
for the variation of the virtual origin, suggested the relation 


dR 2 


n. 


"(€*) 

dR 2 
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with ’n* = 1.0 in the first portion of the flow but falling off 
as the limit of observation (£*/dR 2 = 2.0) was reached. These 
results are consistent in that they imply that the zone of 
curvilinear flow is significantly different than the far field. 

Although the axial component of mean velocity is virtually 
zero, this region of the flow should still show the effect of 
the swirl generated in the curvilinear zone since vorticity is 
remarkably long-lived. The jet in this final phase contains its 
original vertical momentum. It will therefore continue to rise, 
but at an ever-decreasing rate as the momentum is diffused over 
the steadily increasing area of each horizontal plane. Because 
of this, additional vorticity will continue to be generated all 
along the trajectory of the jet. But the rate of generation 
will decrease as we travel further along the jet path and the 
processes of ordinary viscous dissipation will cause the abso- 
lute magnitude of the vorticity to decrease. The limit state 
of this flow will thus be a pair of weak, counter-rotating line 
vortices, being carried along by the main stream. 

The structure of the turbulence in this final phase must 
also be dominated by the viscous decay. The smaller eddies will 
be dissipated most rapidly and as is found with grid-produced 
turbulence, the largest eddies, which are characteristic of the 
generating grid, remain. They tend to retain their identity 
while the smaller scales disappear. In a sense the counter- 
rotating vortex pair in this present system is analogous to the 
large eddy of the grid flow. It is probably more accurate to 
view these vortices, at least in the final phase, as general 
patterns of circulation rather than discrete and perhaps 
measureable line vortices. 

The entrainment mechanism in this final phase must depend 
almost entirely upon the vortex flow since there no longer is 
an effective component of relative axial motion between the 
jet and the free stream. We would expect, therefore that the 
absolute entrainment would decrease and the rate of jet spread 
would fall off. The measurements of Pratte and Baines support 
this. A plot of the width of the jet with respect to the axial 
path gives 


M * ( J-) n 

dR dR 


where * n* is 1.35 in the curvilinear region, 
ruptly to 0.33 for the far field. 


This changes ab- 
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DISCUSSION 


For meteorological applications, the far field region 
would be of most importance. Interest for aerodynamic study, 
however, would centre upon the first two regions of the flow. 
Within the zone of maximum curvature, where the flow has reached 
a developed state of turbulence and where the dynamic effects 
are strong , attempts have been made to formulate model equations 
for the flow. For example, in the work of Fan (ref. 12), the 
presence of entrainment, additional to that supplied by free 
jet theory, is attributed to a pressure difference around the 
jet flow. By employing a suitable drag co-efficient, Fan's 
model appears to predict the trajectories of both buoyant and 
non-buoyant plumes reasonably well. On physical grounds, how- 
ever, this concept is not attractive, since in regions away 
from the source, the processes involved in jet entrainment are 
most likely to result from free stream effects not depending 
upon pressure gradients. 

In contrast, Platten and Keffer (ref. 14) have suggested 
that the extra entrainment must result from the observed circu- 
latory flow, i.e. the so-called vortex pair. As discussed above, 
it may be more reasonable to think of these motions as large 
turbulent eddies within the jet system, rather than a pair of 
strong line vortices such as would be observed from a delta 
wing. This would not affect the analysis, however. The work 
of Pratte and Keffer (ref. 16) on turbulent swirling jets has 
verified that a marked increase in entrainment results when 
swirl is added to a simple jet flow, the spread rate increasing 
by a factor of two for a swirl velocity of the same order 
initially as the axial velocity. Although it is not yet clear 
whether the free jet and vortex contributions to entrainment 
are interactive, the simple additive model of Platten and 
Keffer, predicts the experimental results well. 

Aside from obvious practical applications , the study of 
the deflected jet is of considerable fundamental interest in 
that it represents a particular example of a free turbulent 
shear flow. As an exercise in basic research, the problem may 
be expected to tell us something about the structure of turbu- 
lence generally. Although the simplest of free turbulent shear 
flows (jets and mixing layers) are self-preserving, (ref. 17) 
enabling spread and decay rates to be predicted from the equations 
of motion, the present flow is not. The presence of the co- . 
flowing component of the external stream precludes this particu- 
lar analysis even when use is made of the intrinsic co-ordinate 
system. As mentioned above, however, the results show the flow 
to be roughly self-similar in the zone of curvature. This does 
enable some simplification to be made in the analysis of the . 
problem. Much more experimental work is required in this region 
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to establish the more subtle characteristics of the turbulence, 
however. 
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Figure 2.- Typical velocity contours. 












THE AERODYNAMICS OF THE LIFTING JET IN A CROSS FLOWING STREAM 
By J. E. Hackett and H. R. Miller 
Lockheed-Georgia Aerospace Research Laboratory 


For many methods of treating the problem of the lifting jet in a cross 
flowing stream theoretically, it is necessary to know, in advance, the mean path 
of the jet, which is usually determined experimentally. For the general design 
proDiem, and particularly where novel configurations are being investigated, 
this can prove awkward. The long-term aim of the work currently in progress at 
the Lockheed-Georgia Research Laboratory is therefore to predict the path of the 
lifting jet plume or plumes with sufficient details of the flow structure to 
allow pressures to be calculated on nearby surfaces. 
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Like many other investigators in this area, we started our theoretical 
approach by repeating H. C. Chang's (nee Lu, ref . 1} calculation of the self- 
induced distortion of a two-dimensional circular cylinder immersed in a main- 
stream, with the aim of extending the work into a genuine three-dimensional 
approach, rather than using her time analogue for jet direction. Before 
launching into a three-dimensional perturbation study, however, we felt a need 
to understand more about the validity of any potential-flow type of approach 
and also more about the topology of this particular example of vortex roll-up. 

In 1961 experiments were conducted on a wind-tunnel model of a lifting-fan 
nacelle which generated a 6.4-inch-diameter jet at 200 ft /sec. (See ref. 2.) 

In part of that study, smoke was introduced upstream of the intake to form a 
round jet emerging below the model. The familiar vortex roll-up could be seen 
using the light screen technique and a flow structure was deduced which has 
been confirmed more recently at a very much lower Reynolds number. The mean 
flow seen then strongly resembles the much smaller scale laminar flow shown in 
figure 1. There is therefore some expectation that viscous effects will be 
superposable . 


Figure 2 picks out certain streamlines inside the jet and shows how jet 
fluid becomes part of the trailing vortex. It was pointed out to us by 
Dr. Kroeger of ARO Inc., that the jet (almost) splits into two parts down its 
fore and aft centerline to form one spiral of each trailing vortex. We also 
see in figure 2 that the trailing vortices contain adjacent spirals of jet and 
mainstream fluid. This is further illustrated in figures 3 and 4. Smoke spirals 
of jet fluid in figure 1 and of mainstream air in figure 4 illustrate that, even 
for these essentially laminar flows, mixing is quite efficient. The jet in fig- 
ure 4 is emitted from a representative aircraft model; very similar effects were 
seen using small tubes as for figure 1. 

It is fortunate that the experimentally- found variation of lift interfer- 
| ence force with radius has a maximum value some distance from the cylinder sur- 
I ^ ace * This gives some hope that a meaningful finite vortex representation can 

| be accomplished without the need for an excessive number of vortex elements. 

Nevertheless, we have felt it desirable to study local effects near vortices. 
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In figure 5(a), we see a streamline path for a circle represented by 48 vor- 
tices. Small perturbations, not visible at the scale of the figure, are encoun- 
tered as vortices are passed. If constant length steps are used for tracking, 
the process is quite stable. It is interesting to note that, on starting a 
streamline inside the circle, we find the fluid there is virtually at rest 
(just as Prandtl said it would be). There is slight drift, however, due to the 
fact that the number of vortices is finite, and this allows the streamline to 
penetrate the cylinder surface from inside and join the stable pseudo boundary 
layer" on the outside. 

Figure 5(b) illustrates how the pressure coefficient Cp varies along the 
streamline shown in figure 5(a). The waves in the potential flow curve are due 
to the waviness of the streamline as it passes individual vortices. In fig- 
ure 6 more local errors in the pressure coefficient are examined. It seems 

that a band of about li pitches thick should be avoided. 

We have also looked at local kinematic effects to answer the question "what 
happens to the joining line between vortices?" Figure 7 identifies the area 
we have studied in this regard, which is shown at enlarged scale in figure 8. 
Clearly the thickness of the convoluted layer will decrease as the mesh size 
decreases. One is also tempted to speculate that some sort of mixing process 
is being implied which may be capable of later reinterpretation in viscous 
terms. Even for a two-dimensional model, these convolutions become very com- 
plex; in three dimensions they will be almost unimaginable. 

Much of the work just described was carried out more than a year ago. In 
the interim, we have studied vortex roll-up behind heavily loaded wings, using 
similar techniques (see ref. 3)- Recently we have resumed work on jet plumes; 
a three-dimensional perturbation scheme is being programed in which elements 
with circulation components in the vortex ring direction are used in addition 
to the vertical ones of the two-dimensions-plus-time method. We anticipate 
considerable difficulties in monitoring and understanding the progress of tran- 
sient three-dimensional geometries and have prepared for this by getting 
involved with computer graphics . We are indebted to Dave Smith of our Systems 
Sciences Laboratory, who has been instrumental in the development of our com- 
puter graphics capabilities. Figure 9 shows the familiar H. C. Chang type of 
calculation, using 48 points. In practice, motion is added and the viewpoint 
may be changed almost arbitarily. We believe this will be invaluable when 
interpreting three-dimensional perturbation studies. 
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Figure 1.- A low Reynolds number smoke jet. 
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Figure 3.- Formation of trailing vortices showing entrainment of mainstream 

fluid into trailing vortices. 
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Figure 4.- Smoke visualization of entrainment of main- 
stream fluid into trailing vortices. 
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ONE VORTEX PITCH 



MAINSTREAM 



AREA SHOWN 
IN FIGURE 8 


Figure 7-- Distortion of twin cylinders. 





Figure 9*- Computer-graphic display of the results of an 
H. C. Chang type of calculation. 


48 



EXPERIMENTAL INVESTIGATION OP 


PRESSURES INDUCED ON A FLAT PLATE 
BY A JET ISSUING INTO A SUBSONIC CROSSWIND* 
By Howard M. McMahon and David K. Mosher 
Georgia Institute of Technology 


SUMMARY 


Measurements of static pressure distributions on a flat plate are pre- 
sented for circular and non-circular jets issuing normal to the plate surface 
at various jet exit velocities. Certain trends concerning the pressure 
distribution patterns were noted. It was observed that the pattern for a given 
exit shape and jet velocity is similar to that for a smaller length- to-width 
ratio exit shape at a higher jet velocity. A check was made to see if this 
could be related to the jet plume path. The results indicate that the pressure 
distribution on the plate depends upon the character of the jet as well as 
upon the path of the jet plume relative to the plate. 


INTRODUCTION 


Increasing interest in v/STOL aircraft during recent years has been 
reflected in the reporting of a number of experiments concerned with the 
interaction of a turbulent jet and a subsonic deflecting stream. The inter- 
action of the exhaust from a lift- fan or a lift- jet with the air flowing over 
an aircraft in forward motion can cause large lift losses and moment changes 
during the critical transition phase from vertical to conventional horizontal 
flight. 

Experimental results can be grouped into four categories: 

1) Gross lift and moment measurements on aircraft configurations 
with various jet exit geometries and locations. (See, for 
example, refs. 1 and 2.) 

2) Penetration of the jet into the deflecting stream. (See, for 
example, refs. 3 and 4.) 

3) Measurements of the flow field in and around the jet plume. 

(See, for example, refs. 5 and 6.) 

4) Pressure distribution measurements on the surface from which the 
jet is issuing. (See, for example, refs. 7 and 8.) 


* This work is supported by the U. S. Army Research Office -Durham- under 
contract No. DAHC04 68 C0004. 
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An experimental program is underway at Georgia Tech covering areas 2, 3 
and 4 . Some results have been reported in reference 9 - Further surface- 
pressure distribution results are presented here. 

The interaction between the jet and the deflecting stream is characterized 
by a combination of phenomena - the jet displaces the deflecting stream 
("blockage ") , captures some of the deflecting-stream fluid ("entrainment") , 
and causes a low-energy region to form downstream of the jet plume ("wake") . 
Three different jet exit configurations were used in the present study with 
the aim of achieving large changes in these interaction phenomena, thus leading 
to a better insight into the interference problem. Previously reported data 
on other than circular jets have been confined to penetration or gross lift 
and moment results. The only exception known to the authors is reference 10, 
where pressure distributions on the underside of a wing containing a rectangular 
jet orifice were reported. However, true flat plate conditions were not 
present in this test and no measurements were made away from the surface. Thus, 
a comprehensive study of different exit configurations was felt to be in order. 


SYMBOLS 


a 

m 

V. 

3 

V* 

x 

z 


r j 

Poo 


radius of the circular jet orifice 

ratio of jet exit momentum flux to free-stream momentum flux, 
jet exit velocity 
free-stream velocity 

coordinate in the free-stream direction 
coordinate normal to the plate 
jet exit density 
free-stream density 



EQUIPMENT AND PROCEDURES 


Experiments were conducted in the Georgia Tech 9-foot subsonic wind 
tunnel. Figure 1 is a view of the test section looking upstream and shows 
the plate mounted 12 inches above the bottom of the tunnel on thin support 
struts. The nozzle supply pipe and all of the pressure leads were enclosed 
in a streamlined fairing located under the center of the plate. 

Figure 2 shows the details of the plate construction. The plate is 
3/8-inch-thick aluminum, 48 inches in chord and 66 inches in span, with a 
rounded leading edge. A disc, 32 inches outside diameter and 7 inches inside 
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diameter, contains twelve 0.040-inch-diameter pressure taps along a radius. 

This disc can be rotated remotely to obtain a complete picture of the pressure 
distribution. Inside the disc and bolted to the top of the nozzle supply pipe 
is a replaceable nozzle block. Two nozzle geometries having the same exit 
area were used. The nozzle block containing the 2-inch-diameter circular 
nozzle is shown installed in the plate. The exit geometry of the second nozzle 
(exhibited on the plate surface in figure 2) was a slot 1.0 inches by 2.36 
inches, rounded by a semi-circle of 0.5 inch radius at either end. As shown in 
figure 3, when the longer dimension of this jet exit was aligned perpendicular 
to the free-stream flow direction the exit configuration is termed 'btlunt". 

When the nozzle clock was rotated ninety degrees such that the longer dimension 
was parallel to the free-stream direction the configuration is termed "stream- 
wise . Surface pressures near the jet exit were measured by pressure taps in 
the top of the nozzle block. The components making up the plate surface were 
assembled and aligned within ± 0.010 inches. The velocity profile of the 
boundary layer on the plate was measured on the centerline 15 inches from the 
leading edge with the jet off and was found to follow the 1/7 (turbulent) power 
law, with a boundary layer thickness of about 1 inch. 


Air for the jet was supplied by a 100 
a 6- inch line. The jet exit velocity was 


H.P. centrifugal compressor 
set using the total pressure 


through 

as 


measured in the supply pipe just before the jet nozzle contraction, assuming 
the static pressure at the jet exit to be free-stream static and the nozzle 


discharge coefficient to be 100$. The uniformity of the jet exit flow was 
checked by making total and static pressure surveys in the jet exit plane with 
the tunnel off. Results showed a velocity variation of less than 2$ over the 
exit cross sections. Further details of the test equipment are given in 
reference 9* 


Oil- film flow visualization was done by covering the plate with white 
contact paper and applying a mixture of lamp-black, oleic acid, and diesel oil. 
When a pattern of the oil film became established, the jet and tunnel flows 
were stopped and the result photographed. The main features of the oil pattern 
were found to be repeatable but the intensity of some of the features was found 
to depend on oil mixture and the time allowed for development of the pattern. 

All pressures were measured with a variable-capacitance electric manometer 
calibrated against an alcohol micro-manometer. In general, the plate pressure 
measurements were accurate to within ± 1$ of the free-stream dynamic pressure 
and were repeatable within ±1$. In regions very near the jet and also in the 
near wake behind the jet the pressure readings were fluctuating. The time- 
averaged readings were recorded and were repeatable to within ± 5$ of the free- 
stream dynamic pressure. 


Previous investigators (ref. 7) have shown that for a circular jet the 
primary parameter in the interference problem is the ratio of jet exit momentum 

[p .v . 2 

flux to the undisturbed free-stream momentum flux, / J = m. In order to 


establish the Reynolds number dependence for the blunt jet, plate pressures 
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for this jet were compared at V m = 25 and 50 feet per second while holding 
m = 8. The differences were not significant. All data reported here are for 
a free-stream velocity of 50 feet per second. 


RESULTS AND DISCUSSION 


Two sets of oil-film photographs were taken, one for the jet in the 
crosswind and the other for a long solid body of the same shape as the jet 
exit placed normal to the plate and replacing the jet. The two sets of 
photographs are compared to one another in figures 4 , 5 > and 6 . The free- 
stream velocity in all cases is directed from the top to the bottom of the 
photographs. The oil pattern associated with the jet indicates strong entrain- 
ment of deflecting- stream air into the jet and also into the wake region behind 
the jet. Similar results were observed at m = 4 and m = 12, with the wake 
becoming broader in all cases as m was increased. The oil flow pattern is, of 
course, not necessarily indicative of the flow field off the plate because of 
the presence of the plate boundary layer. 

Figures 7 and 8 show static pressure distributions on the plate surface at 
m = 8 plotted in coefficient form. The static pressure coefficients were 
obtained by measuring the difference in static pressure at a point on the plate 
with and without jet flow and dividing the difference by the free-stream 
dynamic pressure. The static pressure coefficient with the jet off was very 
uniform over the entire plate, the measured values deviating from the average 
value of 0.01 by less than 0.003. Only the pressure distribution over half 
the plate is presented in the figures and is plotted so as to contrast the 
results from the strearawise and blunt configurations with those obtained for 
the more familiar circular jet. The pressure distribution for se”eral cases 
was checked across the span of the plate and the symmetry of the flow with 
respect to the centerline was found to be excellent. 

Figures 7 and 8 show that there is a rearward shift of the low-pressure 
region as the jet exit becomes more blunt. Similar trends were noted at 
m = 4 and m = 12 and the corresponding pressure distributions are shown in 
reference 9 * Additional measurements of plate static pressure have also been 
made at m = 10 for all three jet configurations and at m = 20 for the blunt jet 
configuration . 

When all of the pressure distribution data were assembled in an array of 
contour plots with rows of the array representing increasing m and columns 
representing increasing jet bluntness, a roughly "diagonal" trend was noted in 
the shape of the pressure contour patterns. That is, the shape of the pressure 
contours for the streamwise jet at m = 8 is similar to that for the circular 
jet at m = 12 and also similar to that for the blunt jet at m = 20. A com- 
parison is made in figure 9 to illustrate this trend. It should be noted that 
the observed "diagonal" trend represents a correspondence more in pressure 
contour pattern than in pressure magnitude. 
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Some investigators have suggested that the primary factor determining the 
magnitude and distribution of the surface pressures resulting from the inter- 
ference of a jet with a deflecting stream is the rolling up of the wake into a 
vortex pair downstream of the jet exit. This in turn suggests that the 
position of the vorticity is important in evaluating the interference effect, 
that is, the proximity of the jet plume path to the plate surface may be a 
major factor in the interference as displayed by the surface pressures. One 
possible explanation, then, of the "diagonal” trend noted above is that the 
similar pressure contour patterns for the various jets are a result of similar 
jet plume paths, since it is known (ref. 4) that jets of different initial 
shapes penetrate different distances into the deflecting stream. Measurements 
were made of jet plume centerline paths, here defined as the locus of points of 
maximum total pressure as measured in the jet plume in the vertical p la ne of 
symmetry. The results are shown in figure 10. It is seen that there is no 
correlation between the similar pressure patterns in figure 9 and the corre- 
sponding jet plume paths. To verify this, the jet plume path prediction of 
reference 11 was used to select a value of m for the circular jet which would 
give approximately the same plume path as that for the blunt jet at m = 8. 

The required value was m = 6.8. The two plume paths are shown in figure 10 
and the plate pressure distributions are compared in figure 11. The result 
confirms that the proximity of the jet plume is not in itself a measure of the 
interference pressure distribution on the plate surface. Thus, the combination 
of blockage, entrainment, and wake effects which, in some proportion, specifies 
the character of a given jet plays a major role in the interference problem in 
addition to the role played by the path of the jet plume relative to the surface. 
The character of the jet cannot be well defined nor understood from surface 
measurements alone. Considerable experimental data giving the flow field 
details off the plate are required before a better understanding will be 
achieved. 


CONCLUDING REMARKS 


Oil flow traces and pressure distributions on the plate surface indicate 
large changes in the interference between the jet and the deflecting stream 
when the jet exit configuration is varied. These changes are indicative of a 
complex interplay between blockage, entrainment, and wake phenomena which is 
not as yet well understood. Detailed measurements of the flow field off the 
surface are needed to aid in the understanding. A certain similarity between 
plate pressure contours for the different jets was noted which cannot be 
attributed solely to the proximity of the jet plume to the plate surface but 
which must also depend upon the character of the jet. It would appear that a 
mathematical model must incorporate the features of blockage, entrainment, and 
wake if it is to describe successfully the jet interference problem for various 
jet exit configurations . 
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Figure 1.- Plate installed in test section. Figure 2.- Plate and nozzle block details. 



CIRCULAR STREAMWISE BLUNT 

Figure 5.- Jet exit configurations. 
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Figure k.- Flow patterns on the plate, circular nozzle. 



Igure 5.- Flow patterns on the plate, streamwise nozzle. 







Figure 10.- Measured jet plume centerline paths 





EXPERIMENTAL REACTION JET EFFECTS AT 


SUBSONIC SPEEDS 

By Troy A. Street and Donald J. Spring 

Advanced Systems Laboratory 
of the U. S. Army Missile Command 

SUMMARY 


A brief description of the total overall reaction jet program 
of the Advanced Systems Laboratory is given as well as a brief 
review of the work done to date . A more thorough presentation of 
the subsonic experimental data is made . 


INTRODUCTION 


5 


The purpose of this paper is to give a very brief synopsis of 
the trends in the aerodynamic data that the U.S. Army Missile 
Command has accumulated experimentally in the subsonic regime. 
Predominately, the data collected to date in the subsonic regime 
has been for the sonic jet with underexpanded flow injecting into 
a subsonic freestream. The data presented herein is a collection 
from a test of an axisymmetric body of revolution with a sonic jet 
injecting into a freestream Mach number of 0.285 and a sonic jet 
injecting from a flat plate into freestream Mach numbers of 0.1, 
0.2, 0.4, and 0.6. 

SYMBOLS 

b fin span or slot width, inches 

c fin chord, inches 

Cn normal force coefficient 

Cp pressure coefficient 

D body diameter, inches 

d jet diameter, inches 
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e base of natural logarithm 

F ni normal interference force, lbs 

Ft force due to thrust of the jet, lbs 

K empirical constant 

M freestream Mach number 

m mass rate of flow, lb/sec 

N empirical constant 

P freestream static pressure, psia 

P c total pressure in the jet chamber, psia 

Q jet-to-freestream momentum ratio 

R radius, inches 

V velocity, ft/sec 

Qf angle of attack, degrees 

A increment 

slot length, inches 

6 angle from upstream freestream, degrees 

Subscripts: 
j jet 

00 freestream 

HISTORICAL BACKGROUND 

The U.S. Army Missile Command is currently endeavoring to 
develop those techniques required to predict the aerodynamic 
coefficients necessary for design of reaction jets as control de- 
vices. This stimulates research activities in the dual areas of 
bringing to bear those techniques from other disciplines as well 
as advancing the state-of-art in the solution of fluid mechanics 
problems . 

The current research program of the Advanced Systems Laboratory 
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Of the U.S. Army Missile Command fundamentally draws upon past 
and present analytical and experimental efforts to obtain design 
parameters and prepares for the future to utilize the forthcoming 
generation of computers. Throughout the program the range of 
interest addresses itself to the ultimate objective of developing 
the capability to predict the aerodynamic characteristics of axi- 
symmetric bodies of revolution using multiple reaction jets any- 
where on the body with any jet orientation, jet condition or in- 
jectant, and a wide speed range from low subsonic deep into hyper- 
sonic. Of course, the ultimate objective is so broad that only 
small segments are being examined in logical sequence. 

The current approaches are as follows: 

A. Theoretical Approach 

B . Analytical or Semi-Empirical Approach 

C. Empirical Correlation Approach 

D. Experimental Approach 

The first three approaches will be very briefly described with the 
Experimental Approach in the subsonic regime receiving the majority 
of attention for the interest of this symposium, although con- 
siderably more experimental data has been collected in the super- 
sonic regime by the Advanced Systems Laboratory than has been 
accumulated in the subsonic regime. This is not to de-emphasize 
the Army's interest in the subsonic regime but merely to state 
that the Army, historically, has had more pressing needs in the 
transonic and supersonic regimes. 


APPROACHES 


(a) Theoretical Approach 

In this effort a numerical solution to the Navier Stokes equa- 
tions is being sought for a slot ejecting normally into a super- 
sonic approach flow. Finite difference approximations are used to 
replace the derivatives in the unsteady Navier Stokes equations. 

An initial flow field is assumed and the flow allowed to develop 
with time . The steady state solution is desired and results as the 
limit of the unsteady solution. At present the method is progress- 
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ing normally but steady state has not been reached. An inherent 
drawback in the technique has limited the present effort to ex- 
tremely low Reynolds number and to laminar flow. Future efforts 
will strive to overcome these drawbacks . 

(b) Analytical or Semi-Empirical Approach 

The supersonic analytical approach, as reported previously in 
the literature-*- is based upon the "equivalent body analogy" con- 
cept. This concept is for a sonic, normal, underexpanded jet. The 
jet is assumed to be replaceable by an equivalent solid body that 
will produce a similar effect on the body of revolution. The equiva 
lent body is determined by relating the change in the longitudinal 
component of momentum of the jet to the drag of an equivalent body . 
The equivalent body is assumed to be a sphere-cylinder body with 
the origin of the sphere placed on the jet centerline at the sphere 
radius above the jet exit as shown in Figure 1. The method of 
characteristics is then run for the equivalent body. The inter- 
section of the characteristics net with the missile body provides 
a surface for the integration of pressure to determine the inter- 
ference forces and moments. A comparison of the data obtained by 
this technique with the experimental data for the incremental 
normal force coefficient variation with Mach number for various 
pressure ratios is shown in Figure 2. The trends obtained by this 
equivalent body analogy are correct while the magnitude appears to 
be in error from 11 to 16%. 

The subsonic analytical approach is described in another paper at this 
symposium by Robert Rosen, Norbert A. Durando, and Louis A. Cassel of the 
McDonnell Douglas Corporation, which is under contract to the Advanced 
Systems Laboratory to develop this subsonic approach. 

(c) Empirical Correlation Approach 

The empirical correlation approach used the force data obtain- 
ed during the experimental tests for Mach numbers 0.8 to 4.5 for 
various pressure ratios and axial jet locations on a body of re- 
volution. 

The augmentation ratio, Kjj, was selected as a pertinent para- 
meter to measure the control force . The augmentation ratio is 
defined as 


- 1 + F ni/ F T 
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where F ni is the normal interference force and F T is the jet thrust. 

The procedure for developing the empirical technique was to 
plot the augmentation ratio, K n , as a function of jet momentum-to- 
freestream momentum, Q . The curves were linearized by multiply- 
ing the augmentation ratio by the momentum ratio. The results were 
then curve fit . The curves were then rotated to a common slope and 
shifted to a common origin. 

Two empirical constants, K and N, are required at this point. 

K is the empirical constant required to rotate the curves and N is 
the empirical constant required to translate the curves to a common 
origin. These two empirical constants are functions only of the 
type of jet, slot or circular, and the body axial location. Figure 
3 shows the calculated values of the empirical technique compared 
with the corresponding experimental data to be within an accuracy 
of + 10%. 

(d) Experimental Approaches 

The models used for the experimental reaction jet programs 
consisted of a 4-caliber tangent ogive nose with a cylindrical 
afterbody. Figures 4 and 15, for tests at Mach numbers of 4.5 to 
0.8 and 0.285 and a flat plate at Mach numbers of 0 .6, 0 .4, 0 .2, and 
0.1. Force and moment tests were conducted for the Mach number 
range 4.5 to 0.8 and a flow field survey conducted for the Mach 
number range 1.2 to 0.9? the flat plate was used in a pressure test. 

The force and moment tests were for angle of attack of + 4° 
using the model shown in Figure 4. Either air or nitrogen was used 
as the secondary gas and was injected normal to the surface at a 
sonic exit Mach number. Artificial boundary layer trips were in- 
cluded in all tests to assure a turbulent boundary layer. 

Figure 5 shows a typical variation of the normal force coef- 
ficient with angle of attack for various jet total-to-freestream 
static pressure ratios as a variable. Two features should be noted 
concerning these figures: 

(1) Over the angle of attack range + 4 degrees, these 
data vary linearly with no discontinuities. 

(2) These data at the zero angle -of-attack intercepts 
represent the augmentation or degradation of the 
control force . 

These data are in the transonic region where information is sorely 
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needed . 


Another significant program investigated the jet on stabiliz- 
ing surfaces located aft of the control jets^, with the testing 
performed at Mach numbers of 0.9 and 1.2 for several jet pressure 
ratios. An example of the changes due to jet momentum ratio is shown 
in Fig. 6. A small (%-inch dia.) flow probe and pitot-static tube 
was used to map the flow fields in the immediate vicinity of the 
fin, and the resulting data were reduced in terms of velocity com- 
ponents in the X,Y,Z plane . Figure 7 shows a representative sur- 
vey. The data from this series of tests are being used to determine 
the flow field properties downstream of a control jet, and are be- 
ing applied in the efforts to develop a prediction technique. 

Recently the NASA Langley Research Center and the U.S. Army 
Missile Command's Advanced Systems Laboratory began a joint effort 
to study the effects of a predominantly sonic jet injecting from 
a flat plate into a subsonic freestream. The joint study consists 
of an initial test phase that was completed in August 1969, a final 
test phase to be completed later in 1969, and an analysis of the 
data collected. The initial test phase was a flat plate pressure 
test, while the final test phase will consist of a rake survey of 
the jet wake. 


The subsonic flat plate test data recently collected at Langley 
Research Center has not been fully analyzed. However, a presenta- 
tion of some of the significant aspects so far unveiled is present- 
ed in Figures 8 through 14 . The flat plate was oriented at zero 
angle of attack for the tests and each test point was repeated 
twice. Figure 8 shows the effect of pressure coefficient, Cp, vari- 
ation with radial length, R, for a 0.33-inch diameter sonic nozzle 
on the upstream ray (9 = 0) as well as the variation of C p at each 
orifice with the number of times the data was repeated. A slight 
positive Cp immediately in front of the jet exists at this very low 
freestream Mach number, 0.1, and momentum ratio ( jet-to-freestream 
ratio, Q) of 2.0. When comparing Figures 8 and 9 for the same mo- 
mentum ratio but different Mach numbers (0.1 and 0.6, respectively), 
it becomes apparent that the scatter in the upstream ray of Figure 
8 is that due to Mach number. At the lower Mach number the free- 
stream flow around the jet fluctuates considerably. Most likely 
the fluctuations exist at the higher Mach numbers but at a suf- 
ficiently high frequency that their effect is damped in the length 
of tubing from the static orifice to the pressure transducer located 
outside the tunnel. 
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Additional pressure distribution on rays 9 = 90° and 150° for 
the Mach number 0.6 and the momentum ratio 2.0 conditions are shown 
in Figures 10 and 11 . The 9 = 90° ray shows only a small increase 
m C p near the nozzle. The 9 = 150° ray shows a negative C D near 
the nozzle that is caused by either entrainment or blockage or both. 

Figures 12 through 14 show the effect of higher jet pressure 
which is reflected in the momentum ratio increase to 100. Figure 
12 shows a typical upstream pressure distribution with a signifi- 
cant increase m Cp near the nozzle. This is because of the block- 
age effect at the higher Mach number for this condition. 

Figures 13 and 14 show the pressure distribution for the 90° 
and 150 rays, respectively. For this Mach number, 0.6, the lowest 
value of pressure coefficient exists at the 90° ray at the pressure 
orifice nearest the nozzle . Next examine the 150° ray and note the 
lowest pressure exists an inch from the nozzle and some pressure 
recovery near the jet has been manifested between the 90° ray and 

the 150 ; also the magnitude of pressure coefficient is less on the 
150 ray. 

Comparing Figures 9 through 11 with Figures 12 through 14, respec- 
tively, shows that increasing the momentum ratio effectively moves the low- 
est pressure (originally behind the jet) forward. The lowest pressure remains 
near the jet. As this lowest pressure moves forward, a recovery process 
apparently originates near the jet and moves outward along the rays. 

The above comparison suggests that a tremendous viscous inter- 
action process is displayed in these adjacent areas as the high 
pressure region in front of the jet attempts to rush into low 
pressure regions downstream of the jet. This in turn would suggest 
that considerable mixing occurs not only in the freestream itself 
but between the freestream and the nozzle injectant. 

The subsonic experimental data 3 was accumulated on the axi- 
symmetric body, shown in Figure 15, with the pressure orifices. 

The model is 2.25 inches in diameter, with rows of orifices 30° 
around the body from the jet. Both circular nozzles and radial 
slot nozzles were used. The effects of a sonic, circular jet on 
the longitudinal pressure distribution is shown in Figures 16 
through 19. Figure 16 shows a typical pressure distribution on 
the axisymmetric body. (The rows of pressure coefficient distri- 
bution with longitudinal length X measured from the jet correspond 


to the rows shown in Figure 15. The row through the jet centerline 
is the middle row with the discontinuous distribution at the jet 
location.) There exists a small pressure rise upstream of the jet 
on the jet centerline. However, the dominating factor is the low 
pressure in back of the jet and on the adjacent rows of orifices. 

As shown in Figures 16 and 17 the pressure profile changes very 
little with angle of attack. An increase in momentum ratio fail- 
ed to change the pressure distribution; it only succeeded in re- 
ducing the upstream positive pressure and reducing the value of 
negative pressure as shown in Figure 18. 

The effect of a sonic, slot nozzle was to increase the positive 
C D ahead of the jet and reduce the negative Cp behind and on the 
sides of the jets. Figure 19. The reduction on the adjacent rows 
is because the slots are much closer to the rows than when the jet 
was circular. The same throat area existed for the slot as for 
the circular nozzle. The injectant is spread over a large region 
thereby causing a slightly greater disturbance to the flow field. 
This, in turn, manifests a slightly greater change in aerodynamic 
coefficients . 

It was not felt that these various changes in chamber pressure, 
angle of attack, or slot configurations changed the aerodynamics 
of the body by a significant amount. The conclusion was that, for 
the range of parameters tested, the sonic jet dominated the Mach 
number 0.285 freestream so overwhelmingly that no other test vari- 
able emerged except in a very weak manner. Logically, it would 
appear that, because of curvature effects, the jet would dominate 
the flow over an axisymmetric body more than it would for a flat 
plate . 


CONCLUSIONS 


The following conclusions are drawn: 

1 . The Army Missile Command has been successful in the super- 
sonic regime with a pure blockage model replacing the jet . 

2. In the subsonic regime, based on the axisymmetric body 
and the flat plate tests, the dominating factors appear to be (a) 
the jet in the low subsonic regime and (b) viscous interaction in 
the higher subsonic regimes . 
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Figure 5.- Variation of experimental normal force coefficient with angle of attack. 





INCHES 


CHORD 2 I.5D 

_b_ 

I.OD 
0.75D 
0.50D 



(mVJj/lmVloo 
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momentum ratio for several rectangular downstream fins. 



Figure 7.- Velocity components at downstream fin location obtained 
in flow-field survey. 
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Figure 9.- Flat-plate pressure distribution at M = 0.6 for 0° ray. Q = 2, 
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Figure 14.- Flat-plate pressure distribution at M = 0.6 for 150° ray with increased momentum ratio. 
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A BLOCKAGE- SINK REPRESENTATION OF JET IN TER FERENCE 
EFFECTS FOR NONCIRCULAR JET ORIFICES 
By J. C. Wu and M. A. Wright 
Georgia Institute of Technology 


SUMMARY 


A study of the subsonic jet-in-crosswind problem is made by using a two- 
dimensional potential- flow representation of jet interference effects. Analyses 
are made for two different jet exit configurations: an ellipse with its major 

axis alined with the free stream flow direction, and an ellipse with its minor 
axis in the free-stream direction. In each of the configurations the jet issues 
normally from the surface of an infinite plate into the crosswind. The model 
uses a blockage-sink representation of the jet interference effects on the cross- 
wind. xne blockage element consists of a jet-exit-shaped cylinder with an 
attached afterbody and is used to simulate the effect of the jet-wake combina- 
tion in displacing the crossflow. The sink element is placed in the afterbody 
to simulate the combined effect of entrainment by the jet as well as by the wake. 
Calculated pressure distributions on the plate are compared with ejqoerimental 
data obtained with similar jet exit configurations. The results show that the 
two-dimensional representation is useful for the case of high jet-to-crosswind 
speed ratio. 


INTRODUCTION 


This paper treats the problem of a single subsonic turbulent jet issuing 
normally from the surface of an infinite plate into a subsonic crosswind. It 
has been observed that the presence of the jet in the crosswind leads to a 
separation of the crosswind behind the jet. The combined effects of aerodynamic 
forces acting on the jet and the entrainment of the crossflow fluid into the jet, 
which carries momentum with it, causes the jet to spread, deform, and deflect 
after leaving the plate. In turn, the jet interferes with the flow of the cross- 
wind by displacing it, entraining it, and causing the separation region to form; 
these three factors are usually referred to as the blockage effect, the entrain- 
ment effect, and the wake effect, respectively. Because of the complexity of the 
jet interference phenomena, there is presently no complete theory treating the 
detailed flow in this problem. In lift- jet or lift-fan powered v/STOL applica- 
tions, the change in surface pressure due to the jet interference is of primary 
concern. During the past few years, several analytical flow models (refs. 1, 2, 
3) have been proposed as means of representing the jet interference effects, 
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and from these representations, the interference surface pressure can be calcu 
lated. The models have the common feature of requiring the use of either empir- 
ical or semiempirical information concerning the path of the jet in the cross- 
wind. Potential flow elements, that is, vortices or a combination of doublets 
and sinks, are placed on or near the center line of the jet plume. The velocity 
induced by the potential elements on the surface is used to calculate the inter- 
ference pressure. The models differ in the factors that are accounted for (ot 
neglected) and in the manner in which these factors are represented. In two of 
the proposed models the effect of entrainment is ignored. To the authors 
knowledge, no numerical results have been reported that use these m °3-<~ls for 
other than a circular jet exit configuration. Other models may result from the 
extension of the work of Lu Chang (ref. 4) concerning the time- dependent def or 
mation of a circular column of fluid in a crossflow. This extension is Leing 
considered at the Georgia Institute of Technology as well as other places, but 
results are at present inconclusive. 

The principal objective of this paper is to demonstrate that for the case 
of high jet-to-crosswind speed ratio, it is possible to use a two-dimensional 
approximation with a blockage-sink representation to study the jet interference 
effects with noncircular as well as circular jet exit configurations. The case 
of high-speed ratio is of particular importance for lift-jet or lift-fan powered 
V/STOL aircraft in the transition phase of a flight. Results presented are for 
elliptical jet exit configurations, the major axis of the ellipse being either 
alined with or perpendicular to the free-stream direction. 


SYMBOLS 


constants in conformal transformation, eq. (8) 
integrated pressure coefficient 
pressure coefficient, eq. (l) 
diameter of a circular jet exit 
half -width of Rankine oval 
sink strength for entrainment 
half-length of Rankine oval 
local static pressure 
free stream static pressure 
r plate radial coordinate 

R, R radial limits of integration, eq. (l8) 

S strength of source-sink pair in Rankine oval 
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u 

V 

Uoo 

u j 

X 

s 

z = x+iy, 

P 

C = 5+LT1, 

P 

U) 

Subscripts : 


local velocity in direction of free stream 
local velocity normal to free stream 
free-stream velocity 
jet velocity 

sink location for entrainment 
plate complex coordinates 
plate angular coordinate 
Rankine oval complex coordinates 
density 

complex potential 


b blockage-induced velocity components 

s sink-induced velocity components 

Superscript : 

( ) ' C -plane velocity components 


ANALYSIS 


The blockage-sink representation has the advantage of being relatively sim- 
ple and easily interpretable. In reference 1 an analytical model based on this 
representation was developed. A doublet distribution along the jet center line 
was used to simulate the blockage effect due to the jet plume, and a sink dis- 
tribution was used to simulate the entrainment of crosswind fluid into the jet. 
The authors of reference 1 stated that because of the manner in which the block- 
age effect was simulated, the representation of the flow near the jet was not 
expected to be good. They showed that with a suitable choice of entrainment 
coefficients, the model gives a calculated surface pressure field in good agree- 
ment with experimental data for a circular jet issuing normally from the center 
of a large rectangular wing placed at zero lift and with a jet-to-crosswind 
speed ratio of 10. Neither theoretical nor experimental results were presented 
in reference 1 for a midspan of the wing covering a distance of 10 jet radii. 
Experimental data from other sources (refs. 5 und 6), however, have shown that 
the interference pressure on this portion of the surface contributes substanti- 
al ^he overall interference effects of lift loss and nose-up pitching 
moment. 
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In the present approximation the effect of variation of flow between planes 
parallel to the plate is considered unimportant and ignored as far as the calcu- 
lation of surface pressure is concerned. It is noted that the jet interference 
problem is basically three-dimensional. Nevertheless, the simplifications and 
versatility resultant from its use makes the two-dimensional approach attractive. 
For example, the representation of noncircular jet exit shapes with the use of a 
three-dimensional model involves considerable complications. To the authors 1 
knowledge, no analytical results using a three-dimensional model for other than 
the circular jet exit configuration have been reported. As shall be shown, with 
a two-dimensional model the generalization to noncircular jet exit shapes is 
very straightforward. The question remains as to what is the region of validity 
of the approximation. In reference 5 the authors presented some results using a 
two-dimensional model for a circular jet. Many of the observations presented 
in reference 5 are also pertinent to noncircular jets. These results are not 
repeated here, but the major conclusions are summarized below. 

Experimental observations indicate that if the jet-to-crosswind speed ratio 
is sufficiently high, the jet deforms and deflects significantly only after it 
has penetrated several jet radii into the crosswind. Since the velocity induced 
by a three-dimensional blockage or sink element at a given point in the flow 
field is inversely proportional to the square of the distance between the ele- 
ment and the point, three-dimensional effects on the interference surface pres- 
sure are expected to be small for the high- speed- ratio case, particularly in 
regions close to the jet exit. A surface pressure distribution was calculated 
by using a large number of three-dimensional doublets and sinks placed along an 
empirical jet-center-line path for a speed ratio of 8. The result was shown to 
agree very well with that obtained with an equivalent two-dimensional represen- 
tation. Thus, for the circular- jet case, the effect of jet deflection was con- 
sidered to be unimportant for speed ratios in excess of 8. The same conclusion 
is expected to be valid for noncircular jets, as long as the jet penetration is 
not much less than that for a circular jet at a speed ratio of 8. 

Two other questions were also discussed in reference 5* The first concerns 
the justification of using a potential model to represent the jet interference 
effect on the crosswind, and the second concerns the relative importance of the 
entrainment effect. It was felt that, outside of the jet plume and its wake, 
the crossflow is essentially inviscid and the use of a potential flow representa- 
tion is justified. However, inasmuch as there exists an extensive wake region 
of low total pressure, the representation of the entire crossflow by a constant- 
total-pressure potential-flow model is too strong a requirement on the model. 
Accordingly, a region downstream of the jet covering a total angle of 90° was 
excluded from the surface pressure calculations. With the exclusion of this 
wake region, it was shown that the entrainment effect, contrary to some previous 
conclusions, contributes significantly to the jet interference. 

It should be noted that the two-dimensional approximation used here is not 
related to the two-dimensional slotted jet problem which is also of current 
interest. 

Aside from the simplifications resulting from the use of the two-dimensional 
approach and the subsequent extension to noncircular jet exit shapes, the present 
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model differs from the blockage -sink model used in reference 1 in that, in addi- 
tion to representing the effect of jet plume blockage, an attempt is made to 
simulate the effect of the wake blockage by attaching an afterbody to the cylin- 
der. It has been shown in reference 6 that the pressure distribution aro un d a 
solid circular cylinder calculated from potential theory accounting for only the 
cylinder blockage differs greatly from experimentally observed pressure distri- 
bution on a plate to which a solid circular cylinder is attached perpendicularly. 
Figure 1 presents a comparison of the calculated and experimentally obtained 
constant-pressure contours shown in terms of the pressure coefficient C 
defined by p 



( 1 ) 


The dashed curves are calculated contours and the solid curves are experimental. 
Since the effect of plate boundary layer on the plate pressure is of secondary 
importance (ref. 6), it is suspected that the effect of blockage of the crosswind 
by the wake is of considerable importance. Accordingly, constant-pressure con- 
tours on the plate were calculated by use of a solid cylinder- afterbody repre- 
sentation and compared with the experimentally measured pressure distribution 
around a circular cylinder in figure 2. Again, the dashed curves are calculated 
contours and the solid curves are experimental. The shape of the afterbody used 
in the calculation is also shown in figure 2. The manner in which this particu- 
lar afterbody shape is obtained will be discussed later. The agreement between 
the two sets of contours is satisfactory outside the excluded wake-flow region. 

It is noted that the solid circular cylinder case, aside from the plate boundary- 
layer effect, is truly two-dimensional. A comparison of figure 1 with figure 2 
shows that the addition of the afterbody significantly improved the blockage 
representation. 


In addition, the present analysis considers the entrainment of crosswind 
fluid to take place mainly through the wake rather than along the jet periphery, 
as was assumed in reference 1. There is some evidence supporting this approach 
(ref. 7), although detailed information concerning the detailed entrainment 
mechanism involving nonparallel streams is still lacking. It is recognized that 
the entrainment effect is better represented by a distribution of sinks along 
the jet and wake periphery, but in the absence of empirical information, a single 
sink placed in the afterbody is used to simulate the lumped entrainment effect. 
Preliminary studies using a distribution of sinks along the line of symmetry 
were made. It was found that the resulting interference pressure did not differ 
greatly from that obtained by using an equivalent single sink. As in refer- 
ence 5, no attempt is made here to simulate the flow in the wake or to calculate 
the interference pressure in a region where viscous effects are thought to be 
highly important. Because of the present uncertainty concerning the character 
of the wake, the development of an analytical model representing the wake flow 
must await further detailed experimentation. 

With the potential- flow model, the blockage- induced velocity and the 
entrainment-induced velocity are superpositioned to give the total interference 
velocity. The simulation of the blockage effect of the jet-wake combination is 
obtained from the flow around a transformed Rankine oval. The oval, formed by 
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placing a source of strength S at the point (|p,0) in the |-plane (|,q) and a 
sink also of strength S at fa’ 0 )’ wlth h > ^1> in a uniform stream U^, in 
the positive ^-direction, is shown in figure 3- The complex potential for this 
flow is 


cd = -U^ - S In (£ - |p) + S In (£ - l 2 ) 


and the complex velocity is 


" u b + iv b 


dtD 

d£ 


q, - g + s 

5 - h i - i 2 


( 2 ) 


(5) 


where Up and Vp are, respectively, the |- and ^-components of the flow veloc 
ity, the prime denoting that these quantities are the components in the |-plane 
and the subscript b denoting that the velocity results from blockage. 


It is straightforward to show that the oval is symmetrical about both the 
|-axis and the line | = (|p + iz)/ 2 -- The half-length of the oval is 

1/2 


„ 1 
n = 2 


\h - 5l f + *(*2 ‘ 


(4) 


and thus the upstream stagnation point is ( 1^,0 ), |j being given by 


t 


5 " 


h. + h 

2 - n 


(5) 


The shape of the oval is given by the equation 

iia : Jib = tan M (6) 

|2 + T)2 _ Ip - | 2 + ^1^2 s 


and k, the half-width of the oval, is given by 

" - Si ) 1 ( 

kn 2 - fa - |p) _ 

If any two of the four parameters, | 2 - iq, S/lJ^, n, and k, are given, 
the shape of the oval is determined. 

Consider now a solid line segment extending from the point (I4, o) to the 
upstream stagnation point of the oval ^|^,0) and coinciding with the |-axis, 
with I4 < 0 and ||^| > ||^| as shown in figure 3. The line is taken to be 
attached to the oval. Since the flow past the oval is symmetrical about the 
|-axis, the presence of the solid line does not disturb the flow. The 
transformation 


b'k 2 - ^| 2 - iq)^ - ^l 2 - |p) ctn 
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z = A£ + B ft? 


(8) 



maps the solid line segment in the £ -plane into a portion of an ellipse in the 
z-plane as shown in figure The semiaxes of the ellipse, which has its center 
coinciding with the origin of the z- coordinates, are I4A and I4B. The semi- 
axis 54A coincides with the x-axis and may be either major or minor. The flow 
past the line segment plus the Rankine oval in the £ -plane corresponds to the 
flow past the ellipse with an afterbody in the z-plane, the shape of the after- 
body being that of the Rankine oval transformed in accordance with equation (8). 
The size of the ellipse is determined by the value of I4. By varying the 

parameters k and n, a variety of afterbody shapes are obtained. The points 
of attachment where the afterbody is joined to the ellipse are given by 


x — A|^ 


(9) 


y = ±B /sir - t 5 2 

The complex velocity in the z-plane is given by 


( 10 ) 


" u b + iv b 


dm _ dm / dz 

dz df;/ d£ 


( 11 ) 


where ^ 


is given by equation (3) and 


M = a + _JL_ 
Ji 2 - 


( 32 ) 


The computation of u^ and v-^ in the z-plane from equation (ll) is easily 
performed by first finding the values in terms of a given value of £, and then 
determining the value of z corresponding to that value of £ obtained by use 
of equation (8). 

The entrainment effect is simulated by placing a point sink of strength m 
in the z-plane at the point ^x g ,0^ with x s > I4A. The complex velocity due to 

the sink in the z-plane is simply 


-Ug + iVg 



(13) 


where the subscript s denotes sink- or entrainment- induced velocity components. 
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The total velocity components due to the cylinder-afterbody blockage and 
the sink are given by 


u = u b + u s (l*0 

v = v b + v s ( 15 ) 

For a constant- total-pressure potential flow, the pressure coefficient Cp, 
with It, assigned the value of unity for simplicity in computations, is 

C p = 1 - u 2 - v 2 (l6) 


It is noted that the case of a circular- shaped jet exit is a special case 
of the above. The transformation (8) for this special case becomes, with 
A = l/2, B = l/2, and = d, the diameter of the circle, 

z = + ^ 2 - d2 ) (it) 


Results for a circular-shaped jet have been presented in reference In 
the present paper, interference surface pressure fields are calculated by using 
the two-dimensional blockage- sink model for an elliptically shaped jet exit with 
a major to minor axis ratio of 3-3 6. This ratio is selected so that the calcu- 
lated results may be compared with experimental results using jet exit configu- 
rations approximating the elliptical shapes. The experimental data were obtained 
in the Georgia Institute of Technology 9-foot low- speed wind tunnel. A descrip- 
tion of the experiments is given as a part of this Proceedings, in a paper 
authored by H. M. McMahon and D. K. Mosher. It is noted that the experimental 
jet exits are not elliptical, but are slotted with a total length of 3-36 inches 
and a width of 1.00 inches, as shown in figure 5* In "the calculations, is 

selected so that the major and minor axes of the ellipse have, respectively, the 
values 3*36 inches and 1.00 inches. The area of the ellipse, however, is 
approximately l6 percent smaller than that of the jet exit in the experiments. 

Calculated pressure contours for the elliptically shaped jet with the major 
axis of the ellipse alined with the free-stream direction are compared with 
experimental results for a speed ratio of 12 using the slotted jet exit with the 
length of the slot oriented streamwise. Figure 6 shows this comparison, with 
the calculated contours shown as dashed curves and the experimental contours 
shown as solid curves. It is seen that, with the choice of the parameters as 
given on figure 6, the two sets of contours agree well with each other over a 
large region on the plate. The calculated surface force distribution around the 
jet, represented by the integrated pressure coefficient 
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( 18 ) 



where R Q is the contour of the jet exit, is compared with experimental data in 
figure 7. Again, good agreement is observed. 

Similar comparisons are given in figures 8 and 9 for a jet with its longer 
dimension perpendicular to the free-stream direction. The calculated results 
are for the elliptically shaped jet and the experimental data are for the 
slotted jet at the speed ratio of 12. Again, good agreement is observed. 

Results for speed ratios other than 12 were also obtained, although not pre- 
sented here. In general, the calculated results, with a suitable choice of 
parameters, agreed reasonably with the experimental data over a large region 
on the plate. 

In selecting numerical values for the parameters used in the calculations, 
available experimental contours are used as a guide. Thus the present model, 
like previously proposed models, requires empiricism. The calculation procedure 
required in the two-dimensional model, however, is sufficiently simple to permit 
a parametric study from which the dependence of these parameters on the jet-to- 
crosswind speed ratio and on the jet exit shape may be obtained. This work is 
presently in progress. 


CONCLUDING REMARKS 


The present work is motivated by a need to establish a better basic under- 
standing of the jet interference phenomena. Attempts are therefore made to 
select a simple and yet physically interpretable model to represent the jet 
interference effects. In this regard, the two-dimensional blockage-sink repre- 
sentation appears to be successful. In view of the perhaps drastic simplifi- 
cations introduced, the observed agreement between the calculated interference 
pressure field and the experimental data for noncircular as well as for 
circular jet exits is encouraging. Some dependence on empiricism is presently 
unavoidable in developing any analytical model with the hope of representing 
the very complex interference flow. In the present analysis, this dependence 
leads to the selection of the parameters n and k, which are, respectively, 
related to the length and width of the wake, and to the selection of m and x , 
which are related to the entrainment effects. In the calculation of the s 
interference pressure outside of the wake-flow region, the effect of the wake 
blockage length does not appear to be important. This has been verified by 
computational results. The most pertinent parameters are therefore k, x g , and m. 
The dependence of these three parameters on the jet-to-crosswind speed ratio and 
on the jet exit shape needs to be obtained experimentally. Preliminary results 
based on the parametric study using the two-dimensional model, however, lead 
to certain conclusions which appear to be physically reasonable. For example, 
the trend concerning the entrainment rate is that, with the crosswind speed 
held constant, the amount of crosswind fluid entrained increases with increasing 
jet speed. The conclusion that the entrainment takes place mainly through the 
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wake (ref. 7) is substantiated by the analysis in that the results with the 
point sink placed inside the afterbody led to better agreement with experimental 
data than the results with the sink placed at the center of the jet. It was 
also found, contrary to the indication of oil film photographs (ref. 5), that 
the use of a narrower wake, that is, afterbody, at higher speed ratios gives 
calculated results in better agreement with experimental data. This may be 
explained by the fact that the jet entrains a greater amount of low-energy 
crosswind fluid at higher speed ratios and hence delays the separation. Addi- 
tional parametric studies using the two-dimensional model are in progress. It 
is hoped that these studies will lead to additional information enabling better 
predictions of interference effects and perhaps areas of experimental emphasis 
as well as possibilities of alleviating the adverse interference effects. 

For a jet with a speed ratio significantly less than 8, three-dimensional- 
effects are expected to be important because of the rapid deflection of the jet. 
The two-dimensional model described here provides a basis for the development 
of a more refined three-dimensional model. The extension will require the 
specification of a number of additional parameters, such as the variation of 
the entrainment rate along the jet path, the deformation of the jet, and the 
wake geometry along the jet path. Furthermore, experiments indicate that far 
downstream of the jet exit, the jet plume deforms to a kidney shape containing 
a pair of vortices . It is uncertain at the present time whether the effect of 
this vortex pair is of primary importance to the interference pressure. This 
question as well as the formation of and detailed flow in the wake region appear 
to be areas requiring further investigation. 
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Figure 1.- Constant C p contours around a circular cylinder. 
Potential flow - circular cylinder only. 
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Figure 2.- Constant C p contours around a circular cylinder. 
Potential flow - circular cylinder with afterbody. 





Figure 4.- The z -plane. 


ellipse 
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EXPERIMENTAL JET EXIT 



Figure 5.- Shapes of jet exit. 
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o EXPERIMENTAL DATA, Uj / 12 (REF. 5) 

POTENTIAL FLOW, m = 1.22, x $ = 2.0 

r?= 5.2, k = 0.218, £ 3 = - 1.78 

Figure 7.- Surface force distribution around streamwise jet. 
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m = 2.26, x s = 5.0 
n = 5.2, k = 1.139, £3 


0.674 


Figure 8.- Constant Cp contours around blunt jet exit. 



o EXPERIMENTAL DATA, Uj / = 12 (REF. 5) 

POTENTIAL FLOW, m = 2.26, x $ = 5.0 

n = 5.2, k = 1.139, £3 = 0.674 


Figure 9.- Surface force distribution around blunt jet exit. 
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DEVELOPMENT OF AN ANALYTICAL MODEL FOR THE 

FLOW OF A JET INTO A SUBSONIC CROSSWIND 

By P. T. Wooler 

Northrop Corporation 
Aircraft Division 


SUMMARY 


The interaction between a jet exhausting at a large angle from a surface into a 
subsonic crosswind is discussed. The problem of determining the position of the jet 
centerline is considered first. An analytical model is discussed which represents 
the deflecting mechanism by mainstream entrainment and a pressure force. Two 
methods for obtaining the jet induced velocity field are explored. The first, a 
vorticity model, is shown to be useful for V/STOL preliminary design. The second 
method, a sink-doublet model, describes the flow in more detail and is shown to pre- 
dict wing surface pressure distributions which agree very well with test data. 


INTRODUCTION 


There is considerable interest in airplanes with the capability of taking off and 
landing in a short distance. One method of achieving such performance is to use lift- 
ing jets or fans, installed in the wing or fuselage of the airplane and exhausting at right 
angles to the direction of flight. 

In the transition phase of such flight, important jet on airstream interference 
problems arise. It has been demonstrated experimentally that in the transition phase 
the jet induces an asymmetric pressure loading on the surface from which it is 
exhausting. The pressure loading can adversely affect the lift and pitching moment on 
the airplane, usually in the form of significant losses in the lift and an appreciable 
increase in the nose-up pitching moment. The development of a satisfactory theoreti- 
cal model of this type of flow field is therefore desirable to assist in the design of 
V/STOL airplanes and to enable a satisfactory wind-tunnel correction method to be 
developed for powered-lift wind-tunnel testing. 
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In this paper the development of an analytical model for the flow of a jet into a 
subsonic crosswind is discussed. The development of a method for calculating jet 
centerlines is considered first and then the problem of determining the jet crossflow 
interference velocity field is investigated. 
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SYMBOLS 


A. Jet cross-sectional area 

b Wing span 

C Circumference of jet cross section 

Crossflow drag coefficient of jet 

c Wing chord (root chord for delta wing) 

C Pressure coefficient 

P 

d Length of major axis in elliptical representation of jet cross section 

E Mass entrainment of mainstream fluid into jet per unit length of jet 

E^, E 9 , Eg Entrainment coefficients 

2 

K Jet area factor, = Kd 

L., M. Interference lift and moment (value with jet on - value with jet off) 

m Velocity ratio, ^jo/U 

nf Sink strength of jet element 

R Local jet radius of curvature 

s Coordinate along jet centerline 

S Wing surface area 

T Thrust of jet 

U Mainstream speed 

U. Jet speed 

1 

W Complex velocity potential for ellipse 

x, y, z Coordinate system defined in Figure 1 

Xj x-coordinate of center of exhausting jet 

y Coordinate in complex plane 

d Angle defined in Figure 1 
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P Doublet strength 

p Jet and mainstream density 

r Jet circulation per unit length 

Subscripts : 

0 Refers to conditions at the jet exit 

1 Refers to interference value (value with jet on - value with jet off) 

THE JET CENTERLINE 


When a jet exhausts at an angle into a crossflow, it is deflected into the 
crossflow direction due to the effects of viscous entrainment of the crossflow by the 
jet and also due to crossflow induced pressure forces on the jet boundary (Figure 1). 

A method for calculating jet centerlines and the jet induced velocity field has been 
presented in Reference 1. It is assumed in this model that the mass entrainment E of 
the crossflow into unit length of the jet is given by 

p E„ (U. - U sin 9) C 

E ’ » E 1 Ud C0S + 1 4 E„ U COS W. (1) 

J J 

It is also assumed that the pressure force on a jet element of unit length is 

C D f PU 2 cos 2 0d (2) 

The equations of continuity and tangential momentum for the jet may then be 
deduced by taking IL to be the mean jet velocity. The centrifugal force is equal to 

the mass of air entrained times the normal free stream velocity component plus a 
drag term given by equation (2). This force is equal to the rate of change of the jet 
momentum in a direction normal to the jet centerline. 

The equations of motion for the jet development are then 


P C (AjU.) - E 
Pdf’ (AjU 2 ) = EU sin 9 


( 3 ) 

( 4 ) 
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( 5 ) 





,2 

d x 




p Ud cos 0) U cos 9 


To solve these equations.it is necessary to postulate the cross-sectional shape 
of the jet. It is assumed that the jet changes linearly from the circular cross section 
to an elliptical section of 4 to 1 major to minor axis ratio in a predetermined distance 
proportional to the initial jet diameter and to the initial jet velocity and retains this 
section thereafter. This assumption permits Aj to be written as a function of d, the 
jet major axis . 

Equations (3), (4), and (5) are then a set of differential equations for the 
variables x, U., d as functions of the variable z and the parameters E^, E^, E^, and 


is taken to be 1.8, a representative value for the drag coefficient of an ellipse. 
E_ is set equal to . 08 to give correlation with the data of Reference 2 for the free jet 
case. 


The entrainment coefficients E-, E„ are chosen to give correlation with available 

13 ^ 

test data for jet centerlines. This correlation is demonstrated in Figure 2. Although 
the comparison in Figure 2 is for jets exhausting normally to a wall, equations (3) 
through (5) may be solved for other values of the angle between the jet and crossflow. 

A change in angle simply involves a change in the boundary conditions for the differen- 
tial equations. 


The differential equations (3) through (5) may be solved in closed form in the 
special case when 9 ~ 0 and the velocity ratio is not too large. For, in this case 
equation (1) may then be written 


E = 

p E,Ud 

(6) 

and equations (3) through (5) reduce to 




< A i U j> ' E 

(7) 


> 

II 

o 

(8) 


1 


References 


3 and 4 . 
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( 9 ) 


2 

'Vl = (E+ I C D-> Ud > U 


If we now assume that is proportional to cr and write 


A. = Kd 2 
J 


( 10 ) 


K being a constant, then equations (7) through (9) may be integrated. 

For the case in which the jet is exhausting normally into the crossflow the 
following expression follows for the jet centerline 



(2E X + C D ) K 


2E, 



— i 2 -I 

K md n x 


( 11 ) 


This expression has the interesting property that it implies that x/d 0 is a function of 
z/md 0 . This functional relationship is verified in Table 1 using data obtained from 
Reference 3. 

This functional relationship has been shown to be valid for large values of m in 
Reference 5 where the formula 


(x/d Q ) 2 = 7 (z/md Q ) 5 


( 12 ) 


was obtained for the jet centerline. 
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TABLE 1. TABULATION OF EXPERIMENTAL DATA (REFERENCE 3) 


1 ” 

rn 

x/d o 

4.3 

6.2 

8.1 

z/md 

o 

1 

0.66 

0.65 

0.67 

1.5 

0.79 

0. 79 

0.80 

2 

0.89 

0.88 

0.91 

2.5 

0.98 

0.99 

0.99 

3 

1.06 

1.06 

1.07 

4 

1.17 

1.18 

1.19 

6 

1.36 

1.39 

1.38 

10 

1.63 

1. 64 

1.60 


JET INDUCED VELOCITY FIELD 


Jet induced interference forces and moments can be significant in a V/STOL 
concept and it is desirable to be able to estimate these effects. Two methods will be 
discussed in this paper. The first is a vorticity model which replaces the jet by a 
series of horseshoe vortices. The strength of the vortices are determined from the 
jet momentum flow and the curvature of the jet centerline. 

The second model is a sink-doublet model. The sinks represent the entrainment 
of the crossflow fluid by the jet and the doublet represents the blockage effect of the 
jet on the mainstream. 


VORTICITY MODEL 


When the iet leaves the nozzle it is deflected into the crossflow' direction and 
its cross section changes from its initial circular section into a kidney-type shape as 
illustrated in Figure 1. Two contra-rotating vortices are observed in tests involving 
jets in a crossflow. These vortices are observed to grow in strength as the jet is 
deflected further into the crossflow direction and it would appear that they ought to 
make an important contribution to the jet induced velocity field. 

Consequently, in the vorticity model presented in Reference 6, the jet is replaced 
by a system of horseshoe vortices (see Figure 3). An empirical equation 


o 2 

x/d D = . 19 m cosh (z/. 19 m d 0 ) 


- 1 


(13) 
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is used for the jet centerline. It is assumed that: 

2 

(1) The jet momentum flow p ATL is constant along the jet path. 

(2) The force on a jet element of length <5 s, normal to the centerline, is equal 

2 

to p AjU j 6 s/R where R is the local radius of curvature derived from 
equation (13). 


(3) The mainstream flow near the jet is assumed to be along the jet path 

direction and to have a mean velocity equal to the undisturbed mainstream 
velocity U. 


With these assumptions the circulation, r, per unit length of the jet is given by 


r 

U 


d /U. 

_ 7T Of ]Q 

4 R U 


(14) 


Together with this jet "bound" vorticity there is a trailing vorticity which is a 
consequence of this "bound" vorticity. This trailing vorticity, which increases in 
strength as it moves from the jet exit, accounts for the eventual breaking up of the jet 
into a vortex pair discussed earlier. 


The velocity field due to the interference of the jet on the mainstream is obtained 
by dividing the jet into a number of elements, the vorticity of each element being cal- 
culated from equation (14). Each element is replaced by a horseshoe vortex, the 
trailing arms of the horseshoe vortex being one jet diameter apart. 


Calculations of the velocity field enable the jet induced velocity field on any 
neighboring surface to be obtained. Using a lifting surface theory for a finite wing 
and an image system for an infinite flat plate enables wing loadings and pressure 
coefficients to be derived. Figure 4 shows calculations of the pressure coefficient on 
a flat plate.-*- Calculations of the interference lift and moment for a delta wing are 
shown in Figures 5 and 6. 


The vorticity model is very easy to use and is useful for obtaining the significance 
of the parameters involved. For example, Figure 7 shows how m and x^/c affect the 

interference lift and moment for a 60 degree delta wing. Figure 8 shows the impor- 
tance of the area ratio S/A^. This figure shows the variation of the interference lift 

with area ratio for a given jet static thrust and free stream velocity when either the jet 
area or wing area is kept constant. 


Reference 


7. 
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SINK-DOUBLET MODEL 


The sink-doublet model is a natural extension of the method for obtaining the jet 
centerline discussed earlier. The interference effects of the jet are accounted for by 
distributing sinks and doublets along the jet, as illustrated in Figure 9. The sink 
strengths m*are made proportional to the mass of air entrained by the jet. That is 


The doublets represent the blockage effect of the jet, and their strengths, n, are 
obtained from the coefficient of the l/y term in the complex velocity potential expansion 
W(^) for the two-dimensional flow past an ellipse. 

Thus 


- ”Ud 2 _ 5 z 

^ 4 4 md 0 


(16) 


in the region where the jet is deforming from a circular to an elliptic section, and 





(17) 


after the jet has achieved its final elliptic section. 

The jet induced velocity field can be calculated and wing loadings can be 
obtained by making use of a lifting surface theory. Calculations of wing loadings using 
this model are shown in Figure 10 and show good correlation with test data. 

The normal procedure for calculating surface pressure distributions has to be 
modified in the case of a wing containing a jet exhausting at an angle into a mainstream. 
The reason for this is that besides inducing a velocity normal to the wing, the jet also 
induces a velocity tangential to the wing surface which may be of the same order of 
magnitude as the induced velocity due to the lifting effect of the wing. Thus, although 
this jet induced tangential velocity has only a minor effect on the wing loading, it may 
cause a significant change in the surface pressure field. That there is a significant 
effect is illustrated in Figure 11. Figure 11 demonstrates that the loading is not dis- 
tributed equally on the upper and lower surfaces, which would be the case without the 
effect of the jet. 

Figure 12 shows the theoretical and experimental interference spanwise loading 
for the wing of Figure 10. The theoretical predictions are quite good although it is 
noticed that the tendency is to smooth out the loading which suggests that an improve- 
ment in the numerical procedure employed in the lifting surface program would im- 
prove the correlation. 
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CONCLUDING REMARKS 


An analytical method has been developed for predicting jet centerlines for jets 
exhausting at right angles into a subsonic crossflow. The method calculates center- 
lines which are in good agreement with wind-tunnel test results. This method can be 
extended to include cases when the jet is exhausting into a crossflow at angles other 
than 90 degrees. 


A vorticity method for calculating the jet interference velocity field has been 
discussed. Calculations for the lift and pitching moment demonstrate the usefulness 
oi the method for undertaking a parametric study. 


An alternative model which replaces the jet by a distribution of sinks and doublets 
is also discussed. Calculations of wing surface pressure distributions for a single jet 
exhausting at right angles into a crossflow indicate that this approach will also be a 
useful method for design purposes. This second method has the advantage that it can be 
e. ended to include jet deflection angles other than 90 degrees, multiple iets, changes 
m power and changes in crossflow velocity. 
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Figure 5.- Calculations of interference lift for two-jet delta-wing 

configuration. 



Figure 6.- Calculations of interference moment for configuration 

in Figure 5. 









Figure 7.- Effect of velocity ratio and jet position on interference 
lift and moment for a 60 degree delta wing. 
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Figure 8.- Effect of area ratio on interference lift. 
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Figure 9.- Distribution of sinks and doublets along jet path. 



SPAN WISE 
STATION 


EXPERIMENT 

m JET VELOCITY TUNNEL VELOCITY 
10 1117 FT/SEC 111.7 FT/SEC 

10 430 FT/SEC 43 FT/SEC 

THEORY 


SYMBOL 

A 

□ 



117 











NUMERICAL TREATMENT OF LINE SINGULARITIES FOR 
MODELLING A JET IN A LOW-SPEED CROSS FLOW* 


J. G. Ski fs tad 

School of Mechanical Engineering 
Purdue University 
Lafayette, Indiana 


SUMMARY 


A simple analytical model for the jet in a low-speed cross flow, comprising 
two contra-rotating line vortices and a line sink, is described. The numerical 
formulation of the nonlinear problem and an iterative procedure for its solu- 
tion are discussed. Means for achieving a more complete analysis by modifying 
the simple model are considered. 


INTRODUCTION 


The aerodynamic design of VTOL aircraft would benefit greatly from improved 
means for determining the aerodynamic influence of the lift jets under widely 
varying conditions. Reliable theories would offer highly useful means for that 
purpose and would serve to complement the experimental data. Unfortunately, 
the theoretical problem is essentially nonlinear. Aside from that inconve - 
nience, the problem cannot be accurately formulated in a strict sense because 
the turbulent transport properties for three-dimensional free shear layers may 
only be roughly estimated at best. Experience with semi -empirical theories has 
shown that reasonable and occasionally accurate results may be obtained with 
rather rough analytical models of the jet (ref. 1). It may be anticipated, 
then, that approximate theoretical models of the jet interaction may be suffi- 
ciently accurate and reliable for some aerodynamic purposes. The principal 
concern in this context is the characterization of the aerodynamic effects of 
the jet rather than the structure of the jet proper. 


8 


The object of the study discussed in this paper was to formulate and ex- 
amine the numerical properties of an analytical model of the jet interaction 
problem, subject to the conditions: 


*This study comprises part of the activity sponsored by the NASA Langley 
Research Center under Grant No. NGR 15-005-094, initiated in February 1969. 
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(a) the model should include the principal dynamical aspects of the 
coupling between the jet and the aerodynamic field, 

(b) the initial model should be as simple as possible, but should lend 
itself to more complete representations without major changes in the 
basic formulation, and 

(c) there should be the capability to extend the analysis to include 
multiple jets without major changes in the formulation. 

As indicated in item (b), the initial model was chosen to be as simple as 
possible, while retaining the dynamical aspects of the problem judged to be 
of primary significance. The emphasis was to examine the numerical behavior 
of nonlinear iterative schemes applied to the solution without the burden of a 
more sophisticated representation of the jet. 


SYMBOLS 


A cross-sectional area of the jet 

D diameter of the jet nozzle 

F force acting on the jet 

h 1 1” x r| 

I ^ , 1 2 factors in the integrations for the induced velocity 

M local integrated momentum flux of the jet 

m local mass flow rate of the jet 

rf unit vector normal to the jet trajectory 

R radius of curvature of the jet trajectory 

r,r vector between a field point and a source point and its magnitude, 

respectively (see Fig. 3) 

s distance along the jet 

T unit vector tangent to the line singularity 

U free-stream velocity 

V local velocity in the jet 

V vector velocity at a point on one of the line singularities 
x,y,z Cartesian coordinates of the line singularities 
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Greek Symbols 


e u/v 0 

r circulation of the line vortex 

k factor in the expression for the flow area of the jet 

p mass density of air 

o strength of the sink singularity per unit length 

Subscri pts 

n normal to the jet trajectory 

o evaluated at the nozzle exit plane 

s associated with the sink singularity 

v associated with the line vortex 

x,y,z component in the indicated coordinate direction 

ANALYTICAL MODEL OF THE JET 


The simplified analytical model of a single jet chosen for the initial 
computations comprised a pair of contra-rotating streamwise line vortices and 
a line sink positioned midway between the line vortices. Figure 1 illustrates 
the arrangement of the singularities and the coordinate system for the computa- 
tions. For a single jet as shown in Fig. 1, the two segments of the curves, 
x,y > 0 for the sink, and x,y,z > 0 for the line vortex, are regarded as inde- 
pendent, the other elements of the curves being related to them by symmetry. 

The strengths of the line vortices were presumed constant and known. The 
ratio B = U/V was also assumed to be given. The positions of the curves in 
space and the°distribution of the sink strength along the central curve were 
to be determined in the solution. 

Continuity within the jet flow was simply treated in this initial model 
by relating the local flow rate in the jet to an area arbitrarily assumed to be 
proportional to the square of the spacing between the line vortices. That is, 
a certain similarity of the cross-sectional shape of the jet at different 
streamwise positions was assumed (see ref. 2). Thus, if 

A/A q = 1 + k(z 2 /z 2 - 1) 
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then, since 


nt/m ° °{k ” = 1 + 166 ^ ™ d(s/D) 

one has finally 



The entrainment of mass per unit length is given by 4irpo. 

The line vortices were regarded as free streamlines in the tangent plane 
of the jet; that is, they were considered free to establish equilibrium posi- 
tions in the z-di recti on, such that 

(v x t) z = 0 (2) 


along each curve. 

The integrated streamwise momentum flux in the jet was approximately 
treated by assuming it to be related to the entrainment of mass in the jet 
according to the expression 


dM 

d? 


Ut 


dm 
x d? 


The latter may be rewritten in integrated form as 


M. 


= 1 + 


s/D 
166 2 I 
o 


-2-t 
DU x 


d ( s/D ) 


(3) 


Finally, the force per unit length of the jet acting normal to it was re- 
lated to the integrated streamwise momentum flux of the jet and its radius of 
curvature. 

f n ‘ S (4) 

The force F was evaluated as the sum of the forces per unit length acting on 
the three line singularities at any given position along the jet. For the 
sink singularity (ref. 3) one has 

F ns = 4TTpov*n (5) 
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and for the line vortices (ref. 3) the expression is 


F nv = 2 P r(v x t)-n (6) 

Combining Eqs. (4), (5), and (6), and solving for R/D yields 


R 

D 


= (M/MJ/168 2 


>/' 



• n 


(7) 


The x and y coordinates of the jet trajectory may be expressed (ref. 4) in 
terms of R as 


s s 

x = / [sin / ds'/R(s')] ds 

o o 

s s 

y * / [cos £ ds ' /R(s ' )] ds 


( 8 ) 


Equations (1), (2), (3), (7), and (8) are to be satisfied by x,y,z, and o at 
each position s along the jet. The velocity v in Eqs. (2) and (7) may be 
expressed (ref. 5 ) in terms of the geometry and a as 


v 


17 j. r / dsxr 

U + 47 J r 3 

line 

vorti ces 


-/ 

sink 

line 


(9) 


The analytical model as described above is not considered to be complete 
in the sense that a number of unnecessary, but expedient approximations are 
involved. As noted in the introduction, the object of this initial study was 
to examine the numerical properties of the problem, and the model described 
did include what were considered to be the principal factors having a bearing 
on the numerical problem. A more exact representation of the problem may be 
prepared without introducing significant changes in the numerical scheme. 
Furthermore, the model could rather simply be extended to include several jets. 

The aerodynamic shape of the jet provided even by this simplified model is 
notunlike that observed. For example, the streamlines in a plane cutting the 
jet near its origin appear as shown schematically in Fig. 2. The wake of the 
jet is, of course, entirely neglected by the model in its initial form. 

It may be noted that if r/DU in Eq. (7) is a constant, as might be anti- 
cipated, the formulation of the problem in non-dimensional terms, indicated in 
the format of the equations, yields solutions with the single parameter g. 
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NUMERICAL FORMULATION 


While the line singularities extend to infinity in the analytical model, 
the numerical treatment requires approximating the lines by finite line ele- 
ments, introducing some degree of error. The computational scheme represented 
each of the two independent curves (the sink and one of the vortices) by 40 
points equally spaced along the line. A straight line element tangent to the 
end of each line was found adequate as an "end correction," compensating for 
the truncation of the theoretically infinite line. The computational scheme 
adopted proceeded as follows. 


Initial Trajectory of the Jet 

A trial set of points representing the two curves was first chosen to 
start the computations. Each curve was represented by an even number of se- 
quentially numbered points, not counting the point on each curve in the x-y 
plane. Local interpolation parabolas were determined for sets of 3 points 
centered at each of the odd-numbered points along a given line singularity. 

The length of each segment was then computed, the total length of each curve 
was determined by summing the lengths of the segments, and the increment as 
required for equally spaced points was computed. A new set of points, equally 
spaced along the initial trial curve, was finally determined as starting data. 

Smoothing has been found to be a necessary part of the computational 
scheme to prevent the development of unstable irregularities in the curves, 
particularly following an iterative change in the locus of a curve, but also in 
the initial trial data. A seven-point cubic smoothing routine (ref. 6) was 
employed to smooth the values. The seven-point cubic routine has been found to 
work well with the functions considered; higher difference smoothing has been 
found to be too severe in some circumstances. 

Following the smoothing operation, the trial point sets were re-entered 
as initial data, and the above procedures were repeated, except for the 
smoothing. In addition, the local length of each curve at all points, together 
with the components of the tangent and normal vectors, were computed. 

A trial distribution for a (s) was chosen along with the initial coordinate 
data. Those values were also conditioned and smoothed, as discussed for the 
point sets. 


Computation of the Velocity Induced at Points Along Each Singularity 

The velocity induced at a point on one of the line singularities may be 
determined according to Eq. (9), which involves summing the line integrals over 
all of the singularities present, including, of course, the images of the curves 
considered in the initial trial point sets as shown in Fig. 1. For a given 
"field" point, the integrals were treated numerically by approximating a line 
element about a "source" point by a straight line segment of length as, tangent 
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to the singularity curve, as illustrated in Fig. 3, and summing over all of the 
points on the singularities. For a vortex line segment, one has 


where 


s+As/2 

/ 

S-AS/2 


ds x r _ t x r T 
— r* *1 


h = | t x r| , and 


I 


1 


t* r + As/2 

[h 2 + (t-r + as/ 2 ) 2 f 2 


t»r - as/2 

[h 2 + (t-r - As/2) 2 ]* 5 


For a sink singularity, the expression is 


s+As/2 _ oI.(t x r) x t _ 

-of rds/r 3 = - + a^t 

s-as/ 2 h 2 


where I 2 = t— — r - jzrz r 

1 1* r + as/2 I 1 1* r - as/2 | 

In the computational arrangement, the contributions of a line segment about a 
"source" point and those of its images were computed simultaneously employing 
symmetry relations for the local properties of the curves. 


Iteration Procedures for Adjusting the Positions 
of the Curves and the Sink Strength Distribution 

Before proceeding with the iterations, the free-stream velocity, U, was 
adjusted slightly to ensure that the line vortices at their intersection with 
the x-z plane remained stationary. That option was chosen for the initial 
computations over the alternative of changing the initial positions of the vor- 
tices at that point. The change required in U was generally rather small and 
did not tend to affect the iteration procedure adversely. The x-components of 
all computed velocities were adjusted accordingly. 

The first step in the iteration procedure was to ?llow the line vortices 
to assume their equilibrium positions in the z-direction, corresponding to the 
trial coordinates (x,y) of the jet trajectory. That is, the condition given 
by Eq. (2) was to be satisfied at all points on the line vortex. For rough 
computations it was found that the approximation obtained by taking the equili- 
brium separation of a pair of infinite, parallel, rectilinear vortices tangent 
to the jet at each point was quite reliable. For further refinement, changes 
in the z-coordinate were taken proportional to (v x t) z . The predicted cor- 
rections were smoothed using the 7-point cubic smoothing routine. Between 
iterations the adjusted point sets were conditioned as described for the ini- 
tial data, but without further smoothing. 
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When the condition (2) was satisfied to within an acceptable error, the 
variation of the integrated streamwise momentum flux in the jet was determined 
by numerically integrating Eq. (3). Given those values, the continuity ex- 
pression for the flow in the jet, Eq. (1), yielded the variation of a with 
distance along the jet. A differentiation expression based on five-point 
Lagrangian interpolation (ref. 6) was employed for the computations. 

Following the determination of the new a values, the entire procedure was 
repeated until the condition that the vortex lines assume their equilibrium 
positions was again satisfied, Eq. (2). That normally required no further iter- 
ation because their positions in the z-direction were rather insensitive to 
the sink strength. 

Finally, Eq. (7) was employed to compute the distribution of R along the 
jet. Through Eqs. (8) a projected trajectory for the jet could be determined 
corresponding to the computed R variation. Unfortunately, no straightforward 
or reliable iteration scheme for adjusting the position (x,y) of the jet tra- 
jectory could be established without some experience with the behavior of the 
functions involved. The evidence indicated that a conventional Newton-Raphson 
technique would not be likely to be successful in projecting the coordinate 
positions. Rather, it was determined that efforts should be concentrated on 
changing the R(s) variation iteratively. While a completely satisfactory itera- 
tion scheme was not established, a refinement of the following scheme appeared 
likely to be acceptable. The technique adjusted the local radius of curvature 
proportional to the negative of the difference between the predicted value and 
the value for the trial curve. That is, if the predicted local radius of curva- 
ture was larger than that of the trial curve, the local value was lowered in the 
succeeding iteration. Sample computations following such a procedure are il- 
lustrated in Fig. 4. Figure 4b also shows the variation of o and the position 
of the line vortex corresponding to iteration (B) of Fig. 4a. These results 
merely show the trends in the computations and do not necessarily reflect the 
final values of the solution. It might be mentioned that trajectories have 
been computed appearing to agree much more closely with the trial values, but 
differing more in the local R(s) variation than those shown in Fig. 4. 


DISCUSSION 

The numerical formulation of the problem adopted appears to be accurate 
and acceptable. Further alterations of the iteration procedure may be desirable, 
although the starting curves chosen for the computations were not close to the 
solution positions, and better Initial estimates should improve the performance 
of the Iteration technique. 

Extensions of the analysis to yield a more complete description of the 
problem center on improving the characterization of the jet. Given the aero- 
dynamic field around the jet from the simple model, it is possible to more 
accurately determine the entrained momentum flux, pressure variations around 
the surface of the jet, and so forth- The effects of the wake of the jet 
probably should also be added (ref. 1), accounting for a shedding of streamwise 
vortlcity from the jet and changes in the strengths of the vortices associated 
with the jet proper. 
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Further computations should indicate what degree of approximation of the 
problem is necessary for aerodynamic purposes. 
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ANALYTIC DESCRIPTION OF JET-WAKE CROSS SECTIONS FOR 


A JET NORMAL TO A SUBSONIC FREE STREAM 
By Richard J. Margason 
NASA Langley Research Center 


SUMMARY 

An analytic description is given of jet-wake cross sections for a jet 
exiting normal to the free stream. The results indicate that a reasonable 
description can he obtained for the cross section. However, a conplete descrip- 
tion of the jet wake and associated phenomena requires a more complete model of 
the jet wake than that presented herein. 


INTRODUCTION 


The analytic description of a jet in a cross flow is difficult because of 
the complex nature of the flow in the vicinity of the jet. The primary fea- 
tures of the jet are the solid fluid core near the exit, the rollup by the wake 
into a pair of trailing vortices, and entrainment of free- stream flow into the 
jet. This entrainment can be considered in two parts: the fluid induced into 

the wake region by the swirling flow induced by the vortex pair, and the fluid 
entrained by viscous mixing on the periphery of the jet wake. The jet influ- 
ences the free- stream flow in several ways. For example, near the jet exit the 
solid fluid core blocks the free- stream flow and causes a stagnation point ahead 
of the jet. Beyond the stagnation point, the free stream flows around the jet 
and separates from the downstream side of the jet efflux and, thus, causes a 
separated wake region. 


SYMBOLS 


D jet diameter 

N number of vortex filaments used to describe a jet cross section 

rj jet radius 

— ^ 

s path- element vector 

— > 

V velocity vector 

Vj jet velocity 

V free- stream velocity 
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V^/Vj effective velocity ratio 

X,Y Cartesian coordinates 

z complex plane coordinate 

P circulation strength 

r n circulation strength of the nth vortex filament, where n = 1,2,... N 

0 cylindrical coordinate 

1 natural jet coordinate along the axis of the jet 

$ complex velocity potential 

CALCULATION OF THE ROLLUP IN THE JET WAKE 


Information from reference 1 by Chang-Lu was used in an effort to describe 
the rollup of the jet wake into a pair of vortices. The cross section of the 
jet wake at the exit plane is described by filaments of vorticity as shown in 
figure 1. A discrete number of filaments parallel to the direction of the jet 
velocity are spaced around the jet exit. The circulation strength of these 
vortex filaments is defined by the complex potential function which describes 
flow about a circular cylinder perpendicular to streaming flow: 


* = (z + 

The total derivative of this complex potential function yields the velocity 
vector induced by the presence of the cylinder: 

dz 



Then, this derivative is integrated along a path to determine the circulation 
strength of the discrete vortex filaments, which is given as 


r = 




Integration over a sector of the circle gives the circulation strength for a 
single vortex filament as 


r n = 


sm sin 
N 


f(n 
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These discrete vortex filaments are used to determine the change in the 
cross section of the jet as the flow passes downstream. This is done hy a 
series of computations where the cross-section deformation is treated in two 
dimensions. Each of the vortex filaments is influenced for a small increment 
of time hy the other filaments in the cross section. Since these filaments lie 
on a free surface, they must move to assume a new position where the net force 
induced hy the other filaments is zero. This process is repeated many times as 
the flow moves away from the jet exit. As a result, the wake cross section 
changes shape. This change is indicated in figure 2 for a wake with 12 vortex 
filaments. The circular cross section in section A flattens on the downstream 
face to form the cross section labeled section B. Then, farther down the wake 
at section C, the characteristic kidney shape is formed. In section D and 
farther downstream, this evolves toward a very tightly wound pair of vortices. 

To form the surface of the wake boundary, a three-dimensional lattice of vor- 
ticity is constructed around the jet path. The jet path has been obtained from 
an empirical equation (ref. 2). This series of computations represents one 
feature of the jet, the rollup into a pair of vortices. 

The precise form obtained analytically for this deformation is shown in 
figure J, where the cross sections are described by 96 vortex filaments. The 
cross sections are presented at one-diameter increments along the jet axis. 

They represent the rollup of the wake cross section where the effective velocity 
ratio is 0.25. Previous attempts at using this procedure to describe the rollup 
have failed to form the spiral pattern because of truncation error. 

The accuracy of this representation is indicated by figure 4. This photo- 
graph shows the cross section of a jet wake in a water tunnel at a point 
approximately 6 nozzle diameters along the jet axis measured downstream from 
the nozzle exit. 


FLOW INDUCED BY THE JET WAKE 


The usefulness of this three-dimensional vortex- lattice model of the jet 
wake will be determined by using it to estimate the pressure distribution 
induced on a plane flush with the jet exit. A simplified version of this model 
is represented by potential flow around a solid circular cylinder. (See fig. 5.) 
Ahead of the jet, there are positive pressures which are a result of the blockage 
of flow by the clinder; on the sides, there are negative pressures caused by the 
accelerating flow around the cylinder. Since the potential- flow model does not 
represent flow separation, there are positive pressures at the rear of the jet 
exit. 


If this circular cylinder is allowed to deform along the jet path and roll 
up into a pair of vortices as shown earlier, the pressure distribution presented 
in figure 6 is obtained; there is little change from the pressure distribution 
shown in figure 5- Poor correlation was obtained between the data of figure 6 
and those of figure "J, which are experimental results from reference 3* It 
should be noted that the present representation describes the rollup of wake 
caused by the external free stream and the blockage caused by the presence of 
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the jet. Comparison with experimental results shows the present representation 
of the jet wake to be incomplete. 


CONCLUDING REMARKS 


The results presented indicate that a reasonable description of the cross 
section of a jet wake normal to a free stream can be obtained by an analytic 
procedure. However, a complete description of the jet wake and associated flow 
phenomena will require a more complete model of the jet wake than that presented 
herein. 
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Figure 1.- Cross section of jet wake at exit 
plane described by filaments of vorticity. 
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Figure 4.- Photograph of the cross section of a jet wake in a water tunnel exiting 
perpendicular to the free stream at a point approximately 6 nozzle diameters 
along the jet axis downstream from the nozzle exit. (Photograph from O.N.E.R.A. 
Film Wo. 575 entitled "Flows with Large Velocity Fluctuations.") 




Figure 5*- Pressure distributions around a circular cylinder. 



Figure 6.- Pressure distribution induced by a vortex- lattice jet wake with 

changing cross section. 
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CROSS WIND EFFECTS ON TRAJECTORY AND CROSS 

SECTIONS OF TURBULENT JETS 

By Gerhard W. Braun and 
Jack D. McAllister 

University of Tennessee Space Institute 


SUMMARY 


Results in two loosely connected research areas concerning 
jets in cross flow are reported. An analysis of specific mass- 
momentum flux data showed that the entrained mass flux, Q e , 
carries the undisturbed cross flow velocity Voo into the jet in 
such a way that it adds the momentum Q e Voo to the jet flow. This 
relation permits a calculation of the jet trajectory, when the 
mass entrainment is given, or a calculation of the entrainment, 
when the trajectory is known. 

If specific entrainment data are not known, the jet trajectory 
can be calculated by means of the assumption that the sidewise 
spreading rate of the jet and the nondimensionalized distribution 
of the excess velocity are the same as known from submerged jets. 
Closed form solutions are obtained for simple cases. Measured 
trajectories are curve fitted by adjusting a free spreading rate 
parameter. 

The second part deals with the roll-up process of a jet in 
cross flow. Computer results obtained from potential flow cal- 
culations are shown for single and double jets with circular or 
elliptic exit nozzles. 


INTRODUCTION 


The theory of the submerged turbulent jet of circular cross 
section, i.e., the jet issuing into still air is, after Tollmien's 
pioneering work (1), sufficiently well known. The flow pattern 
of the jet in cross flow, however, is much more complex, and it 
is unlikely that, in the near future, its theory would be brought 
up to the level of knowledge on the submerged jet. 

On the other hand there exist empirically derived formulae 
for the jet trajectory in cross flow. Such formulae have been 
published by Heyson (2), Monical (3), Wooler (4), and others. But 
these formulae tell too little about the physical reasons Jor the 
jet curvature. 
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More insight is provided by the jet theories which explain 
the jet bending as the effect of a drag force acting along the 
jet. Abramovich's (5) formula (12.160) is of this type. If 
one analyzes what drag coefficients would be needed, then one 
finds values up to 35 which is about twenty times the value one 
would expect. 

Crowe and Riesebieter (6) recognized this difficulty and 
tried to eliminate it by prescribing a drag coefficient of 1.5 
and making the necessary adjustment by selecting the needed 
frontal area of the jet. The second author of this paper found 
(7) that one needs a jet width of up to three times the width 
one would estimate from experimental observations. 

Because of these difficulties it was tried to derive the 
trajectory equation of a jet in cross flow on the basis of the 
momentum theorem. The problem was solved by the second author 
in his thesis (7) to which the bulk of this paper is devoted. 

As an appendix some potential flow aspects of the roll-up 
process of jets are added. 

For the sake of clarity it is here mentioned that the fol- 
lowing analysis is confined to the similarity region of the jet 
(Fig. 1), where the unmixed ("potential") core has disappeared 
and the transition to pure free turbulence has been established. 
The effects of the jet in front of this similarity region are 
taken care of by introducing an apparent origin of the jet. 


SYMBOLS 


a 



K(R) = || 


M 


location of the apparent jet source measured in orifice diameters 
upstream of the orifice 

jet exit area 

area of the jet projection onto the x-y plane, see figures 7 an( f 8 

diameter at jet exit 

width of the jet in y direction 

ratio of jet mass flux at potential core end to value at jet exit 
spread rate of the jet width 


Q 


r 


142 


jet momentum flux 
jet mass flux 

polar coordinate in jet cross plane 





5 

n 

9 

I 

p 

q> 


characteristic ratio for jet in cross flow defined by equation (l8) 

local flow velocity 

_ c)M 

dQ 

Cartesian coordinate system, see figure 7 

_ j 8 K(R) see e q ua tion ( 25 ) 
yn sin 0 Q 

nondimensional radial distance, r/rx 

angle between local jet velocity and free stream velocity 

distance along curved center line of jet 

density 

angle of polar coordinates in cross-sectional plane 


Subscripts: 

value at center line of jet 
j pertaining to the jet 

n center of jet 

o value at jet exit 

r reference station at end of "potential core" 

x component parallel to x-axis 

z component parallel to z-axis 

§ value corresponding to a velocity excess one-half the maximum value 

oo free stream value 


MOMENTUM THEORY OF JET TRAJECTORY 


Analysis of Mass-Momentum Flux Data 

The first goal was to establish assumptions about the jet 
momentum flux and the mechanism of interaction between the jet 
and the external flow. 

For this purpose mass-momentum flux data published by Keffer 
and Baines (ref. 8) were analyzed. These values are reproduced 
in figure 2. They were obtained by integration of measured 
velocity profiles using as the jet boundary the surface where the 
velocity excess is one-tenth the maximum value at a given axial 
distance . 
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Figure 2 illustrates that the momentum flux increases lin- 
early with mass flux for each jet to free stream velocity ratio. 
Thus a constant velocity V c = dM/dQ characterizes for each ve- 
locity ratio the axial momentum per unit mass being entrained in- 
to the jet. To obtain quantitative values for this characteristic 
velocity, V c , the following assumptions were introduced. 

1. The linear relation between mass and momentum fluxes 
is also valid for the initial mixing region between 
the orifice and the end of the potential core. 

2. The reference mass flux values occurring at the potential 
core end of a jet in cross flow may be reasonably es- 
timated from data measured for the submerged jet. 

3. The velocity is constant over the orifice cross section. 
(Confirmed by private communication from Professor 

Kef f er . ) 

4. The jet and external stream have the same density. 


The mass-momentum flux data of fig. 2 may be expressed in 
the form of equation (1) with F representing the experimental 
slope of the linear relationship. 


M/M r = 1.0 + F(Q/Q r - 1.0) . 


( 1 ) 


From fig. 2 one reads the values for F: 


R 

4 

6 

8 

F 

.412 

.304 | 

.269 


Assumption (1) provides a second expression between the mass and 
momentum fluxes: 


M = M 0 + V C (Q - Q q ). 


( 2 ) 


A relation for the characteristic velocity, V c , is then obtained 
by combining equations (1) and (2) 


M r 

V = F — 
c Q 

^r 


(3) 


The desired expression for the characteristic velocity is found 
in the form of equation (6) after introducing the abbreviations 


k = Q /Q 
r v « r / 


v j E M <A> ’ v j s R v < 


(4) 

(5) 


R 


V„/V_ = 


k r /F - (k r - 1) 


( 6 ) 


144 


According to assumption (2) the mass flow ratio, k r/ does not de- 
pend upon V m so that for it the value k„ =2.2 based upon measure- 
ments for a submerged jet (given by Albertson, et al. (ref. 9)) 
can be substituted. Equation (6) then gives the values 0.97, 

0.99, and 1.14 for the ratio of the characteristic to free stream 
velocity . 

This result suggests that, to the first approximation, mass 
entrained into the jet brings with it the associated free stream 
momentum. Further it suggests that deflection of the total jet 
momentum flux is more adequately treated as an effect of momentum 
entrainment than as the result of pressure forces normal to the 
jet . 


Fearn took the relation V c = for granted and 
used it already in early 1968. 


Approximate Integral Analysis 


In fig. 3 the nondimensional velocity excess above the free 
stream value published in ref. 8 is plotted versus the non- 
dimensional lateral coordinate. The nondimensionalizing distance 
ri is the radius at which the velocity excess is half the maximum 
value. The curve is universal, i.e., does not depend on R or £ . 
Schlichting ’s (ref. 10) expression 


U 


obtained for the 
curve . 



= f (rj) = 


1 

[1 + (V"2 - 1) T ) 2 ] 2 


(7) 


submerged axially symmetric jet fits well this 


The approximation for the jet axial velocity distribution is 
completed by using measured values of the decay in center line 
velocity a s published in ref. 8 together with a suitable ex- 
pression for the angular and axial variations in ri. The latter 
was determined by approximating the jet cross section with an 
ellipse with the eccentricity k and the spread rate a, giving 



cos 


<p. 


( 8 ) 


The factors k and a have to be determined experimentally. 


The obtained axial variations in mass and momentum flux are 
presented as dashed lines in figs. 4 and 5 for velocity ratios of 
4 and 8. Note that the mass entrainment is materially increased 
when the free stream velocity becomes comparable with the jet 
velocity. The result for the submerged jet was obtained through 
the limit process of R-*» with k->0 and a taken from submerged jet 
data of Forthmann (ref. 11). The limit value is in excellent 
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agreement with the direct measurements of Ricou and Spalding 
(ref. 12) and data analysis of Albertson (ref. 9), both works in- 
dicating Q/Q 0 =0.32 i/d Q for submerged jet. 

Applying the same limit process to the momentum flux ex- 
pression with the conventional assumption of constant momentum 
flux when the external medium is at rest yields the correct ex- 
pression for the center line velocity decay even to the numerical 
constant; U _ y 

G(|) = — — = — ' — for submerged jet. (9) 

v j - ? /d o 

The analysis of the following section will show to what ex- 
tent the approximate integral analysis is in agreement with the 
principle that the entrained mass carries its free stream momentum 
into the jet. 


Relation of Jet Trajectory to Mass Flux 

The principle of momentum entrainment may also be applied to 
link the jet center line inclination to the mass flux at a specif- 
ic axial distance* . The analysis is again based upon three 
assumptions : 

1. The jet trajectory has the direction of the jet momentum 
flux vector. 

2. The momentum flux at a station (* of the trajectory is 
the vector sum of the momentum flux through the orifice 
and the momentum of fluid entrained along the trajec- 
tory up to station £. 

3. The entrained fluid momentum is the product of the un- 
disturbed cross flow velocity and the entrained mass. 


Using the coordinate system of fig. 6, the momentum flux 
has the components: 


M x - M o sin 0 O’ 

M z = M 0 COS 9 0 + V oo (Q " V • 


( 10 ) 

(ID 


Hence the following relationship exists between the trajec- 
tory angle 0 and the mass flux: 


M V Q 

7 00 O 

cot 0 = — = cot 0 + 

M ° M 

X 


sin 0^ 
o o 



- 1 


( 12 ) 


* Fearn also used this principle in ref . 13 
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(13) 


Solving for the mass flux ratio gives 

q/q = 1 + R sin 0 (cot 9 - cot 0 ) , 

^ o o 

where R = V ./V . 

J 00 

Mass flux values based upon equation 13 and jet trajectory 
data of ref. 8 are compared in fig. 4 with the corresponding 
results of integral analysis. The good agreement not only sub- 
stantiates the momentum entrainment hypothesis but also suggests 
that mass entrainment for a jet in cross flow may be readily 
obtained from center line trajectory data via equation (13). 

A simple relationship for the momentum flux may also be ob- 
tained by using the fact that the x component is constant 

sin 9 

M/M = - . (14) 

sin 9 

Equation (14) is compared with momentum flux results based 
upon the integral analysis in fig. 5. 

The comparison is satisfactory especially at the more down- 
stream positions where the axial velocity should be quite close 
to the assumed similarity profile. 


Development of the Trajectory Equations 

At first thought the back-flow produced by the vortex street, 
observed in the wake of the jet, seems to eliminate the possi- 
bility that the entrained masses carry in the average the full 
cross flow velocity. Platten and Keffer (ref. 14) developed 
about simultaneously with the reported work a jet entrainment 
theory which has the contribution of the shed vortices as a 
specific feature. The following pages will show that the as- 
sumption of a wake with zero velocity can serve the same purpose. 

The analysis of a circular jet entering at an angle 0 O to a 
uniform stream is initiated by considering the balance of mom- 
entum flux for the control volume presented in fig. 7. Ag is 
part of a rigid wall through which the jet exhausts. The shear 
force acting on A 3 is assumed to be negligible in comparison to 
the momentum terms. As is the surface normal to the jet trajec- 
tory at the general axial distance, The momentum flux in the 

x and z directions are established under the following conditions: 

1. The velocity vectors in each jet cross section are 
parallel to the jet trajectory. 

2. Average pressure differences between corresponding sur- 
faces of the control volume are relatively small and 


147 


may be neglected in comparison to the momentum terms. 

3. The flow through the plane A* can be divided in three 
parts: ^ 

a) A jet region, A., bounded by the jet-external 
stream interface. 

b) A wake region in the "shadow of the jet” with 
essentially zero velocity. 

c) An exterior region comprising the remainder of At 
and having the undisturbed free stream velocity. 

Using these assumptions the x momentum equation is 

M q sin 9 q = sin 6 J pU 2 dA, (15) 

A . 

J 

which is a mathematical statement that the component of jet mo- 
mentum normal to the free stream is invariant. 


The z momentum equation has the form 

+ » 0 cos e 0 = cos e / p u2<1a + ( Wp „ v2 <16) 

A j 

where A is the projection of the jet onto the x-y plane, see 
fig. 7. p 


This expression is further simplified to equation (17) by 
using equation (15) and slight rearrangement. 

p ^ 

cot e = — = cot e + — 2L-S2 — E_ . ( 17 ) 

dx ° M sin 9 

o o 


Equation (17) suggests that the effective velocity ratio be de- 
fined in terms of the orifice momentum flux as given by 




A 

o 


A p V 

O r oo oo 


M 

o 

A p V 2 

o 1 OO oo 


(18) 


Equation (17) is thus expressible in terms of the jet geometry 
and the velocity ratio, R: 


- = cot 0 + ^ -5 . (19) 

dx A sin 9 

o o 
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This trajectory equation may be solved as soon as A p is expressed 
in terms of x or z . 

By generalizing Vakhlamov 's (ref. 15), approach, the projected 
jet cross section A p (fig. 8) was built-up of a parabolically 
flaring out stem of width, h, and a head in form of half an 
ellipse. The effective jet width is taken to be a linear function 
of axial distance with an empirical proportionality parameter, 

K(R) called the spreading factor and a set-back factor, a, which 
corrects for the fact that the apparent origin is located a 
orifice diameters upstream of the orifice. This gives 

h = (2ad 0 + 4 ) K(R) (20) 

Evaluation of tests gave for a jet entering normal to the 
free stream 

— = K(R) at 1.0, for 4 < R < 8, (21) 

and lim 3h ^ 0 32 (22) 

R+~ ^ - 

for a submerged jet. 


Solution of the Trajectory Equation 

There is no difficulty in solving the trajectory equation 
(19) on a computer. For basic investigations, however, an 
analytical solution is desired. 

In order to make such solutions easier, the elliptic head of 

x 

the jet shadow A p is neglected so that Ap —J h dx . One has 
then to substitute in equation (19) 0 

= 1 JilSi f (2 ad + ?) dx (23) 

A ir < 0 ° 

o o 


— | 2 dx (24) 

dx/ 

Equation (24) causes certain analytical difficulties. However, 
if one plots the nondimensionalized penetration distance x/d Q R 
versus the nondimensionalized trajectory arc length |/d 0 R then 
one sees that a good approximation of the arc length can be 
obtained by assuming 


where 


4 = 


AA * 


4 


x 


for 


— — < 1 . Region A 

d R ' 
o 
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i = 2x for > 1 Region B 

V 

Integration gives in Region A the trajectory equation 

t — 

.2 


Z = X cot 0+5' 
o 


where 


Z = 


Rd 


X = 


Rd 


— X 2 + -L X 3 
2R 

; s 2 = 


12 
8 K (R) 
7T sin 0 


3tt 

For the case a = o and K = — equation (25) is identical with 
Heyson’s <ref . 2) empirically found formula 


( 25 ) 


Z = X cot Q + 


4 sin 0 


(26) 


For Region B, i . e . , X > X-^ one finds 


Z = (Z 1 + -) cosh[5(X-X 1 )] + 5 -1 z' sinh[ 5(X-X 1 ) ]-a/R 
R 


(27) 


Z^ and are the nondimensionalized z-coordinate and its 
derivative at the matching point of Region A and Region B. 

If one would consider a jet issuing perpendicular to the 
cross flow and if one would expand the validity of ( 27 ) down to 
the jet orifice, i.e., if one would assume X^ = Z-. = Zi = 0, 
then (27) would be simplified to 

Z = - [cosh (5X) -1] (28) 

R 


This formula is similar to an empirical relation given by Wooler 
in ref. 4. 


Correlation of Computed Trajectories with Measurements 

The primary substantiation for the simplifying assumptions 
used in the development of the trajectory equations is the abil- 
ity to match measured trajectories over a wide range of R using 
values of K(R) and the apparent origin set-back distance, a 
that are in agreement with experimental observations. The re- 
sults of such matchings are presented in this Chapter using 
previously published data (16), (l7), (18), as well as measure- 

ments obtained during the present investigation. 

It was found that for all trajectories the same value a=0.8 
may be used for the apparent origin set-back distance. The 
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coefficient a/2R in equation (25) provides for the observed change 
in trajectory shape from quadratic to cube as R increases. 

Figure 9 illustrates the close degree to which the air jet tra- 
jectories of references 16, 17, and 18 may be matched. Values 

of the spread rate parameter, K(R), used in these computations 
are plotted in fig. 10. It is important to note the close agree- 
ment between the values required to match measured trajectories 
and the independent experimental values based upon a ten per- 
cent velocity excess (equations (21) and (22)). This is a 
strong substantiation for the assumptions used in developing the 
trajectory equation. 

The original thesis presents also matched water jet trajec- 
tories which show the same good agreement. The thesis further- 
more investigates the effect of a plate from which the jet issues. 

Finally, the effect of having the jet enter at angles other 
than normal to the external stream was investigated. Trajec- 
tories measured for entry about thirty degrees each way from the 
normal are shown in fig. 11. 

The numerical predictions of fig. 11 were accomplished with 
the K(R) and a used to match for the normal jet entry. Note 
that, as 0 O decreases, the jet turns more slowly toward the 
free stream indicating decreased mixing and hence the desirability 
of a smaller value for K(R). However, for 0 Q between sixty and 
one hundred twenty degrees the predictions presented are con- 
sidered to be quite satisfactory. 


INVESTIGATION OF JET ROLL-UP DUE TO CROSS FLOW 


A typical and very conspicuous phenomenon of the jet in 
cross flow is the formation of the already mentioned vortex 
system. This system consists of two types of vortices. One 
type stays within the jet and causes it to roll-up into two 
vortices; the other vortex system separates from the jet and forms 
a type of von Karman vortex street. 

Chang (ref. 19) has investigated the roll-up process on the 
basis of two-dimensional potential flow theory. She represents 
the deformation which takes place on a jet slice while it moves 
down the jet by assuming that it would be identical with the de- 
formation a cylindrical jet would undergo in two-dimensional 
cross flow while time goes on. Mrs. Chang’s investigations are 
based on potential flow theory, but her results show such a 
striking similarity to the actual phenomenon that it was decided 
to continue her investigation. This was done by Mr. Stephen M. 
Soukup (ref . 20) , who included in his analysis also jets with 
elliptical outlets and added investigations on couples of jets 
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in side by side and tandem arrangements. 

The analytical approach is as follows: The two-dimensional 

flow around the jet outlet configuration is established by dis- 
tributing a vortex sheet of such strength around the openings 
that it would become a stream line. For computation purposes 
the continuous vortex sheet is replaced by 32 discrete vortices, 
and the velocity components, v x and v , imparted to each vortex 
by the rest of the vortices , are calculated . Selecting a time 
interval At and permitting the vortices to move freely, the dis- 
placement of each vortex in the direction of the two coordinate 
axes is given by ax = At-v x and Ay - At-v The vortex dis- 
placements are different for the different vortices and represent 
the deformation of the jet cross section. Two examples of cross 
section deformations are shown in fig. 12. The upper row is for 
a circular exit nozzle and the lower row for an elliptical exit 
nozzle with the numerical eccentricity E = 0.80 and the same 
area as the circle. 


The calculations for single jets cover elliptical nozzles 
with the numerical eccentricities of 0.0, 0.2, 0.4 0.6 0.8, and 0.95. 
For a more compressed presentation the deformed cross sections 
were stacked as it is shown in figs. 13 and 14. The side by side 
configurations consisted of elliptical jets with the following 
eccentricities, e: 


right jet 

e = 0.00 

0.0 

0.95 

left jet 

e = 0.00 

0.95 

0.95 


Circle and ellipse have the same outlet cross section area of 
Ao = it* The distance between the centers of the two cross 
sections was in all cases 3 units. 


Finally the following two configurations of jets in tandem 
arrangement were investigated: 


Front jet 

circle 

circle 

Rear jet 

circle 

ellipse e = 0.95 


The results of the first case are shown in fig. 14. 


The calculations reveal that the jets of elliptic cross 
section roll-up slower than the circular jets. For side by side 
configurations the cross section deformation is much larger on 
the inside than on the outside of the jets. In tandem arrange- 
ment the rear jet is pulled into the front jet. 

Unfortunately no relation between the time sequence of the 
deformation of the hypothetical two-dimensional jets and the 
deformation of the actual jets in space has been found. It is 
obvious that the potential flow theory cannot describe the mix- 
ing phenomenon of the jet. However, it is hoped, that Soukup's 
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calculations may help in finding the effect of the cross flow on 
the enlargement of the jet surface which obviously is important 
for the entrainment process. 


CONCLUDING REMARKS 


The presented studies lead to the following conclusions: 

1. The deflection of a turbulent jet in cross flow is 
primarily due to a mixing process which entrains free 
stream momentum into the jet. This process permits 
determination of the jet mass flux from the more easily 
measured trajectory angle. 

2. Momentum mixing considerations allow jet trajectory 
predictions over a wide range of velocity ratios, R, 
and jet entry angles, 0 Q . 

3. These favorable results should not raise the hope that 
refinements of this approach, like a more careful 
representation of the jet shadow area, A p , could produce 
an essential improvement of the accuracy of the results. 

4. The momentum theory is confined to cases where flow 
restrictions are of minor importance. If one wants 
to include ground or wall effects one would have to 
introduce new features like vortices which have a 
far-field effect and are subject to a feedback from 
flow obstructions. 

5. A challenging extension is seen in including the vortex 
system especially the shed vortices in the mathematical 
model of the jet in cross flow. 
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Figure 6.- Coordinate system for relating jet trajectory 

to momentum flux. 











• NUMERICAL SOLUTION 
A ANALYTICAL SOLUTION 



T, 


Figure 10.- Spread rate parameters, based on published air jet data. 











Figure 12.- Rolling-up of two jet cross sections, one with circular, the other with 

elliptic exit nozzle. 





Three-Quarter View 


Three-Quarter View 


Figure 13.- Stacked cross sections Figure 14.- Stacked cross sections 
of a circular jet. of two circular jets in tandem 

arrangement. 
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A GENERAL JET EFFLUX SIMULATION MODEL 


11 


By Fred L. Heltsley and Richard A. Kroeger 
ARO, Inc. 

SUMMARY 


The complex flow fields induced by V/STOL aircraft in transition flight 
give rise to an interest in the problem of a jet issuing into a crossflowing 
stream. The greatest contributor to the induced flow is the roll up and ensuing 
vorticity of the viscous jet as it penetrates the free stream. This paper 
describes recent work done both to understand the forces of interaction of the 
jet and the field and to present a satisfactory mathematical modeling technique 
which will be useful in aircraft engineering. The results suggest a need for 
extensive experimental and analytical studies of the jet in a cross flow. 


INTRODUCTION 


Classical air vehicles produce an exhaust jet efflux which is oriented 
nearly parallel to the vehicle flight path. The literature contains many 
approaches to the solution of the flow resulting from the interaction of the 
jet and the surrounding field. For analytical simplicity, the flow may be con- 
sidered in three parts: 

1. The internal jet flow 

2. The flow external to the jet 

3. The boundary of interaction of these fields 

Regions one and three require viscous modeling; however, upon application of 
suitable boundary conditions region two may be effectively represented by 
singular solutions to Laplace's equation. 

V/STOL lift systems in transition flight offer an additional complication 
in analytical modeling. This is known as the "jet in a crosswind" problem. 

The three regions of interest are shown in figure 1. The purpose of this paper 
is to introduce a semiempirical method for providing solutions to region number 
two, which is of importance in the performance, stability and control of the 
air vehicle. Special emphasis is placed on generality and engineering 
simplicity. 


SYMBOLS 


A area 

b span of vortex element 

d jet diameter 
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du/dn rate of change of velocity with respect to a surface normal 
e entrainment per diameter 

F force acting on jet 

M momentum flux 

s distance along jet trajectory 

V velocity 

X,Y,Z jet coordinates in diameters 

1 unit vector 

p viscosity 

p density 

r shearing stress 

\| t angle between vortex element and impinging velocity 

F vortex filament strength 

Subscripts : 

e exit conditions 

i entry conditions 

j referring to jet 

n normal component 

p referring to potential solution 

t axial component 

oo referring to remote stream 

o initial condition 

1,2,3, 4, 5 order of shed vorticity 

PHYSICS OF THE JET 

The high degree of complexity of the flow created by a jet issuing normal 
to a free stream has led to a search for suitable variables for describing its 
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performance. For example, variables such as jet temperature, exhaust distor- 
tion, turbulence and total pressure distribution, Mach n umb er and Reynolds 
number may be used. On the other hand, these variables may be implied in an 
empirical set such as entrainment, wake separation, trailing sheet roll up, 
trajectory and cross-section geometry. Once a sufficient amount of empirical 
evidence has been compiled, these characteristics may further be condensed to 
the following set: 

1. Jet to free- stream velocity ratio 

2. Geometric distortion due to 

Vehicle interference 
Ground interference 
Multiple jet interaction. 

Some variables not expressed in this arrangement have disappeared as a result 
of compromise for simplicity or are implicit in the mathematical formulation. 

As a starting point, the simple free jet in a cross flow was studied. The 
trajectory (ref. l) was approximated by 

X = 1^) 2 (z)5 

where X,Z are the jet coordinates nondimensionalized by diameter from fig- 
ure 1 and V^/Vj is the ratio of the free stream to jet velocity ratio. The 
jet entrainment was found by integration of experimental data from reference 1 
to be 


( 1 ) 


e = 0.19 per diameter (2) 

which is the decimal rate of mass flow increase along a specified jet surface 
per exit diameter. Due to the lack of experimental data, the accuracy or vari- 
ation of this parameter with V^Vj was not known. 

Figure 2 shows the distortion of the total pressure profile along the path 
of the jet. Two major counterrotating vortices appear to be formed which 
increase in strength as the jet grows. These have been observed in many exper- 
imental studies and are of great importance to the flow. 

This has been a cursory description of the jet in light of its extreme com- 
plexity. It remains to be shown that such naive observations can be useful in 
defining the flow region external to the jet. 


VORTEX SINGULARITY MODELING 


A good approximation to the flow in the region where viscous effects are 
small may be made by the potential theory. This was applied by direct super- 
position of singular solution to Laplace’s equation with the appropriate boundary 
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conditions. Vortex singularities are very useful for three-dimensional lifting 
bodies. Other types of singularities including doublets, bicirculation, sources 
and sinks may be more appropriate for particular applications. Though some com- 
promise in numerical simplicity was paid, vortex singularities were used for 
the modeling for total system simplicity. 

Figure 3 shows the application of the vortex lattice method to a simple 
fan-in-wing configuration. Vortex singularities were used throughout with the 
following boundary conditions: 

1. V = uniform 

00 

2. The flow component normal to the wing surface vanishes 

3. The flow into the inlet = V^ 

J 

4. The flow normal to the jet efflux simulation tube surface must 

compensate for entrainment. 

Figure k shows an isometric view of the simulation of a propeller with its 
axis normal to the free stream. Some calculated streamlines are shown. The 
flow asymmetry at the prop inlet was induced by variation of the inlet flow 
between the individual boundary control points. 

Figure 5 shows how the jet roll-up was simulated by the trailing vortex 
elements. Two cross sections are shown for comparison of the model flexibility. 
Note that the downstream cross section (B-B) has grown in size and the number 
of trailing vortices used in its composition have increased. Figure 6 shows 
how these vortices were dispatched. The similarity with the low aspect ratio 
lifting surface analogy was noted. As the jet axis is traversed away from the 
exhaust nozzle, the number of trailing elements may be seen to continuously 
increase . 


INTERACTION BETWEEN THE JET AND THE SURROUNDING FLUID 


As the jet issues into the crossflowing stream, forces of action and 
reaction develop. The jet is deflected and distorted by the stream with a set 
of forces whose reaction causes an influence to be felt over the entire field. 
Figures 1 and 2 show the two primary results of these forces. It seems reason- 
able that if the surrounding fluid acts to form such distortions, and if that 
field may also be represented by singular solutions to Laplace's equation, then 
vortex elements may be employed to simulate the forcing elements in that field. 

In the model described here vortices were used for simulating wake block- 
age, entrainment and the trailing sheet. The latter was shown to be the greatest 
contributor to the flow distortion behind the jet. As in the case of wing theory, 
the strength of the trailing sheet was calculated by application of the correct 
boundary conditions. However, the presence of the viscous jet made the corre- 
sponding boundary difficult to define. 
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The objective of this section is to present an analytical model for the jet 
which will provide a solution to the field away from the jet. In effect, this 
is the inverse of the jet trajectory problem as the kinematics of the jet are 
a priori empirically specified. The vortex model described in the previous 
section was applied to the jet in a cross-wind problem. The trailing system 
was arranged as in figure 6. The results are shown in figure 7* Though the 
trailing system was very roughly approximated by the straight vortex lines, the 
experimental data shown in the plane of symmetry were matched quite nicely. 

The importance of the trailing system can be seen by comparison with figure 8 
wherein the same jet was simulated by entrainment and blockage only. 

The boundary condition sought in the modeling of the trailing vortex sys- 
tem was that no flow could cross the defined sheet. A one pass semiempirical 
computer solution was found which simultaneously solved all the vortex strengths 
as a consequence of the geometry and the boundary conditions. The strengths of 
the individual vortices in the trailing system were distributed along the jet 
axis as shown in figure 9« A cut was taken across the sheet, normal to the 
filaments, and the total wake vorticity on one side of the jet was obtained. 

This strength reflects the forces acting on the fluid field. A plot of the 
trailing element strength is shown in figure 10. 

It must be pointed out that the vortex strengths calculated in this exer- 
cise are effective strengths and would correspond to the viscous case only if 
the latter were 100 percent effective. It would be expected that if the actual 
vorticity could be experimentally measured it would appear higher than those 
which the mathematics requires in the simulation model. This argument must 
preclude the existence of a complex wake system with additional counterrotating 
filaments . 


As a means of further understanding the meaning of the results of the ana- 
lytical modeling, a more in-depth look was taken at the development of the 
trailing system. This approach required a detailed experimental knowledge of 
the internal jet flow. Since such data were nonexistent, a particular jet 
boundary was defined and information was extrapolated from the experiments of 
Shandorov. The goal was to compute, for comparison, the circulation strength 
from the calculated induced forces of reaction between the jet and the sur- 
rounding fluid. 

Figure 11 shows the jet induced force F_ which was to be found. The com- 
ponent Fn is normal to the jet axis while is tangential. The momentum 

flux rate along the jet is 


M = 



pV 2 dAl 


(3) 


where p is the fluid density, A is the jet cross-sectional area, V is the 
local jet velocity component and T is a unit vector in the direction of local 
flow. If the momentum flux was known along the jet, the interaction force 
could be easily calculated. In finite element form, the force acting on a jet 
control volume element is given by 
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-F 


= J pV 2 e d^lg - J pV i 2 dA i l i (4) 

The subscripts i and e denote entry and exit, respectively. The control 
volume was defined as an axial element of the jet whose extremes are planes 
normal to the jet axis. The peripheral surface coincided with the locus of the 
velocity profile where the shear stresses become sufficiently small to divide 
the jet from the field which was called potential. This roughly matched the 
outer boundary of the profiles shown in figure 2. The need becomes apparent 
for the definition of a maximum allowable shearing stress 


T 


P 



(5) 


which would not invalidate a potential solution. The derivative in equation (5) 
is the rate of change of velocity with respect to a surface normal. 

The next step was to calculate the strength of the bound vortices which 
would produce the force described in equation (4). From elementary three- 
dimensional lifting theory the force on a bound vortex is given by 


F = pV X fb 


( 6 ) 


where T is the strength and b is the span of the lifting element. Next, 
the force on the jet was assumed equal to that induced by a group of distributed 
bound vortices, each with a pair of trailing elements. This arrangement is 
shown in figure 9 • Their strengths were found by equating equations (4) and (6). 
Thus , 


r = 


F 

pVb sin \|r 


1 


( 7 ) 


Since the jet velocity vector imposed on T was normal to the bound element, 
sin \|r was taken as unity and the scalar magnitude of vorticity written 


r 


F 

pVb 


( 8 ) 


The rationale for using the jet velocity in this equation was that the bound 
element was assumed to span the jet. In the potential model, the average inter- 
nal jet flow velocity could be matched with that in the real jet. As each vor- 
tex strength was determined by the induced velocity over the bound filament, 
the jet velocity fulfilled the requirement. 
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Upon breaking the jet into longitudinal control volumes, equation (8) was 
calculated and integrated along the trajectory s for the total sheet strength. 

The propeller model introduced earlier was used as an example of the appli- 
cation of this technique. This example was chosen at the risk of making the 
foregoing developments appear somewhat precocious in light of the experimental 
data at hand. 

First, the jet characteristics were calculated using the data from refer- 
ence 2. Some guess work was necessary to arrive at the integrals of the flow. 
Figure 12 shows the axial momentum distribution. It is int eresting to note 
that the momentum starts to decrease quite rapidly upon exit to the free stream 
and then begins to increase. Though the validity of the data application was 
somewhat in question, it is felt that the resulting flux could be reasonable. 

The mechanism for such distortions was felt to be the trailing sheet rapidly 
rolling up at the knee of the trajectory causing some of the axial momentum to 
be changed to rotational momentum. Later, the rotational component is appar- 
ently reintroduced in the axial direction as the vortices are bent. 

3y applying equation (4) to this data, figure 13 was obtained. The force 
vectors shown are the negative of those forces which act on the fluid field. 
Figure Ik shows the strength of the trailing sheet calculated from the forces 
on the fluid field. The width of the jet was taken from figure 2 to be the span 
of the bound vortex filament . 

The forces of interaction between the jet and the free stre am have been 
defined and a method of semiempirically evaluating the resulting remote poten- 
tial flow field described. The interesting comparison between figures 10 
and 14 shows the existence of a method for imposing vorticity in the trailing 
jet wake either by direct potential calculations or by derivatives of experi- 
mental data. 


CONCLUSIONS AND RECOMMENDATIONS 


A semiempirical solution of the flow field induced by a jet in a crossflow 
has been presented. It is yet somewhat premature to speculate on the quantita- 
tive effectiveness of the technique even though figure 7 is very well ma tched 
in the plane of symmetry. 

The comparative trends of figures 10 and l4 show the logic of the simula- 
tion to be in order. The shortcoming in the proof is in the limitation of 
applications due to the availability of adequate experimental information. 
Experiments should be conducted which define the jet boundary and momentum flux 
so that the induced forces may be computed accurately. Flow-field information 
away from the jet should be compiled to complete the verification of the semi- 
empirical model. Finally, computer experimentation should be conducted to 
define the boundary conditions which provide adequate analytical representation. 
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Figure 1.- Areas of interest. 



Figure 2.- Total pressure per cross section. 








Experimental Data 



Figure 8.- Comparison of analytical and experimental results with only blockage 

and entrainment encountered for. 
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Figure 9*- Trailing vortex nomenclature. 



Figure 10.- Analytical distribution of shed vorticity 





CALCULATION OF JET INTERFERENCE EFFECTS ON V/STOL 
AIRCRAFT BY A NONPLANAR POTENTIAL FLOW METHOD 
By Paul E. Rubbert 
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SUMMARY 


A method of calculating the flow field about V/STOL aircraft with a 
potential flow model is reviewed with emphasis on the flow modeling require- 
ments. Comparisons with experimental data are used to demonstrate the impor- 
tance of proper jet simulation. The dominant flow features associated with a 
jet are reviewed and a theoretical model is proposed which recognizes these 
features. An experimental program is suggested as a means for determining 
various parameters appearing in the jet model. 


INTRODUCTION 


The calculation of jet interference effects on V/STOL aircraft is con- 
siderably more complex than the classical treatment of interference in the 
context of linearized lifting line or lifting surface theory based on small 
disturbance concepts. Within the linearized small disturbance framework, 
interference effects are interpreted in terms of superimposed downwash fields 
which are akin to camber and angle-of-attack changes, and buoyancy effects 
associated with the longitudinal interference velocity field. Such concepts 
are based on the requirement that all disturbance velocities are small com- 
pared with the freestream speed, and that the rate of variation of the inter- 
ference velocity components within the confines of the airplane dimension re- 
main small. Such is not the case for V/STOL aircraft, which typically operate 
with high inlet inflow velocities, are subject to interference effects from 
jet entrainment that vary widely over distances comparable to the aircraft 
dimensions, and can operate at low forward speeds. 

Thus, one is led to abandon the simplified small disturbance concepts 
and adopt a less approximate scheme. The present paper describes such a 
scheme (refs. 1 and 2) and focuses attention on its limitations and on the 
particular problems associated with the efflux jet. 
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SYMBOLS 


an 



c 


m 

n 

P 

q 

r 

s 

S 

Si 

s 2 

Voo 

VlNFLCW/Voo 

x,y,z 

a 


n 

e 

p 

<p 


entrainment coefficient (n = 0, 1, . . .) 

drag coefficient 

lift coefficient 

pitching moment coefficient 

pressure coefficient 

chord 

source strength 

distance normal to a surface 

unit vector normal to centerline in plane of trajectory curve 

velocity 

distance 

distance along the jet centerline 

boundary in a potential flow model 

portion of S consisting of the body surface 

portion of S encompassing the jet 

freestream speed 

jet velocity ratio 

Cartesian coordinates 

angle of attack 

span fraction 

angle around the jet 

doublet strength 

singularity strength on panel with index j 
velocity potential 
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POTENTIAL FLOW CALCULATION METHOD 


Modeling the Real Flow 


The present scheme is based upon a capability for solving Laplace's 
equation (irrotational incompressible flow) numerically in three dimensions 
subject to arbitrary but properly posed boundary conditions; namely, Neumann, 
Dinchlet or mixed boundary conditions, or variants thereof, on closed 
boundaries whose special orientation is known. The initial task, when apply- 
ing this method to a V/STOL problem, is to formulate a theoretical model of 
the real flow that meets the above requirements while yet retaining the 
essential features of the flow. 


Figure 1 shows such a theoretical model that was developed (refs. 1 and 
2) for a fan— in— wing aircraft. The assumptions inherent in its formulation 
are as follows: 


1. The outer flow is assumed to be incompressible and irrotational. 

2. Boundary layer displacement effects are ignored. 

3. Regions of separated flow are not simulated. 

4. The special orientation of the trailing vortex sheet is prescribed 
on the basis of a priori knowledge or experience. 

5. The inflow distribution at the fan face must be prescribed. 

6. The trajectory and distribution of entrainment on surfaces en- 
compassing the efflux jets must be prescribed. For this model 
the entrainment was assumed to be zero. 

7« The fuselage has been deleted for computational economy. 

These assumptions serve to reduce the problem to manageable form. They 
also impose certain limitations on the usefulness and domain of applicability 
of the present calculation procedure as discussed in refs. 1 and 2. 


Numerical Procedure 

The numerical solution procedure is based on the concept of source and 
doublet (or vortex) distributions on the boundary surfaces. Distributions of 
these singularities are approximated numerically by networks of small quadri- 
lateral panels with the singularity strength of each panel set eq ual to an 
arbitrary constant value, say cfj . Figure 2 shows the paneling used for the 
model of Fig. 1. Certain portions have been omitted for clarity. This 
numerical panel approximation reduces the problem to the determination of a 
finite number of discrete singularity strengths, one for each panel. 
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One point in each panel, the "boundary point", is selected as the point 
where the boundary conditions are to be satisfied. The requirement that the 
boundary conditions be satisfied at all boundary points produces a system of 
linear algebraic equations whose solution yields the desired values of CJj . 
From these the entire velocity and pressure field are calculated. References 
1 and 2 contain a more detailed description of the procedure. 

Experience has shown that one of the key points in using the method is 
the selection of the singularity type (either source or doublet) to represent 
the various surfaces. There are an infinite number of possible choices, and 
the user must be careful to select a combination of source and doublet panels 
that guarantees a unique solution. The basic concepts contained in Chapter 
III of Lamb (ref. 3) are invaluable for providing the proper insight. Another 
guideline for the choice of singularity distributions is that the panel ap- 
proximation is limited accuracy- wise in a manner related to the change in 
singularity strength from panel to panel (refs. 1 and 2). This limitation 
imposes restrictions on both panel size and spacing and on the initial selec- 
tion of singularity types. Experience with the numerical method together with 
a firm understanding of the fundamentals of potential flow theory serve to 
minimize the importance of these restrictions. These words of caution are 
offered to prevent the uninitiated from becoming overly enamored at the out- 
set by the apparent simplicity and power of the method. Successful applica- 
tion of this method to new and different problems is strongly dependent on 
the user's ability to model the real flow properly and to select the appro- 
priate singularity and panel arrangements. 


Initial Results and Problem Areas Uncovered 

The present method, when properly used, provides a finely detailed solu- 
tion of the inviscid model problem. The solution provides surface and flow 
field pressures and velocities, and forces and moments obtained from inte- 
grated surface pressures. Figure 3 » which has been reproduced from ref. 2, 
shows a comparison between the measured and calculated forces on the fan-in- 
wing model described in Fig. 1. These were obtained with a moderate inlet 
velocity ratio of ViNFLOW/^oo = 1.667 and with the efflux jet vectored aft by 
20°. No entrainment was assumed on the periphery of the jet boundary. The 
various theoretical curves denote different assumed inflow distributions at 
the fan face and demonstrate the sensitivity of the calculated forces to this 
parameter. The differences between the curves arise primarily from the com- 
puted internal force on the fan assembly calculated by a momentum balance 
across the fan. This force is added to the integrated pressure force on the 
external surfaces to give the total force acting on the configuration. The 
calculated drag levels are lower than the measured values, which is to be 
expected since the profile drag is not included in the calculated values. 

The variation between the pitching moment curves can be attributed to the 
absence of the fuselage on the theoretical model and to boundary layer dis- 
placement and local separation effects. 
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A comparison of the calculated surface pressures with experiment often 
shows a variation between them which increases with inlet velocity ratio. 
Figure 4 , taken from ref.. 2 , shows an example with an inlet velocity ratio of 
3.636. Large differences are apparent, particularly in the vicinity of the 
fan and flap regions. Gross deviations of upper surface pressures near the 
f an are probably caused by inaccuracies in the assumed inflow distribution. 
Large deviations on the lower surface near the fan exit can be attributed to 
jet entrainment effects which were not simulated theoretically. The discrep- 
ancies in the flap pressures, which are greatest at the station directly aft 
of the jet, may also be due in part to the strong influence of the entraining 
jet as well as to viscous effects which are important in this region. This 
example serves to illuminate the importance of both inflow distribution and 
proper jet representation in the theoretical treatment of V/STOL problems. 
Lesser discrepancies that exist in regions more remote from the fan must also 
be due partly to these effects, because experience with the method applied to 
conventional airplanes has demonstrated an order of magnitude improvement in 
detailed surface pressure comparisons when lift fans are not present. 

Thus it appears of utmost importance for V/STOL applications to improve 
our capability to treat these effects properly. The question of what is the 
proper inflow distribution can be answered partly by further experimental 
measurements and correlation with the basic parameters, and involves no change 
in the basic theoretical or numerical representation. The representation of 
the jet efflux, however, poses a more difficult problem. 


REPRESENTATION OF A JET IN A POTENTIAL FLOW 

Let us examine the jet representation in some detail and attempt to 
identify those parameters which can have the most significant effects on the 
calculated flow near the surfaces of the airplane. 


Basic Concept 

The representation of a viscous, turbulent jet by means of a potential 
flow model is based on the following concept. Consider the example of a 
jet emerging from a body into a flow field. The flow external to the jet is 
considered to be inviscid and irrotational. This allows the definition of a 
velocity potential (p associated with the inviscid external flow. Let us 
define a closed boundary S consisting of two parts; one part being the 
body surface and the other part S2 being any arbitrary surface encompassing 
the turbulent jet flow, Fig. 5. Green's theorem applied to the exterior 
potential flow field (ref. 3) leads to the following expression for the 
potential (p . 


185 



( 1 ) 


where 


* <*>».*> ’ - ds + -Ir (-B d * 


SOURCE 

DISTRIBUTION 


DOUBLET DISTRIBUTION 


m(s) = source strength on S 
|i(s) = doublet strength on S 

r = distance from S to a field point (x,y,z) 
n denotes the direction normal to S 


The terms on the right-hand side of Eqn. (l) can be interpreted as dis- 
tributions of source and doublet singularities on the boundary surfaces. The 
singularity strengths, m(s) and ja(s), can be obtained as the solution of 
the boundary value problem obtained by letting the field point (x,y,z) in 
Eqn. (l) approach the boundary surface and differentiating the entire equation 
with respect to the surface normal, giving 


9£fS)_ 

an 



rf v + 




( 2 ) 


The numerical method previously described furnishes solutions for m(s) and 
H(s) as a function of the normal velocity component on the boundary 

surfaces. Thus the entire flow field described by <p Cx,y,z) is ultimately 
dependent upon a knowledge of the normal velocity component ^-(S) on the 
boundary sur faces. n 

In the usual problem, is known on the body surface S]_, which in- 

cludes the external surface of the airplane through which no flow can pass, 
plus the portion of across an inlet on which the value of , which 

provides the inflow distribution, must be specified. If 4j~ were known on 
S 2 as well, then the influence of the jet on the external flow would be 
accurately established. However, this is not the case, and so our ability to 
simulate the effect of a jet on the external flow ultimately reduces to a 
need for estimating the distribution of the normal velocity component -^-(^ 2 ) 
on any arbitrary boundary surrounding the turbulent jet. 

This needed information must be obtained from experimental observation. 
The recent literature (ref. 4) indicates that experimentalists are acquiring 
greater insight into this problem. Ultimately, it is hoped that their work 
will lead to a set of empirical relationships which provide the information 
necessary for an accurate jet representation with a potential flow model. 
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Important Features of a Jet 


These experimental efforts lead one to try and identify the most 
important characteristics of the theoretical jet model in order to reduce the 
task from one of predicting the detailed 4^!L distribution on an arbitrary 
boundary surrounding the jet to that of establishing those gross features of 
the jet which produce the greatest influence on the flow near the configura- 
tion surfaces. Let us consider first the jet trajectory. There is no need to 
place S 2 exactly at the jet boundary, but it is necessary to estimate the 
approximate path of the jet. Its initial path is usually known from the 
nozzle orientation and so the problem reduces to the establishment of the 
trajectory some distance downstream from the exit where the jet is turning in 
response to its interaction with the cross flow. The required accuracy of the 
trajectory estimation is strongly dependent on the proximity of the jet to any 
of the configuration components (wing, tail, etc.). If the jet passes close 
to any parts of the configuration, its path must be known with greater accu- 
racy than if it is directed away from the configuration, because the error in 
the velocity at the configuration surface caused by an error hr in the 


location of a segnent of the jet located a distance r from the configuration 
surface behaves as . 


Experimental observations of a single jet emerging into a uniform cross 
flow have led to the formulation of a well-known empirical trajectory equa- 
tion by R. J. Margason (ref. 4). This equation is adequate for those con- 
figurations where the jet does not pass near any of the configuration com- 
ponents. In cases where the jet does pass in proximity to the configuration, 
it may be necessary to take into consideration the effect of the local air- 
plane-induced flow field on the jet, and also perhaps the detailed properties 
of the jet flow which govern to some extent the rate of entrainment and the 
subsequent turning force applied to the jet by the outer flow. These effects 
greatly complicate the problem, and it thus appears that in the near future 
the application of potential flow methods will be restricted to configurations 
having jets directed away from any configuration components. 


Having determined the trajectory of the jet centerline, the remaining 
problems involve the determination of an appropriate approximation for the 
location of the surface S 2 with respect to the jet centerline, and the 
specification of the velocity component on S 2 . Consider first the 

limited portion of the jet immediately adjacent to the nozzle exit and ex- 
tending down the jet perhaps one or two nozzle diameters, as shown in Fig. 6. 
The jet boundary in this limited region must retain quite closely the known 
cross-sectional shape of the jet nozzle and expands in a manner that could be 
approximated by known relationships governing the rate of spreading of axi- 
symmetric turbulent jets. Such an approximation could easily be checked and 
possibly improved by means of flow visualization experiments. Actually, it 
is probably not necessary to attempt that degree of sophistication in order to 
locate S 2 « A simple approximation consisting of a nonexpanding jet cross 
section retaining the shape of the jet exit. Fig. 6, may be adequate as dis- 
cussed below. Referring to Fig. 6, the normal velocity components on the 
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actual jet boundary can be analytically continued by Taylor series expansion 
to the proposed surface S 2 by an equation of the form 


9 (r»Ar) _ _25£_<r) . a- d 2( P(r) 
9n an f Ar Fn'ar 


+ 0(Ar 2 ) 


(3) 


where 


<W (r+Ar) 

an 


= normal velocity component on the actual jet boundary 


_ 2 £(r) 

an 


= normal velocity component at a corresponding point on 
the surface S 2 


r = radial distance from the jet boundary to the surface 


If Ar « r , which is the case in this limited region near the jet exit, 

then |£(rfAr) may be replaced by |£( r ) . Or, one could seek to refine the 

value of specified on S 2 by including the second term on the right-hand 

side of Eqn. (3). Actually, as will become apparent, it is probably as easy to 
establish the value of on S 2 needed to produce the correct jet influence 

as to establish the correct normal velocity on the jet boundary. If that is 
true, there is no need for a jet boundary representation in this region that 
is more sophisticated than the simple one proposed above. 

Let us continue our examination of this limited portion of the jet and 
turn to the property most important and most difficult to determine; namely, 
the value of on S 2 . Figure 7» taken from ref. 4, shows the pressure 

distribution measured on the surface of a plate which has a jet exhausting 
normal to the plate in the presence of a freestream. The pressure pattern 
displays fore and aft asymmetry, which strongly indicates that the flow 
around the jet is not similar to attached flow about a cylinder. Reference 4 
also shows oil flow patterns on the plate which display the features sketched 
in Fig. 8 . The gross features indicate that the flow is deflected about the 
jet as with a cylinder. However, the details show flow into the jet along 
the sides, a stagnation point on the plate aft of the jet with upstream flow 
ahead of the stagnation point into the jet, and a well-defined oil accumula- 
tion line streaming aft of the jet. There is obviously entrainment of 
exterior flow into the jet, both on the sides and into the aft portion. Thus 
we are faced with the prospect of having to determine the strength and dis- 
tribution of this inflow around the jet. In a potential flow model this in- 
flow, which appears as the value of to be specified on the boundary S 2 , 

is the quantity which most strongly controls the detailed features of the 
flow near the jet exit. 


Let us now consider the downstream portions of the jet. Flow visualiza- 
tion studies have revealed that in the presence of a cross flow, the jet pro- 
gressively rolls up into a vortex pair, such as sketched in Fig. 9» Photo- 
graphs of the cross section of a jet wake such as given in Fig. 21 of ref. 4 
reveal a strong induced inflow into the aft portion of the jet (see Fig. 9 ). 
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Thus it appears that the downstream portions of the jet can be characterized 
as a vortex pair that is continuously entraining fluid from the outer in- 
viscid flow into the viscous jet. 

The main objective in attempting to model this jet flow is to reproduce 
the proper induced effects in the vicinity of the configuration surfaces. 

The detailed features of the flow in the immediate vicinity of the jet wake 
itself do not matter as long as its long-range effects are correct. This is 
true in so far as the distance of the jet from the configuration components 
is large compared with the mean jet radius. Hence let us attempt to identify 
those features of the jet wake which dominate its long-range influence. 

One can speak of the influence of a jet in terms of the velocities in- 
duced by a distribution of singularities on a surface S 2 surrounding the 
viscous jet. This representation is exact provided the singularity strengths 
are known. Alternately, the surface distribution of singularities can be re- 
placed by a multipole expansion on the jet centerline, which will be of the 
form 


<?(x,y,z) * - -±-f ds ♦ ofo) W 


where 


r = distance from jet centerline to the point (x,y,z) 

P = unit vector normal to the centerline in the plane of 
the trajectory curve. Fig. 10 

The first term on the right hand side contains the long range effect of 
the entrainment, or inflow into the jet. It represents a distribution of 
sinks along the jet centerline, with their net strength per length As approx- 
imately equal to the integral of the inflow velocity — 1-£. around a strip of 

width As encompassing the viscous jet, as sketched in Fig. 10. At very 
large distances from the jet, this term alone will produce the primary effect, 
because all other terms die off more rapidly with increasing r . 

At intermediate distances, the second term may also be important. It 
represents a distribution of doublets with axes normal to the jet centerline 
in the plane of the paper, Fig. 10, and produces the major long range effect 
due to the distribution of entrainment or inflow around the jet circumfer- 
ence. Furthermore, its effect at large distances will be the same as that due 
to a vortex pair oriented in the manner of Fig. 9 . This latter feature can be 
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derived by expanding the potential of a two-dimensional vortex pair spaced a 
distance 2d apart in powers of d/r . The result is 



<p(x,z) s 
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27 T 
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2 7T 


2cosfl 
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( 5 ) 


The first term on the right represents a two-dimensional doublet located 
at the origin. The difference between the influence of the doublet and the 
vortex pair at distances larger than the jet radius ( ~ d) is 0(d^/ r^). 

Thus one is able to associate with a sink and doublet distribution along 
the jet axis the major features of the jet that have yet been identified by 
flow visualization; namely, a net entrainment into the jet with the distribu- 
tion of entrainment around the circumference dominated by the vortex pair flow 
pattern. Figure 11 shows a two-dimensional velocity pattern about a super- 
imposed sink and doublet which clearly displays these characteristics. It 
should be noted that an actual distribution of sinks and doublets on the axis 
is not necessary to produce these effects. The jet can be represented by any 
convenient surface cross section having boundary conditions producing the sink 
and doublet effects. As an example, a circular cross section of radius R with 
boundary conditions of the form 

(R.0) _ Jj , Cg cosfl (6) 

an “ 27 tr + 2ttR2 


will produce the same effect as a sink and doublet at the circle center having 
strength and C 2 respectively. 

For configurations where the jet does not pass in close proximity to any 
aircraft component it is unlikely that the finer details contained in the 
higher order multipoles will have a significant effect because of their rapid 
attenuation with distance. In any event, the estimation of the jet trajectory 
itself with sufficient accuracy becomes very difficult with close proximity 
jets, as previously noted, so it appears that, for the present, the use of 
potential flow methods shall be restricted to configurations with remote jets 
for which the sink-doublet representation is applicable. 
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In summary, the most important properties of the jet representation are 
the following: 

1. Trajectory. At present there exists an empirical trajectory 
equation which is adequate for jets not in close proximity to 
the configuration. 

2. Near the jet exit, a surface S 2 is required. It can have the 
cross-sectional shape of the jet exit and need not be tapered 
to correspond to the expansion of the jet. 

3. A detailed distribution of 4^- on S 2 near the jet exit is re- 
quired. This probably has the greatest influence of any of the 
jet properties. Our present knowledge concerning the distribution 
°f -|jj- on S 2 is only qualitative. 

4. The influence of the downstream portions of the jet must contain 
sink and doublet effects. The magnitude of these effects is pre- 
sently unknown. 


PROPOSED POTENTIAL FLOW MODEL OF A JET 

In view of these properties, I propose the model shown in Fig. 12 as 
suitable and convenient for use with the potential flow calculation procedure 
described herein. The cross-sectional shape of the jet exit is retained all 
along the jet as a matter of convenience. The trajectory is established by 
Margason's equation (ref. 4). The remaining quantity to be established is the 
distribution of on the surface, and this must come from experimental 

observation. Ideally, one would like to express the ^ distribution in 
terms of a few parameters and attempt to establish a correlation between these 
parameters and the jet deflection angle and velocity ratio. One possible 
parametric representation is the expression 


s y; q n (s) cosnfl 


(7) 


The first terra, a Q (s) , provides the longitudinal distribution of the 
total inflow or entrainment along the jet. The second coefficient, a^(s) gives 
the longitudinal variation of the doublet or vortex pair effect. Higher order 
terms, if needed, will be important only near the jet exit. Their influence 
on the configuration rapidly disappears with increasing s, and their presence 
or absence for large s will not affect the flow in the vicinity of the con- 
figuration surfaces. This latter behavior suggests that the variation of 
with s for n > 2 can possibly be neglected, leaving only constants correspond- 
ing to the value of a n (0) . This possibility should certainly be explored. 
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SUGGESTED EXPERIMENTAL PROGRAM FOR THE 
DETERMINATION OF THE JET PARAMETERS 

With this type of representation, one is faced with the necessity of 
determining, as a minimum, the functions a Q (s) and a-^(s) plus perhaps some 
constants for the next few a^s- These will have to be established from 
experimental observation, but how? 

One possible approach would be to examine the results of an experiment 
conducted with a jet emerging into a cross flow at various angles from a 
large flat plate as sketched in Fig. 13a. The plate must be large enough to 
avoid boundary layer separation at the edges due to jet induction effects 
and should contain a detailed static pressure survey. The basic idea is to 
adjust the a n 's in a theoretical solution of this problem, Fig. 13b, to re- 
produce the measured pressure distributions on the plate. The final result 
would be a number of curves or fitted equations describing the variation of 
the a n 's with jet velocity ratio, thrust vector angle, and a jet dynamic pres- 
sure decay parameter (see ref. 4). An attractive feature of the scheme is 
that theoretical solutions for severed different distributions of a n (s) can be 
obtedned nearly as economically as for one, since the influence coefficient 
matrix in the numericed solution procedure is independent of the a n 's. Or, 
one could attempt to mechanize the process by expressing Sq(s) and a^(s) 
as power series in s and calculating, with a potential flow model, the 
separate influence of each term in the series. This could also be done 
economically due to the constancy of the influence coefficient matrix. These 
individual solutions, or flows, would serve as basic building blocks which 
could be linearly superimposed to produce any desired flow pattern. Armed 
with these, an investigator could readily construct all possible theoretical 
flows. 


A point to note in attempting to match theoretical and experimental plate 
pressures is that the pressures away from the jet exit are governed mainly by 
Sq(s) and a^(s), whereas the pressures adjacent to the exit are strongly in- 
fluenced by the higher order a n *s . 


CONCLUDING REMARKS 

The results of previous potential flow calculations for V/STOL applica- 
tions have demonstrated a need for improved modeling of the jet flow. With 
this objective in mind, a theoretical model for a jet has been proposed which 
is based on the observed major features of real jet flows. An experimental 
program is suggested as a means of obtaining a correlation relating the jet 
boundary conditions in the theoretical model to the jet parameters. 
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• Arbitrary planform, thickness, camber, inlet geometry, and fan inflow distribution 

• Source panels on wing and centerbody surfaces 

• Vortices on efflux tube and trailing sheet 

• Internal and fan face vortices omitted for clarity 


Figure 2.- Panel arrangement on fan-in-wing configuration. 
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PRESSURE COEFFICIENT, C 
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Figure 4.- Surface pressure comparison on fan-in-wing configuration. 
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Figure 4 .- Concluded. 


199 


ARBITRARY BODY. 



Figure 5 • - Potential flow boundaries. 



Figure 6.- Initial portion of jet. 
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Figure 11. - Velocity pattern about a superimposed sink and doublet in a 

uniform stream. 



Figure 12.- Proposed model of jet. 





INVISCID MODELS FOR THE PRESSURE INDUCED BY A 
JET TRANSVERSE TO A SUBSONIC STREAM* 

By Robert Rosen, Norbert A. Durando, and Louis A. Cassel 
McDonnell Douglas Astronautics Company- -Western Division 


SUMMARY 


Methods are presented for approximately predicting the pressure 
distribution on a flat plate in the vicinity of a transverse jet exhausting into 
a subsonic uniform stream. A method is formulated for empirically fitting 
the available data to enable interpolation and extrapolation. In addition, sev- 
eral inviscid flow models are postulated, based on incompressible potential 
flow theory. Some models were found reasonably successful while others 
were not. This paper presents the formulation of the empirical fit, analysis 
Oi die m«_,<_iels, and comparison of results from the models with experimental 
data. 


Th® first subject discussed is the formulation of an empirical technique 
for fitting the data. Following this, results of two approaches to potential 
modeling of the flow are presented. In the first approach, all of the known 
significant phenomena, such as entrainment and blockage, are included in an 
approximate way. Singularities are distributed to account for these effects 
but their strengths are left arbitrary. A pressure-matching scheme is 
employed to quantitatively determine the strengths. The potential flow 
methods used in the second approach are more direct for the purpose of the 
study. They consist of postulating a flow field which will give the proper 
qualitative pressure distribution and adjusting constants to fit this distribution 
as closely as possible. 


INTRODUCTION 


Interaction of a highly under expanded sonic or supersonic jet with a sub- 
sonic transverse flow represents a jet interaction (JI) problem encountered in 
the design of many reaction-controlled atmospheric flight vehicles. Conse- 
quently, an effort is being made to develop analytical methods for predicting 
the pressure distribution, and therefore the forces and moments produced by 
the jet, and the interaction flow field for various flow conditions and vehicle 
geometries. Some results of such an effort, being conducted by the McDonnell 
Douglas Astronautics Company-- Western Division under the sponsorship of the 
U. S. Army Missile Command, are described in this paper. The work dis- 
cussed here is limited to an analysis of the pressure on a flat plate in the 
vicinity of a jet exhausting transverse to the plate. Both empirical methods 
and incompressible potential flow modeling will be described. 


*This research supported by the Advanced Systems Laboratory, U. S. Army 
Missile Command under Contract number DAAH01-68-C- 19- 19. 


Detailed descriptions of the flow surrounding a subsonic jet exhausting 
transverse to a subsonic uniform stream, based on the results of careful 
experiments, are available in references 1 and 2. Similar detailed experi- 
ments for sonic or supersonic underexpanded jets do not exist. In the opinion 
of the authors, however, adding a sonic or supersonic jet only complicates 
the flow field further. In both cases, the complexity of the flow field is such 
that an exact analytical description of it is beyond current developments in 
fluid mechanics. Consequently, empirical methods and potential flow anal- 
ogies to the fluid mechanics, as known from available experimental data, 
have been employed in the initial efforts to predict pressure distributions. 

The objectives of effort described in this paper are to devise techniques for 
approximately predicting the forces and moments on bodies from which the 
jet exhausts and not necessarily to derive detailed analytical descriptions of 
the flow. The direct problem of interest initially hasj, therefore., been taken to 
be the analytical modeling of the pressure distribution on a flat plate in the 
plane of the jet exit. 

In the analysis that follows, it will be assumed that the fluid behaves as 
if it were incompressible and ir rotational. Viscosity, if taken into account 
at all, will be introduced through its bulk effects, such as entrainment, and 
will not enter any of the equations. Thus^a potential exists for this flow and 
solution will be sought based on incompressible potential flow theory. Some 
other authors (refs. 3 through 6) have also used inviscid flow theory to 
describe the problem but with limited success. 

Although some flat-plate pressure data do exist for an underexpanded jet 
exhausting into a subsonic stream (refs. 7 and 8), they are rather incomplete. 
Recent experiments have been carried out by the U, S. Army Missile Command 
which are quite comprehensive, but data from these tests are not yet available. 
For lack of data more appropriate to the present study, the subsonic jet data 
of Vogler (ref. 9) will be used. It is felt that theoretical and empirical 
models which can be made to agree with these data can also be made to agree 
with sonic and supersonic jet data. 

The coordinate system used in the analysis of the pressure distribution 
on a flat plate is shown in figure 1. The polar coordinate, 0, is measured 
counterclockwise from the mainstream velocity. The radial coordinate, 
indicated as r* when dimensional, is measured from the jet center. When non- 
dimensionalized by the jet diameter, the radial coordinate is indicated as r . 


SYMBOLS 



singularity strengths^ 

C 

P 

pressure coefficient 

c 

Fourier coefficients, 

n 


d. 

J 

jet nozzle diameter 


where 


where 


n = 0, 1, . . . 


n = 0, 1 ? • . . 
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F 

k 

M 

P 

<1 

r 

r 

U 


* * * 
x ,y ,z 


z 

9 


P 

Subscripts : 

i 

o 

CO 


force 

scale factor 
moment 
pressure 
dynamic pressure 

radial distance normalized by jet diameter 

dimensional radial coordinate 

velocity 

complex potential 

dimensional Cartesian coordinate 

complex variable normalized by dj 

azimuthal angle measured positive in the counterclockwise 
direction from the leeward ray 

density 

jet 

vortex position 
free stream 


A bar over a symbol denotes a complex conjugate. 


DISCUSSION 
Empirical Data Fit 

A typical plot of pressure coefficient, Cp, versus 8 at a given value 
of r is shown in figure 2 based on the data of Vogler. Because of the smooth- 
ness of these data for 0 > 30 ° and the fact that the data are even, it appears 
that the data can be approximated quite well at a given radius by the first few 
terms in a Fourier cosine series. 


m 



n=0 


cos 


(n9) 


( 1 ) 


where the coefficients c n are given by their usual integral definitions. Fig- 
ures 3 and 14- show comparisons of two- and a three-term series to the data at 
two values of r. It can be seen that the truncated Fourier series does indeed 
provide a good representation of the 0 dependence of the data. The same 
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results have also been found to hold true at other values of r and other velocity 
ratios (U /U.). 

To complete the empirical fit, the radial dependence of the coefficients c n 
must be found for the various values of the velocity ratio. This dependence is 
shown empirically in figures 5 and 6. It can be seen that the variations of cq 
and Cj as functions of r and the velocity ratio are well behaved. It would be 
a simple matter to obtain functional relations for this dependence to com- 
pletely determine the empirical fit. However, a better way has been found to 
use this information through potential flow models. 


Phenomenological Models 

Two models will be considered under this heading. The first, called the 
doublet model, consists of the superposition of a free stream, a source (or 
sink), and a doublet. The second flow model, called the doublet- vortex model, 
is the same as the first except'that two counterrotating vortices are added 
downstream of the origin. A sketch of this model is shown in figure 7. Both 
of these models are two dimensional, in that it is assumed that there are no 
variations in flow properties in the direction normal to the plate. 

The complex potentials of these simple models are well known. The only 
quantities to determine are the strengths of the singularities and locations of 
the vortices for the doublet- vortex model. In the present analysis, the com- 
plex potential for the flow field is written in terms of the unknown singularity 
strengths. From this complex potential, an expression for the pressure 
coefficient in terms of cosine functions is obtained and put in the form of 
equation (1). By setting each coefficient in the derived pres sure- coefficient 
expression equal to the corresponding one in the Fourier series obtained from 
the data at a given value of r, the strengths of the source, doublet, etc. , are 
found. This is done for several values of the velocity ratio, and curves of the 
source strength, etc.^ are obtained as a function of the velocity ratio. The 
benefits of this technique are twofold: First, with its use a limited amount of 

data can be greatly extended; second, and perhaps more important, by seeing 
which singularities are strongest, the dominant factors in the flow field can be 
inferred. This information could be important in trying to relate the results 
obtained for a flat plate to an axisymmetric body. 

The model will yield results that are accurate at one value of r, the one 
where the data and model were matched. The accuracy of results at other 
values will depend upon how closely the r dependence of the model depicts 
the actual situation. There is an assumption implicit in the above analysis 
which must be satisfied to obtain reasonable results. The quantity which will 
ultimately have to agree with experiments is the expression for the complete 
pressure coefficient, not just the first few terms in its Fourier series repre- 
sentation. In order for theoretical and experimental pressure distributions to 
agree, the first few terms in the series derived from the complex potential 
must be the dominant ones. Terms which appear later in the series and are 
not used in determining the strength of the singularities, must be small. For 
the models considered^ this has been the case. 
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There is one distinct advantage to this truncated Fourier series technique. 
The force or moment on the plate due to the interaction can be calculated from 
the formula, 


F or M = 


it 

q 

V 0 


r ( 1+n )c cos(n0)drd9 

0° p 


where the force is obtained if n = 0 and the moment results if n - 1. If the 
pressure coefficient is written in the form of equation (1), the integral 
becomes 


or 


-oo i* Ztt 

' k 

I cos(n9) 

y c m cos ( m0 ) 

'o •'o 

m=0 


drd0 


k r oo r 2tt 

j, 



y J ^°o r ^ +n ^ C m cos ( n0 ) cos(m9)drd0 


0 "0 


Due to the orthogonality of the cosine, this now becomes 


F 



q c rdrd9 

00 0 


and 


M 



cos 0 r drd0 


Thus, the total contribution to the force comes from the first term of the 
Fourier series (truncated or complete)^ and the total contribution to the 
moment comes from the second. 

Matching the first couple of terms of a Fourier series takes on more 
significance in the light of these results. If models could be found that would 


209 


have very close agreement in these terms for all values of r, then the two 
quantities of specific interest could be predicted quite well. 

Although comparison of forces and moments between the theoretical 
models and experimental data will not be made, an approximate indication of 
how the forces compare can be obtained from the following. Consider equa- 
tion (1) at 0 = tt/2, 


C (it/ 2) 

r 


o 


C 2 + C 4 + 


If 


c » c_ 

0 2 


» 


which is true in at least some cases, then 


C p (w/2) = c Q 

Thus, since the force is equal to the area under the c Q (r) curve, weighted by 
r, the force will also be equal to the area under the Cp(Tr/2) curve, again 
weighted by r. In this case if the model predicts accurate pressure coeffi- 
cients along the ray 0 = tt/2, then the force predicted will also be accurate. 
Therefore by just comparing pressures at 0 = tt/2, a good indication of the 
quality of force prediction of the model can be obtained. 

The best way to illustrate the above ideas is through an example. The 
complex potential for the doublet model is 


a 2 

w(z) = a Q z + a 1 log (z) + — 


( 2 ) 


where z is normalized by dj. Usually ag is chosen equal to and the pres- 
sure coefficient defined 


C 

P 


P - 


P 


oo 


0/2) pu£ 
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is given by 


C = 
P 


1 / dw\ / dw\ 

y2 \ dz / \dz / 


However, this leaves only two arbitrary constants in the analytical expres- 
sion for the pressure coefficient which can be matched to only the first two 
terms of the Fourier series representation. Since three terms in the series 
are considered necessary to approximate the wake-like distribution for 
0 £ 30°, a third constant is required in the analytical expression. To provide 
this third constant, a^ is chosen equal to U^/U,*,, and the pressure coefficient 
is defined as 


C 

P 


P' 


(1/2)PU^ 


where the primed quantities are for the two-dimensional model and are not 
necessarily equal to the unprimed, or actual quantities. With these defini- 
tions, the pressure coefficient becomes 


C = 
P 


2 /dw\ /dw\ 
a 0 ‘ \dz^ / Vdz/ 


(3) 


These relations can be obtained by considering the static pressure at infinity 
in the model to equal the actual static pressure, while the total pressure (and 
therefore free-stream velocity) of the actual flow differs from that of the 
model. These manipulations leave a^ freeto be used in the pressure matching. 

Substituting equation (2) into equation (3) gives 


C 

P 



zz 


2_2 
z z 


a o a i 


z + z 
zz 


+ a 0 a 2 


2 -2 

z+z z + z 

2_2 + a l a 2 2_2 

z z z z 


(4) 


Equation (4) is now to be expressed in the form of equation (1). This can be 
done directly for the present example by the substitution 

i0 — -i0 

z = re , z = re 
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or for more complex models by the integral definitions of the coefficients 
c n (r). 


if* 


c o (r) = ¥ I V r<e)de 


c (r) 
n 


sf 

fk 


C^(r, 0)cos(n0)d0 


After some algebraic reduction the constants for the present example are 
found to be 


c o (r) 


Cj (r) 

c 2 ( r ) 



(5) 


These expressions for the coefficients c n are to be set equal to experimentally 
determined values at a fixed value of r. The location chosen was the circle 
r = 1, a choice which is somewhat arbitrary. For a velocity ratio (U^/Uj) of 
0. 4, the values of the coefficients c n are 


c 0 (l) - -0.619 
CjU) = -0.843 
c 2 (l) = +0.287 


Numerical solution of equations (5) then gives 
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a Q = 0. 745 ^ 

a 1 = 0. 762 J 

a 2 = 0. 192 j 


( 6 ) 


These results indicate that for the proper behavior, the source is indeed a 
source and not a sink. Further, the blockage effect of the doublet, which 
opposes the free stream, is considerably smaller than the effect of the source. 

Expressing equation (4) in the form of equation (1) gives a summation 
which contains just three terms. Thus, the constants given by equations (6) 
will reproduce the three- term cosine series exactly at r = 1. An idea of the 
accuracy possible through use of this doublet model can be obtained from 
figure 3 in which results from a three-term series are compared to the data 
at r = 1. It can be seen that the agreement of the doublet model with the data 
is quite good for most values of 0. The real test of the model, however, is 
the agreement at other values of r. Comparison of the model with Vogler's 
dataforfive rays, 0= const, is shown infigures 8 through 12. Examination of 
these figures shows consider able disagreement around the upstream and down- 
stream rays and fairly good agreement near the ray 0= tt/2. 

The accuracy of other results obtained from this model at different 
velocity ratios is similar to that described above. Thus, it can be concluded 
that this choice of singularities will not yield pressure distributions suffi- 
ciently accurate to be of practical use. 

The other phenomenological model employed in this study is shown in 
figure 7, as mentioned previously. In addition to the source (or sink), dou- 
blet, and free stream, there are two vortices of equal strength (but opposite 
sign) located downstream of the origin. These are included to represent the 
vortex motion that actually exists in the flow. It was at first expected that by 
taking the reverse flow downstream of the jet into account, increased accu- 
racy would result. 

The complex potential for this flow is 


w (z) 


z 4 a ^ 


log(z) + — + ia 3 log 



(7) 
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Here z is, again, normalized by dj, and w is normalized by U M , Since the 
free- stream velocity in this model is U^, the usual definition for the pressure 
coefficient is chosen in the form 


C = 
P 



1 


/ dw\ / dw\ 

\dz/ \dz / 


Calculating C from the above expression results in 
P 


= 4. & Jl _ b a 2 \ U + z \ *1 

1 -2 " a l \ " zz) \ zz / zz 


2 2 
2 


2.2 

z z 


2-2 2 . 2 
zz + a z z - a, z 

+ a (z - z ) ( 9 + - 

3 o o I 2_2 , w _ , z 


z + a j z - a^ 


z z (z - z )(z - z ) 
o o 


z z / zz \ / zz \ 

•'h h 


+ a 


0 (z - z ) 

2 o o 

3 z z 

o o 


(z - z ) (z - z 
o v o' 


(■-?)("?)) 


( 8 ) 


This expression contains five unknowns, aj_, a^, a 3 , r Q , and 0 O , where the 
last two are the coordinates of the location of the vortices. These unknowns 
can be determined in several ways. One is to express the pressure coeffi- 
cient in a five-term cosine series and equate the coefficients to a five-term 
series of the data. Another is to write a three-term series and equate the 
coefficients at two values of r (i. e. , equate three coefficients at one loca- 
tion and two at the other). A third option is to use physical considerations 
to determine two unknowns, and use the three-term series to determine the 
other three. This third option was chosen as the one to use, although it is 
obviously an arbitrary choice. The physical consideration is the so-called 
vortex equilibrium criterion; namely, the velocity at the location of each 
vortex due to all of the other singularities is zero. This criterion insures 
that there is no force tending to convect the vortices downstream (ref. 10). 

It yields the two relations, 


cos9 a, 

o _ i_ 

r “ 2a 0 

o 2 
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and 



After some algebraic reduction, equation (8) can be written in the form of 
equation (1) and the coefficients c n set equal to those determined from the 
data. The expressions are rather lengthy and will not be presented here. 
Solutions of these three equations along with the two stability equations have 
been obtained numerically for a velocity ratio of 0. 4. Since the equations 
are nonlinear, several extraneous roots appear. However, there is little 
doubt as to which ones are physically realistic. The results are 


a 1 = -0.121 

a 2 = -0.048 

a^ = 0. 668 

0 = 36.3° 

o 

r = 0. 64 
o 

These again show that the doublet is considerably weaker than the source but 
in this case, the source is really a sink. Physically, this result seems more 
realistic than the result of the doublet model since the sink represents 
entrainment. Also, the vortices have the sense of rotation expected physically. 

The value of the pressure coefficient determined from these results as a 
function of r for various values of 0, compared with the data of reference 9, is 
again shown in figures 8 through 12. Evidently, the agreement is unsatis- 
factory. For much of the surface of the plate, the agreement is worse than 
that obtained with the simpler doublet model without vortices. 

The failure of the doublet-vortex model to correctly predict pressure 
distributions on the leeward side points to a fundamental difficulty of all 
inviscid models which attempt to simulate the observed velocities near the 
plate surface. It is known from experiment (ref. 11) that immediately 
behind the jet the fluid velocity is inward toward the orifice. (The presence 
of vortices in the doublet-vortex model was supposed to account for this fact. ) 
Considering figure 7, it is evident that because of symmetry and the condition 
of no flow through the plate, the velocity vector along the 0=0 ray must be 
aligned with that ray. Consequently, if far downstream perturbations are to 
decay and the velocity along the 0=0 ray is to become equal to U^, the flow 
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must reverse direction and have a stagnation point in between. If no account 
is taken of viscous dissipation, the pressure at this stagnation point will 
recover to free- stream stagnation pressure, and the value of the pressure 
coefficient will be unity. Apparently, this difficulty will be encountered with 
any inviscid model that attempts to reproduce the inward velocity observed 
experimentally, no matter how complicated. 


Pressure Models 

For these models the point of view is taken that since the pressure dis- 
tribution is known qualitatively, a flow can be constructed that will yield 
approximately this pressure distribution. At this point, free parameters (if 
there are any) are adjusted to obtain the best agreement with the data. 

With this viewpoint, it is clear that the results obtained from the model 
should have a high-pressure region on the windward side of the jet and a low- 
pressure region on the leeward side. These both asymptotically decay to 
free-stream pressure at infinity. The nature of potential flow is such that to 
obtain this sort of pressure distribution, the velocity of the fluid due to what- 
ever represents the jet must be in the same direction as the free stream on 
the leeward side of the jet and opposing the free stream on the windward side. 
Further, the singularities making up the disturbance to the free stream must' 
be restricted to the neighborhood of the origin. 

The simplest model which will meet those criteria is a source in a free 
stream. A source, doublet and free-stream model is more general and was 
tried above. However, by considering the 0 = 90° ray and the fact that the 
sign of the doublet is positive it can be seen that agreement will improve if 
the doublet was removed. Another model was tried which consisted of the 
superposition of the axisymmetric flow out of an orifice in a flat plate and a 
uniform free stream. Both the source model and this orifice model will be 
described in this section. 

The complex potential for the source model is well known. It is 


w(z) = a Q z + ajlog z 


( 9 ) 


Calculating ag and a^ by using the method described above for the doublet 
model yields the pressure coefficient 


C 

P 


2 

a, , - 

1 z + z 

— =r - a a. — = — 
zz 0 1 zz 


( 10 ) 
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which can be written as 


C 

P 


2 

a l 

~2 ' 


2a 0 a i 


COS0 

r 


(ID 


The coefficients of the cos(n9) terms are now equated to the c n determined 
from experimental data. Since these Fourier coefficients are independent 
of each other, c 0 and c^ for a velocity ratio of 0. 4 should have the same value 
as given above. However, for these last two models an adjustment has been 
made in the data of Vogler. At large distances from the orifice, the pressure 
coefficient should approach zero but the data do not show this. To correct 
this, the simplest approach was taken. A constant was added to the value of 
the pressure coefficient, Cp, on each ray so that it approached zero as r 
grew large. Of course, a different constant was used for each value of 0. 

With this, the values of cq and c 


C 1 " 

where the matching is again at r = 1. 


become 

-0. 725 
-0. 823 

Solving for ag and aj yields 


a Q = 0.483 
a, = 0.851 


Comparison of this model with experimental data is shown in figures 13 
through 17. Like figures 8 through 12^these are plots of the pressure coeffi- 
cient as a function of the radial distance from the center of the nozzle for 
several values of 0. As can be seen, this model gives quite good agreement 
in the vicinity of 0 = 90°, but the agreement around 0=0° and 0 = 180° is not 
as good. This is typical of the results obtained for other velocity ratios. It 
can therefore be expected that while this model will not predict pressure dis- 
tributions accurately for all 0, it will be good for 0 near 90 and will give 
good results for forces on plates. 

The discrepancy between the data point at r = 1 in figures 10 and 15 is due 
to difficulties in reading Vogler' s data at these values for r and 9. 


Curves of the values of 8 q and aq for this model as functions of the veloc- 
ity ratio are shown in figure 18 . 
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The final model to be discussed has been mentioned before. It is 
composed of the superposition of two simpler flows, but this time one is axi- 
symmetric and the other is added to it in such a way that the result is fully 
three dimensional. The axisymmetric portion is the flow out of a circular 
orifice in a flat plate. Reference 12 has a discussion of this solution. The 
other flow is a uniform free stream aligned so that its velocity vector is tan- 
gent to the plate. Carrying out the analysis and obtaining an expression for 
the pressure coefficient on the plate gives 


C - Uj /Up cos6 x (U./UJ 2 

p " ’ T (¥ [(if - >] 


where the pressure coefficient retains its usual definition and k is a function 
of the velocity ratio, which is used to adjust the accuracy of the results. A 
curve of k as a function of the velocity ratio is shown in figure 18, along 
with the dependence of a q and aj on Uco/Uj. 


Comparison of the results of equation (12) with data is shown in fig- 
ures 13 through 17* It can be seen that there is good agreement except near 
0 = 90°* Since the model is inaccurate near 9 = 90°, it cannot be expected 
to give good results for forces (the singularity in the model is excluded 
from the integration), but its accuracy along the 0=0 and l 80 ° rays indi- 
cates that the moments derived from it should agree closely with experiment. 


CONCLUSIONS 

A summary of the conclusions which have been reached in this paper are 

given below 

1. By using a truncated Fourier series, a relatively simple empirical fit 
of the data is possible. While it is clear that an exact fit, in the mean, 
is always possible using an infinite series for a fixed r, it is shown 
that good results are possible by using only two or three terms. 

2. Reasonable results for pressure distributions can be obtained using invis- 
cid models. However, good agreement is usually restricted to a partic- 
ular range of 0. For instance, the source model gives good results 
around 0 = tt/2j while the orifice model is best near 0=0 and 0 = tt. 

3. As a corollary to conclusion 2 , it is possible to obtain forces accurately 
from an inviscid modeLand it is possible to obtain moments accurately 
from an inviscid model. At the present time, it does not seem to be 
possible to obtain both from the same model. 
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4. Adding greater complexity to the phenomenological models does not 
appear to improve their results. This is due to the stagnation point 
occurring along the axis. 

5. Effects of viscosity must be included in any analysis which is intended 
to predict the flow properly, especially downstream of the jet. 
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Figure 3-“ Comparison of two- and three- term series with 
data from reference 9 at r = 1. 



Figure k.- Comparison of two- and three-term series with 
data from reference 9 at r = 2. 
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Figure 7. - Sketch of doublet —vortex model. 
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8.- Results of doublet -vortex and doublet models compared 
with data from reference 9 . 9=0°. 
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Figure $).- Results of doublet-vortex and doublet models compared 
with data from reference 9» 0 = 60°. 
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Figure 10.- Results of doublet-vortex and doublet models compared 
with data from reference 9. 0 = 90°. 
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Figure 11.- Results of doublet -vortex and doublet models compared 
with data from reference 9* 9 = 150°. 
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Figure 12.- Results of doublet-vortex and doublet models compared 
with data from reference 9* 0 = l80°. 
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Figure 13* - Results of source and orifice models compared 
with data from reference 9* 9=0°. 
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Figure 14.- Results of source and orifice models compared 
with data from reference 9* 9 = 60°. 
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THE USE OF MATCHED ASYMPTOTIC EXPANSIONS AS AN APPROACH 
TO THE PROBLEM OF THE JET IN A CROSSFLOW 
By Jack E. Werner 
New York University 

SUMMARY 

The method of matched asymptotic expansions is proposed as an ap- 
proach to the problem of a circular jet issuing from a plane surface 
into a crossflow. "Inner" and "outer" expansions are written in powers 
of the parameter r 0 /R. Introduction of these expansions into Laplace's 
eq. for the velocity potential results in a series of quasi two-di- 
mensional equations for the terms of the inner expansion. The matching 
principle is illustrated and the manner in which different models of 
the flow may be introduced into the theory and compared is discussed. 

SYMBOLS 

r characteristic "diameter" of jet cross section 

o 

R radius of curvature of jet trajectory 

U undisturbed upstream flow velocity 

coordinate in downstream direction 
coordinate perpendicular to x x and jet axis 
coordinate perpendicular to plane of jet orifice 
velocity potential 

outer expansion for velocity potential 
inner expansion for velocity potential 
Q complex variable x + iy 

Q complex variable X + iY 

Re|f(£)j- real part of f (Q ) 


z i 

cp 

0 

cp 

1 

CP 
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INTRODUCTION 


One of the important objectives of an analytical treatment of the 
jet in a crossflow is the determination of the induced flow particularly 
near the aerodynamic surfaces of jet VTOL craft. Unfortunately the 
phenomenon is too complex to admit of a complete theoretical treatment. 
The effects of viscosity and turbulence in the jet present formidable 
enough obstacles but to this must be added the difficulties of solving 
the free-boundary problem of the three-dimensional flow around the jet 
to arrive at the induced flow field. Among the initial efforts in this 
field were investigations of the jet itself and its trajectory, refs. 

[1] to [4]. More recently the results of these investigations have 
been used to infer the shape of the trajectory at the outset and thus 
reduce the induced flow problem to a fixed rather than free-boundary 
problem. Such a treatment has been carried out by Wooler, ref. [5] and 
by Rubbert and Saaris, ref. [6], to evaluate the flow field near lifting 
surfaces. 

Presently the state of development is such that we must accept 
admittedly crude models of the physical situation in order to make 
analytical headway. This being the case it is still worthwhile to 
propose and compare different analytical models in an effort to arrive 
at one or a combination of useful ones. It is the purpose of this 
paper to propose a theoretical framework in which various analytical 
models may be explored and compared. Basically the method centers 
about the use of matched asymptotic expansions, Van Dyke ref. [7], 

FORMULATION OF THE PROBLEM 

We consider the problem of a jet issuing into a crossflow from an 
orifice in a flat plate. Outside the jet we seek solutions to the 
inviscid flow governed by the Laplace eq. for the potential. Far from 
the jet, in the so-called "outer 11 region, we expect the flow to be un- 
disturbed away from the wake. On the jet itself the flow is to be 
tangential to the jet surface and for small values of the parameter 
r 0 /R we may expect this flow to be almost two dimensional but not com- 
pletely so. 

To accomodate these last considerations the "inner" and "outer" 
expansions of the potential are introduced in the form: 


CO 

cp°(x,y, z) = Y 

O 



(x, y, z) 


cp 1 (X, Y, Z) =Y 

O 



(X, Y, Z) 
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where x, y, z are nondimens ional outer variables 
x=x 1 /R , y = Yi/R » z=z 1 /R 
and X, Y, Z are the dimensionless inner variables 


X = 


x 

r c /R 


Y = 


y _ 

r 0 /R ’ 


Z = z 


co° is a solution to the Laplace eq. satisfying the boundary conditions 
far from the jet. cp 1 is a solution to the Laplace eq. satisfying 
boundary conditions on the jet surface. Each of these solutions requires 
another condition to determine them uniquely. The remaining condition is 
the requirement that the two expansions match asymptotically over a 
region intermediate to the inner and outer regions. By this we mean that 
over the matching region the Nth partial sums of cp° and cp 1 exhibit the 
same functional behavior and agree within an order of (tq/r)^. Thus the 
condition supplying the intermediate boundary condition is for each value 

of N 


lim 

(r o /R)~*0 


N m N m 

[1(f) ^ -[1(f) 4 


(r 0 /R) 


N 


= 0 


Where [ ] represents the behavior of the inner solution near the 

outer bound of the matching region, i.e. , as X, Y -* co with x, y finite this 
limiting process is equivalent to letting Tq/R approach zero with x,y 
finite. Similarly [ '] i represents the behavior of the outer solution 

near the inner bound of the matching region and is found by allowing r 0 /R 
to approach zero with X, Y finite, i.e., letting x,y tend to zero. 

EVALUATION OF cp° , cp^ 

The governing eqs. for the functions cp° and c0* are obtained by 
substitution of the expansions into the Laplace eq. and setting the coef- 
ficient of powers of (r 0 /R) equal to zero. Thus the complete eqs. 
governing cp° and cp 1 are 


& 


+ f^ + 

. c^cp 1 , 
+ FT2 + 



= 0 
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Substituting the inner expansion into the second of these yields 


- + = o 


C^cp 1 c^cp?” 


S 2 cp 1 ^cp 1 d 2 cp 1 

m _ __m-s 

dF + dF dF 


m 2: 2 


The terms of the outer solution all obey 

O N8 o \s o 

ycp m & cp y® B 

+ W~ ‘ 0 


m > 0 


Thus cp^ is a solution of the Laplace eq. satisfying the outer boundary 
condition of uniform flow at infinity. A solution for cp° will therefore 
be ° 


Cp Q = Ux 


From the matching condition with N = 0 


(r o /R)-0 


ML - K]J 


Where cp 


H r r -i 

lim [Ux] = lim U X = 0 

-*1 f-r fv\-+Ci /r>V .n L R 


J i (r o /R)-0 
X finite 


(r o /R)-0 
X finite 


Therefore we have [ cp * ) - 0 which together with the governing eq. for m*" 

L oJq M o 

and the inner boundary condition yields 


CD o = 0 


This result, although very simple, serves to illustrate the rudiments of 
the matching technique. Turning to the next order cp^ satisfies the two 
dimensional Laplace eq. together with the condition that its derivative 
normal to the jet surface is zero to within terms of order (r Q /R). A 
general solution satisfying these conditions will be assumed in the form 
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! 


KI * c ? + lh} 


This form assumes a general two-dimensional symmetrical flow about a 
closed streamline enclosing the origin Q = 0. The closed streamline 
may be completely on the jet surface, i. e. : fully attached flow, or may 
follow the upstream portion of the jet surface but then separate from 
it to enclose a wake region and then close at some point downstream. 

It is at this point that we may introduce any one of a number of possi- 
ble models for the flow around the jet. 

Applying the matching condition with N = 0 we arrive at the con- 
ditions 

X C = U 

n C = 0 n s 2 



Thus we shall have cp° = 0 and 

00 

O ^ 

in which the constants C are determined by the requirement that co 1 
represent the potential for a uniform flow about the closed streamline 
(tangential to jet and wake surfaces) at the station z. In other words, 
the third dimension z enters by means of the coefficients determined 
by the boundary condition at the station z. 

Continuing the process to N = 2 yields the result 



% ' Re { l Vc"; 
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The interpretation placed on the latter result is that to within the 
order of (r Q /R ) 2 the outer solution consists of a uniform flow plus a 
distribution of doublets along the jet axis. The origin of the doublet 
flow is found in the "blockage" of the jet and the presence of sheets 
of vorticity occurring symmetrically over the jet and wake surfaces. 
Since the strength of these vortex sheets are determined by the jet 
cross section at a given station z along the jet axis, the vorticity 
will vary with z. To maintain the balance of vorticity according to the 
laws of Helmholtz we must have a trailing vortex system breaking away 
symmetrically from the sides of the jet and wake surfaces. Far from the 
jet these trailing vortices appear as a distribution of vortex doublets 
inducing a downwash behind the jet. 

CONCLUDING REMARKS 

The above considerations represent a rough sketch of the more 
important features of the proposed method as they are discernible at 
this stage. Principal among these is the quasi two-dimensional nature 
of the equations fqr cp*, cp 1 and at least the homogeneous parts of the 
subsequent terms cp*. This 2 makes available to us the growing body of 
knowledge regarding two-dimensional wake flows. The method also con- 
tains a built-in system of "accounting" in the matter of obtaining 
solutions which are valid to within a given order of magnitude. Thus, 
in obtaining say the first two terms of cp 1 it is only necessary to 
satisfy an approximate boundary condition to within an order of (r Q /R) 
rather than inconsistently posing an exact boundary condition for an 
approximate formulation of the solution. 
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MASS ENTRAINMENT OF A CIRCULAR JET IN A CROSS FLOW 


By R. L. Fearn 

University of Florida 
Gainesville, Florida 


SUMMARY 


By use of conservation of mass and momentum, an equation is derived 
relating the mass entrainment of a round turbulent jet exhausting into a cross 
flow to the equation for the jet axis. If an empirical equation for the jet 
axis is used, one finds that the mass entrainment approximates a linear depend- 
ence on distance along the jet axis. If, on the other hand, a linear entrain- 
ment of undetermined slope is assumed, a differential equation is obtained for 
the jet axis which can be solved analytically. Several members of the resulting 
family of one-parameter curves are fitted to empirically determined jet axes to 
obtain the entrainment as a function of the ratio of free-stream velocity to jet 
velocity. 


INTRODUCTION 


An approximate technique to determine numerically the aerodynamic charac- 
teristics of a wing in subsonic flow is to use lifting- surface theory. This 
method is a generalization of Prandtl's lifting-line theory in which the wing 
is represented by a superposition of singular solutions to Laplace's equation. 
Boundary conditions that the flow be tangent to the surface of the wing and that 
the Kutta condition be satisfied are imposed at appropriate points. This method 
has also been used in an attempt to calculate the aerodynamic characteristics of 
the fan-in-wing configuration which is one design approach to VTOL aircraft 
(ref. 1). The slipstream from the wing fan constitutes a round turbulent jet 
exhausting into a cross flow. Since this is a boundary- layer phenomenon, one 
would expect the singularity method to be applicable in the region exterior to 
the wing and jet. However, knowledge of the jet in a cross flow is required to 
establish the boundary conditions for the potential flow problem. The lack of 
such knowledge is illustrated by the work of Rubbert and associates who found 
that the only usable information available to them about the jet in a cross flow 
was an empirical equation for the jet axis. In the absence of any other informa- 
tion, they used a constant-cross-section jet with no entrainment. (See ref. 1.) 

Despite this lack of information, considerable effort has been expended to 
describe this phenomenon both experimentally and theoretically. Experimental 
studies have established the location of the jet center line for various condi- 
tions at the jet orifice and have provided some information about velocity and 
pressure distributions. (See refs. 2 and J. ) Because of the complexity of the 
theoretical problem, however, little more has been attempted than to describe 
the location of the jet center line. One theoretical approach has been to 


assume that the interaction of the jet with the free stream could be represented 
by an appropriate drag coefficient for the jet (refs. 4 to 6) which presupposes 
that a pressure distribution is a mechanism for jet deflection. This is possibly 
significant in the vicinity of the jet orifice, but it is probably negligible 
compared with the momentum mixing of the jet with the free stream for distances 
along the jet axis greater than a few jet diameters. 


SYMBOLS 


a 


Cl 

E 


e s 

A A 

i,k 


M 

Mz 




m 


Q 

Qo 

A %: 


s 


x,z; 

x o> z 


X* 


parameter, M z /mV oo , feet 


constant, a.\jl + (x Q ^ - x Q , dimensionless 
entrainment coefficient, dimensionless 

Q 0 ds 

unit vector tangent to jet axis 
Cartesian unit vectors 


momentum flux across jet cross section, slug- feet /second^ 

component of jet momentum flux perpendicular to free-stream velocity, 
slug-feet /second^ 

momentum flux entrained across increment of jet boundary, 
slug- feet/ second2 

slope of curve of mass entrainment as a function of distance along 
jet axis, constant, slugs/foot-second 

mass flux across jet cross section, slugs/second 

mass flow across jet orifice 

mass flux entrained across increment of jet boundary, slugs/second 
distance along jet axis, measured in jet diameters, dimensionless 
initial jet velocity, feet/second 
free-stream velocity, feet /second 

Cartesian coordinates of jet axis, measured in jet diameters, 
dimensionless 

derivative of x with respect to z, dimensionless 
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z' ,z" first and second derivatives of z with respect to x, dimensionless 
a angle between jet axis and free-stream velocity, degrees 

Arrow over symbol indicates vector. 


MASS ENTRAINMENT AND THE EQUATION FOR THE JET AXIS 


The purpose of this paper is to investigate the possibility of obtaining 
some information about the jet axis and entrainment of a circular turbulent jet 
exhausting into a cross flow by considering conservation of mass and momentum. 
Simply stated, the approach is to assume that the mechanism for deflecting the 
jet is an inelastic collision of the entrained mass with the mass of the jet,-^ 
that is, the pressure field as a mechanism for deflecting the jet is neglected. 
Additionally, the mass entrained by the jet is assumed to have a velocity 
approximately equal to the free-stream velocity V^. 

To formulate this assumption quantitatively, consider the conservation of 
mass and momentum for a small control volume defined by two cross sections of 
the jet, located at s and s + As from the jet orifice along the jet center 
line, and by the jet boundary. The conservation equations for the control 
volume may be written as 


Q(s+As) = Q(s) + AQgts) 

(1) 

m(s+as) = m(s) + m E ( s) 

(2) 


where Q(s) and M(s) are the mass and momentum flux across the jet cross sec- 
tion at s, and AQe(s) and AMjj;(s) are the mass and momentum flux entrained 
across the jet boundary between the cross section at s and s + As. Assume 
that the momentum entrained across the jet boundary is the entrained mass with a 
velocity equal to the free-stream velocity V^, that is. 


AMg(s) = AQg(s) 


(3) 


The change in the momentum of the jet with respect to distance along the jet 
axis may thus be written 


dM _ dQ =» 
ds ds 00 


00 


■'•This approach was first brought to the author's attention by Dr. Gerhard 
Braun in December 1967. 


241 



If dM/ds is related to the equation for the jet axis, then equation (k) 
can be used to obtain a relation between mass entrainment and the equation for 
the jet axis. Define the jet axis as the curve passing through the center of 
the jet orifice with the direction at every point on the curve t an gent to the 
momentum flux at that point. If the momentum flux is written as 


M(s) = M(s) e s 


and the unit vector e s is written in terms of the Cartesian unit vectors, then 


dM 

ds 


1 3M z 
z ' ds 


z"M„ 


. ,21 ^ N 1 / 2 

z ' \1 + z ’ ) 



(5) 


where i and k are the unit vectors parallel and normal to the free stream, 
respectively, z(x) defines the equation for the jet axis, and M z is the com- 
ponent of the jet momentum flux perpendicular to the free- stream velocity. 
Equating equations (k) and ( 5 ) gives the desired relation 


^ V = - 
ds 00 


z"Mz 


,2 


(l + z' 2 ) 


lj2 


(6) 


where M z is constant. 

Using an empirical equation for the jet axis to calculate the mass entrain- 
ment as a function of s, one finds that the dependence is approximately linear 
except in the vicinity of the jet orifice (fig. 1). 


A SEMIEMPIRICAL EQUATION FOR THE JET AXIS 


If the mass entrainment is assumed to vary linearly with distance along the 
jet axis in the region of established flow, that is, 

dQ/ds = m 


where m is a constant, then the jet axis is described by the differential 
equation 


mV = 


X 


z ,2 (l + z’ 2 ) 


1/2 


(7) 


242 



Letting a = M^mY^) and integrating equation (7) results in a semiempirical 
equation for the jet axis 

[(x + Ci) 2 - a 2 ] ^ + x + C! 


z = z Q + a In 


Qxq + C1) 2 - a 2 J + x Q + Ci 


( 8 ) 


where Cp = a|l + (x q ') 2 j - x Q and the parameter a must be determined 

experimentally. The position (x 0 ,z 0 ) and the slope x Q ' of the jet axis at one 
point in the region of established flow must be specified as initial conditions. 


For a jet exhausting normally into a cross flow, 


M z = 

where Qq is the mass flux across the jet orifice. One may define an entrain- 
ment coefficient as 


1 _ dQ 
Qo ds 


which in terms of the parameter a is 



aV 


The parameter a, or equivalently E, may be varied to obtain a best fit to 
experimentally determined jet axes for various velocity ratios X»/^j* 

The entrainment coefficient for a number of velocity ratios was obtained 
by fitting equation (8) to the following empirical equation for the jet axis 
developed by Ivanov (ref. 4 , p. 5^) : 


(9) 



where a is the initial angle that the jet axis makes with the free- stream 
velocity. Margason confirms that this equation provides a good fit to experi- 
mental results (ref. j). It was found that a good fit between the semiempirical 
equation (8) and the empirical equation (9) could be obtained by assuming that 
the jet does not deflect from the Z-axis until it reaches a particular distance 
from the jet orifice. The point at which deflection begins is referred to as 
the effective source of the jet, and the distance from the jet orifice to the 
effective source varies with velocity ratio. This probably marks the point 
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where the potential core has disappeared and where the region of established 
flow begins. Under these assumptions, the semiempirical equation simplifies to 


( 10 ) 


where z G locates the effective source , and R is the velocity ratio VyVh. 

The entrainment coefficient E as a function of velocity ratio is 

shown in figure 2. This entrainment coefficient is relatively insensitive to 
the matching point, that is, to whether equation (lO) or equation (8) is used. 

The entrainment for the submerged jet (ref. 8) has been included for comparison. 
The location of the effective source z Q , which is obtained in addition to the 

entrainment coefficient when equation (lO) is fitted to the empirical curve, is 
shown for several velocity ratios in figure 3* An example of the fit between 
the semiempirical and empirical equations for the jet axis is presented in 
figure 4. 

From consideration of the integral equations for the conservation of mass 
and momentum, it becomes apparent that equation (4) is an oversimplification. By 
assuming that the velocity on the jet boundary does not differ appreciably from 
the free- stream velocity, it can be shown that the entrainment presented in fig- 
ure 2 may be too large by as much as a factor of 2. 



CONCLUDING REMARKS 


It has been shown that the observed jet axis can be accounted for by con- 
sidering momentum mixing to be the only mechanism for jet deflection. Recourse 
has not been made to a drag coefficient which in the region of established flow 
has the disturbing connotation of flow over a solid body. The values predicted 
for the mass entrainment necessary to account for the observed jet deflection 
are physically reasonable and should be considered as upper limits. 
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